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1.0  INTRODUCTION 

The  work  reported  herein  was  performed  by  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Materiel  Command  (AFMC),  under  Program  Element  65807F.  The  Air 
Force  Program  Managers  were  Capt.  D.  G.  Burgess,  Capt.  H.  Martin,  and  Lt.  P.  Zeman, 
DOT.  The  work  was  performed  by  Calspan  Corporation/ AEDC  Operations,  technical  services 
contractor  for  the  Aerospace  Flight  Dynamics  testing  effort  at  AEDC,  AFMC,  Arnold  Air 
Force  Base,  TN.  The  work  was  performed  in  the  Aerospace  Systems  Facility  (ASF)  and  the 
Technology  and  Development  Facility  (TDF)  under  AEDC  Project  Number  DD01,  Job 
Number  0115.  The  work  was  conducted  during  the  period  between  1  October  1990  and  16 
October  1992. 

Water  cooling  is  used  extensively  at  AEDC  to  prevent  the  failure  of  high-enthalpy  arc 
heater  components  at  elevated  heat  loads  where  cooling  by  radiation,  conduction,  and  natural 
convection  from  air  circulation  are  ineffective.  Other  cooling  techniques  such  as  the  use  of 
liquid  metals,  transpiration  or  film  cooling  pose  their  own  set  of  disadvantages  or  limitations 
(Ref.  1)  but  are  certainly  possibilities  for  providing  cooling  at  much  more  severe  heating 
conditions  where  water  cooling  alone  may  prove  to  be  inadequate.  The  current  interest  is 
in  exploiting  backside  water  cooling  to  the  limit  of  its  capability,  especially  for  arc  heater 
nozzles.  Such  an  application  requires  an  understanding  of  the  parameters  that  affect  the  cooling 
processes  and  the  determination  of  the  limiting  point*  where  the  surface  fails  due  to 
inadequate  cooling. 

The  purpose  of  this  report  is  threefold:  (1)  to  review  the  different  processes  encountered 
in  backside  water  cooling  and  identify  those  processes  that  are  desirable  for  effective  cooling 
of  high-enthalpy  facility  components;  (2)  to  identify  those  parameters  which  affect  the  cooling 
processes  of  interest  and  the  trend  of  cooling  performance  when  the  parameters  are  individually 
varied;  and  (3)  to  present  a  summary  of  applicable  theoretical,  analytical,  and  experimental 
work  in  surface  cooling  performed  in  the  past  50  years. 

2.0  THE  BASIC  BOILING  CURVE 

The  characteristic  temperature  response  of  a  surface  cooled  by  an  adjacent  fluid  as  heat 
flux  is  increased  was  first  proposed  by  Nukiyama  (Ref.  2)  in  1934.  Figure  1  shows  the  typical 


*  The  limiting  point  is  also  identified  as  the  burnout  point,  boiling  crisis,  departure  from  nucleate  boiling  (DNB),  or  critical 
heat  flux  (CHF)  condition,  among  other  terms  by  various  authors.  Terminology  for  the  limiting  point  will  be  discussed  in 
Section  2.0. 
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shape  of  the  curve,  known  as  the  boiling  curve,  obtained  when  heat  flux  from  a  cooled  surface 
is  plotted  as  a  function  of  the  liquid  superheat  (the  difference  between  the  surface  temperature 
and  the  liquid  saturation  temperature).  At  low  heat  flux  (segment  A-B  in  Fig.  1),  the  surface 
temperature  is  below  the  saturation  temperature  of  water  and  similar  liquids  at  or  above 
atmospheric  pressure.  Cooling  at  low  heat  flux  is  by  pure  convection  heat  transfer  and  is 
fairly  well  understood.  In  this  region,  analytical  solutions  and  experimental  correlations  (for 
more  complex  flow  geometries)  have  been  derived  which  allow  for  accurate  prediction  of 
heat  transfer  in  various  cooling  configurations  (see,  for  example,  Chapters  8-9  of  Lienhard 
(Ref.  3),  Chapter  6  of  White  (Ref.  4),  or  Kays  and  Crawford  (Ref.  5),  Chapters  13-14).  As 
heat  flux  from  the  surface  is  increased  beyond  the  point  where  the  surface  temperature  becomes 
significantly  greater**  than  the  water  saturation  temperature  (point  B  in  Fig.  1),  vapor 
generation,  better  known  as  boiling,  occurs  from  discrete  nucleation  sites  on  the  surface  or 
in  the  bulk  fluid  depending  on  the  bulk  fluid  temperature.  The  region  of  boiling  from  discrete 
sites  is  termed  nucleate  boiling  and  is  shown  as  segment  B-C  in  Fig.  1. 

As  the  upper  limit  of  nucleate  boiling  is  reached  (point  C  in  Fig.  1),  a  typically  undesirable 
transition  to  another  form  of  boiling,  known  as  film  boiling,  begins.  Because  the  boiling 
near  the  limit  of  nucleate  boiling  becomes  so  violent,  numerous  vapor  bubbles  begin  to  combine 
near  the  surface.  Since  the  vapor  has  a  much  lower  thermal  conductivity  than  the  liquid, 
the  patches  of  vapor  effectively  insulate  the  surface  from  the  efficient  heat  transfer  to  the 
liquid  as  experienced  during  nucleate  boiling.  Unlike  a  temperature-driven  system  which  would 
proceed  toward  point  D+  in  Fig.  1,  a  heat  flux-driven  system  such  as  a  high-enthalpy  arc 
facility  would  proceed  directly  toward  point  Eff.  The  processes  involved  when  a  boiling 
system  undergoes  the  path  from  the  limit  of  nucleate  boiling  (point  C)  toward  stable  film 
boiling  (point  E)  are  highly  unstable,  and  failure  of  the  surface  due  to  melting  occurs  in  most 
cases  because  the  melting  temperature  of  the  surface  material  is  lower  than  the  temperature 
associated  with  stable  film  boiling  at  point  E.  The  point  at  which  transition  to  film  boiling 
begins  is  often  termed  the  burnout  point,  even  though,  in  some  cases,  failure  does  not  occur 
(e.g.,  a  temperature-driven  system  or  a  surface  material  with  a  very  high  melt  temperature). 


**  A  certain  degree  of  liquid  superheat  is  required  to  establish  a  flow  of  heat  to  create  vapor  bubbles.  The  superheat  can  become 
quite  large  for  very  smooth  surfaces  or  highly  wettable  fluids  (e.g.,  refrigerants  or  liquid  metals)  which  can  lead  to  explosive 
boiling  (known  as  bumping)  that  can  lead  to  structural  damage  to  hardware. 

t  Transition  boiling  occurs  between  the  limit  of  nucleate  boiling  (point  C)  and  point  D,  which  is  known  as  the  Leidenfrost 
point.  The  path  between  the  points  is  dependent  on  numerous  factors  and  the  heat  flux  becomes  a  multivalued  function  of 
the  liquid  superheat  (e.g.,  see  Witte  and  Lienhard,  Ref.  6).  This  region  is  inaccessible  with  a  heat  flux  controlled  system. 

tt  Point  E  lies  on  a  portion  of  the  boiling  curve  known  as  film  boiling.  Beginning  at  the  Leidenfrost  point  (point  D),  heat 
transfer  across  the  vapor  film  adjacent  to  the  surface  increases  due  to  radiation  from  the  wall,  through  the  vapor,  to  the  bulk  fluid. 
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The  burnout  point  is  known  also  as  the  departure  from  nucleate  boiling  (DNB),  fast 
burnout,  and  burnout  of  the  first  kind  in  systems  which  involve  a  subcooled*  liquid.  In 
saturated  liquid  boiling  systems,  another  type  of  burnout  occurs  (usually  at  lower  heat  flux 
than  subcooled  boiling),  and  the  terms  slow  burnout,  dryout,  or  burnout  of  the  second  kind 
are  used  (see  Ref.  8).  Boiling  crisis  and  critical  heat  flux  condition  are  used  in  a  general  sense 
to  describe  all  burnout  situations.  Because  of  a  lack  of  standard  nomenclature,  hereafter 
the  point  at  which  transition  to  film  boiling  begins  is  designated  the  burnout  point  or  boiling 
crisis,  and  the  peak,  ultimate,  or  maximum  heat  flux  is  termed  the  critical  heat  flux  (CHF), 
except  where  designated  differently  in  the  referenced  material. 

Design  of  the  cooled  components  of  a  high-enthalpy  facility  for  safe  operation  necessitates 
the  cooling  processes  to  lie  to  the  left  of,  and  well  below,  point  C  in  Fig.  1 .  High  heat-transfer 
coefficients,  hence,  efficient  energy  transfer,  can  be  realized  when  operating  in  the  nucleate 
boiling  regime  rather  than  the  pure  convection  regime  at  low-to-moderate  coolant  pressure 
and  velocity.  Any  method  used  to  increase  the  CHF  will  provide  a  factor  of  safety  when 
operating  well  into  the  fully  developed  nucleate  boiling  regime.  The  determination  of  the 
CHF,  along  with  the  shape  of  the  nucleate  boiling  curve  for  a  given  cooling  system,  remains 
one  of  the  most  challenging  areas  of  research  despite  the  considerable  attention  it  has  received 
in  the  last  50  years  (Ref.  9).  Complexity  of  the  boiling  mechanisms,  their  unsteady  nature, 
and  the  small  scale  of  the  heat -transfer  processes  have  hindered  the  development  of  theoretical 
models  and  understanding  of  the  phenomena  from  experimental  studies.  Moreover,  the  cooling 
flow  structure  and,  hence,  the  heat-transfer  processes  vary  depending  on  the  heating  condition, 
fluid  properties,  and  surface  geometry.  For  instance,  boiling  heat  transfer  that  occurs  in  a 
high-enthalpy  facility  (subcooled  flow,  surface  boiling  optimized  for  effective  wall  cooling 
over  short  lengths  at  high  heat  loads)  can  be  quite  different  than  that  which  occurs  in  a  nuclear 
reactor  (saturated  flow  with  possible  slug  or  annular  two-phase  flow,  bulk  boiling  optimized 
for  efficient  energy  transfer  over  great  lengths  at  low-to-moderate  heat  loads). 

The  objective  of  a  large  portion  of  the  experimental  work  undertaken  in  recent  years  in 
the  area  of  boiling  heat  transfer  has  been  to  obtain  an  empirical  or  semiempirical  correlation 
for  a  given  cooling  configuration  and  set  of  test  conditions.  Many  of  the  correlations  have 
limited  application  and  cannot  be  generalized  for  use  in  analyzing  other  cooling  configurations 
and  conditions.  The  experiments  do  provide  a  qualitative  evaluation  of  what  parameters  affect 
the  boiling  processes  of  interest  and  the  trend  in  cooling  performance  when  each  of  the 


*  If  the  cooling  fluid  bulk  temperature  is  well  below  the  saturation  temperature  of  the  fluid,  then  the  fluid  is  known  as  subcooled. 
Boiling  that  occurs  on  a  surface  adjacent  to  such  a  fluid  is  known  as  subcooled  local  or  surface  boiling.  If  the  bulk  temperature 
is  at  or  near  the  saturation  temperature,  the  fluid  is  known  as  saturated.  Bulk  boiling  can  occur  when  the  fluid  is  near  saturation, 
and  undergoes  various  flow  regimes  in  forced  flow  such  as  slug,  annular,  and  drop  flow.  See  Collier  (Ref.  7)  for  further  discussion 
of  saturated  flows. 
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parameters  is  varied.  Optimizing  the  influential  parameters  for  high-enthalpy  facility 
component  cooling  will,  in  turn,  aid  in  identifying  those  boiling  mechanisms  which  play  a 
part  in  the  cooling  process  of  interest.  Once  the  mechanisms  are  identified,  appropriate 
analytical  models  and  experimental  data  may  be  assembled. 

3.0  PARAMETRIC  EFFECTS  ON  THE  BOILING  CURVE 

More  than  20  parameters  that  affect  the  cooling  process  have  been  identified.  Some 
parameters  play  an  important  role  from  pure  convection  to  boiling  crisis  while  others  play 
intermittent  roles.  For  instance,  velocity  has  a  large  effect  on  cooling  performance  in  the 
pure  convection  regime  and  at  boiling  crisis,  but  it  has  virtually  no  effect  in  the  nucleate 
boiling  regime.  Therefore,  the  contribution  of  each  parameter  in  each  of  the  three  regions 
of  interest  (pure  convection,  nucleate  boiling,  boiling  crisis)  will  be  discussed.  The  effects 
of  various  parameters  on  the  microscopic  boiling  scale  (i.e.,  bubble  characteristics)  have  not 
been  studied  in  great  detail  and  are  presented  here  for  the  more  influential  parameters.  The 
cooling  process  is  not  without  parametric  distortion;  that  is,  the  change  in  one  parameter 
causes  a  change  in  another  parameter.  Because  of  the  complexity  of  the  boiling  process,  many 
of  the  parametric  distortions  are  ill-defined  or  unknown;  therefore,  a  discussion  of  these 
distortions  included  here  is  limited.  Several  excellent  reviews  of  various  parametric  effects 
on  boiling  are  available  and  are  drawn  upon  in  the  following  review,  most  notably:  boiling 
in  general — Tong  (Ref.  10),  Westwater  (Refs.  11  and  12),  Rohsenow  (Refs.  13  and  14),  Addoms 
(Ref.  15),  and  Katto  (Ref.  16);  convective  boiling — Collier  (Ref.  7);  subcooled  flow  boiling — 
Boyd  (Refs.  1  and  17),  Bergles  (Ref.  8),  Macbeth  (Ref.  18),  Hewitt  (Ref.  19),  Fiori  and  Bergles 
(Ref.  20),  and  Hughes  (Ref.  21). 

3.1  IMPORTANT  PARAMETERS  AND  THEIR  EFFECTS  ON  HEAT  TRANSFER 
Velocity 

Velocity  (or  mass  velocity,  G,  when  combined  with  the  fluid  density),  as  mentioned 
previously,  plays  an  intermittent  role  in  the  cooling  process.  Figure  2  (from  Kreith  and 
Summerfield,  Ref.  22)  illustrates  the  significant  effect  of  fluid  (butanol)  velocity  as  heat  flux 
from  an  adjacent  surface  increases  in  pure  convection  heat  transfer.  However,  as  transition 
to  boiling  occurs,  the  effect  of  velocity  diminishes.  Clark  and  Rohsenow  (Ref.  23)  and  Tong 
(Ref.  10),  among  others,  have  noted  that  velocity  has  little  or  no  effect  on  heat  transfer  in 
vigorous  (or  fully  developed)  nucleate  boiling.  Analysis  of  boiling  processes  is  divided  into 
two  categories  with  regard  to  velocity.  Boiling  which  occurs  on  a  surface  or  in  a  fluid  initially 
at  rest  where  only  free-con vection  motion  is  present  is  known  as  pool  boiling.  Boiling  which 
occurs  on  a  surface  or  in  a  fluid  where  the  fluid  has  a  bulk  (imposed  by  some  external  source 
such  as  a  pump)  and  free-convection  motion  is  known  as  forced-convection  boiling  or  simply 
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flow  boiling.  The  advantages  of  flow  boiling  over  pool  boiling,  when  present  in  a  cooling 
process,  become  evident  as  boiling  crisis  is  approached.  McAdams,  et  al.  (Ref.  24)  showed 
that  the  CHF  could  be  increased  with  increasing  velocity  (Fig.  3).  Grace  (Ref.  25)  generalized 
the  velocity  effect  on  CHF  for  various  qualities  of  water  liquid/vapor  systems  as  shown  in 
Fig.  4.  An  inversion  in  velocity  effectiveness  occurs  at  higher  vapor  qualities  (i.e.,  saturated 
fluid)  as  illustrated  in  Fig.  4;  however,  as  will  be  discussed  next,  subcooling  (vapor  qualities 
less  than  zero)  provides  the  largest  values  of  CHF,  where  increasing  mass  velocity  increases 
the  CHF  for  a  given  degree  of  subcooling  (Fig.  4). 

Analyzing  the  boiling  process  on  a  microscopic  scale,  Abdelmessih,  et  al.  (Ref.  26)  found 
that  both  the  maximum  bubble  radius  and  the  bubble  lifespan  decrease  as  the  liquid  velocity 
increases  (Fig.  5).  Gunther  (Ref.  27)  found,  in  addition  to  decreasing  the  bubble  radius  and 
lifetime,  that  increasing  velocity  also  decreases  the  bubble  population.  Tippets  (Ref.  28)  found 
that  when  high  velocity  is  combined  with  large  subcoolings,  the  flow  adjacent  to  a  surface 
sustaining  boiling  is  an  irregular,  frothy  layer  made  up  of  bubbles  that  grow,  slide,  and  then 
collapse  along  the  surface.  These  bubble  characteristics  make  their  detection  and  measurement 
at  high  flow  velocity  difficult. 

Subcooling 

Subcooling,  or  the  amount  that  the  fluid  bulk  temperature  is  below  the  fluid  saturation 
temperature,  has  been  shown  to  have  an  effect  similar  to  velocity  on  the  cooling  process. 
Researchers  have  correlated  experimental  data  using  a  fluid  bulk  temperature  measured  at 
various  locations  in  a  given  apparatus.  Generally,  for  larger  length-to-diameter  ratios  where 
the  bulk  fluid  may  enter  at  subcooled  conditions  but  reach  saturation  conditions  somewhere 
along  the  heated  length,  the  inlet  subcooling  (inlet  saturation  temperature  minus  the  inlet 
bulk  temperature)  is  typically  used  in  analyses.  For  constant  diameter  configurations  and 
uniform  heating  where  burnout  occurs  at  or  near  the  exit  of  the  heater,  the  exit  subcooling 
(exit  saturation  temperature  minus  the  exit  bulk  temperature)  is  typically  used.  Some  researchers 
have  used  a  mean  average  bulk  temperature  (average  of  the  inlet  and  outlet  bulk  temperatures) 
in  some  instances  to  define  subcooling.  In  some  cases,  subcooling  is  calculated  using  an  inlet 
bulk  temperature  and  an  exit  pressure  (exit  saturation  temperature).  For  a  configuration  such 
as  an  arc  heater  nozzle  where  nonuniform  heating,  velocity,  and  pressure  exist  along  the  length 
of  the  nozzle,  a  local  subcooling  using  a  local  bulk  temperature  and  saturation  temperature 
is  more  appropriate.  A  discussion  of  the  general  subcooling  effects  follows,  with  an  effort 
made  to  identify  the  type  of  subcooling  used  in  the  experimental  results  cited. 

McAdams,  et  al.  (Ref.  24)  showed  that  inlet  subcooling  has  a  slight  effect  on  heat  transfer 
in  the  pure  convection  and  the  nucleate  boiling  regions  of  the  boiling  curve  because  of  small 
variations  in  physical  properties  of  the  fluid  at  different  temperatures.  Subcooling,  like  velocity, 
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has  its  primary  advantages  at  boiling  crisis.  Van  Huff  and  Rousar  (Ref.  29)  and  Kreith  (Ref. 
30)  indicated  that  subcooling  and  velocity  are  the  most  important  variables  on  burnout  (CHF). 
Green,  et  al.  (Ref.  31)  measured  burnout  heat  flux  for  vertical  water  flow  in  a  rectangular 
channel  at  high  pressure  and  found  that  increasing  both  velocity  and  inlet  subcooling  increased 
the  burnout  heat  flux  (Fig.  6).  Van  Huff  and  Rousar  (Ref.  29)  combined  velocity  and  an 
average  subcooling  and  derived  a  linear  fit  to  previous  data  (Fig.  7).  Macbeth  (Ref.  18)  found 
that  increasing  inlet  subcooling  increased  CHF  linearly  at  constant  mass  velocity;  however, 
Katto  and  Yokoya  (Ref.  32)  indicated  that  nonlinear  regimes  can  occur  at  high  mass  velocity 
at  low-to-moderate  inlet  subcooling,  and  Lienhard  (Ref.  33)  found  that  the  maximum  heat 
flux  (CHF)  versus  degree  of  subcooling  curve  “flattens”  out  for  any  velocity  above  a  threshold 
degree  of  subcooling.  In  any  event,  it  is  likely  that  the  CHF  generally  can  be  increased  by 
increasing  velocity  and  subcooling. 

On  a  microscopic  scale,  Gunther’s  work  (Ref.  27)  with  a  flow  boiling  apparatus  found 
that  the  bubble  lifetime,  surface  coverage,  and  maximum  radius  decreased  with  increased 
exit  subcooling,  while  the  bubble  population,  which  initially  decreases,  increases  for  subcoolings 
above  72°C  (130°F  as  shown  in  Fig.  8).  As  reported  by  Westwater  (Ref.  12),  the  motion 
pictures  of  Dew  (Ref.  34)  confirmed  Gunther’s  findings  except  bubble  population,  which 
decreased  with  increased  average  subcooling.  Ibrahim  and  Judd  (Ref.  35)  confirmed  that 
the  growth  period  decreased  with  increased  subcooling  and,  in  addition,  found  that  the  waiting 
period  between  bubble  generations  from  an  active  nucleation  site  also  decreased  with  increased 
subcooling  except  at  low  subcooling  (below  22°C).  It  should  be  noted,  however,  that  Ibrahim 
and  Judd  performed  their  subcooling  studies  on  a  pool  boiling  apparatus. 

Pressure 

The  primary  contribution  of  pressure  on  a  cooling  process  is  its  effect  on  the  saturation 
temperature  of  the  fluid,  although  an  interesting  inversion  of  pressure  effect  occurs  at  boiling 
crisis.  Kreith  and  Summerfield  (Ref.  22),  who  studied  flow  boiling  with  butanol,  found  that 
increasing  pressure  had  no  effect  on  heat  transfer  in  the  pure  convection  region,  but  increased 
wall  temperature  for  a  given  heat  flux  in  the  nucleate  boiling  region,  as  shown  in  Fig.  2. 
An  increase  in  wall  temperature  is  not  particularly  beneficial;  however,  some  improvement 
to  the  CHF  can  be  achieved  by  increasing  pressure  -  up  to  a  point.  Cichelli  and  Bonilla  (Ref. 

36) ,  while  studying  pool  boiling  of  organic  liquids,  found  that  heat  flux  increased  with 
increasing  pressure  to  a  maximum  value  which  occurred  at  about  one-third  of  the  critical 
pressure  (Fig.  9).  Addoms  (Ref.  15)  obtained  similar  results  from  pool  boiling  of  water; 
however,  the  maximum  peak  heat  flux  varied  with  heater  wire  size  (Fig.  10).  Chang  (Ref. 

37)  suggested  a  similar  maximum  CHF  in  flow  boiling  that  occurs  at  various  pressures  between 
500  and  1,000  psia,  depending  on  mass  velocity  and  subcooling.  Boyd  (Ref.  1)  suggested 
a  similar  “optimum  pressure”  as  shown  in  Fig.  11. 
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Tolubinsky  and  Kostanchuk  (Ref.  38)  found  that  increased  pressure  decreased  the 
maximum  bubble  diameter  and  the  generation  frequency  (Fig.  12),  although  the  former  appears 
to  remain  constant  or  to  begin  increasing  above  a  pressure  of  8  to  10  bar  (120  to  150  psia). 
The  bubble  life  span  is  inversely  related  to  the  generation  frequency,  and,  therefore,  increases 
with  increased  pressure. 

Surface  Roughness 

Surface  roughness  plays  a  complicated  role  in  cooling  efficiency  and  is  divided  into  micro¬ 
roughness  (typically  smaller  than  3  /un)  and  macro-roughness  (artificial  protrusions  such  as 
knurling,  fins,  or  ribs).  Corty  and  Foust  (Ref.  39)  performed  an  extensive  micro-roughness 
study  of  various  fluids  in  pool  boiling  on  copper  and  nickel  surfaces.  They  found  that  increased 
roughness  had  no  effect  in  the  pure  convection  region,  but  heat-transfer  coefficients  were 
more  than  doubled  in  the  nucleate  boiling  region  for  the  roughest  surface  when  compared 
to  the  smoothest  surface  (Fig.  13).  Brown  (Ref.  40),  however,  found  a  modest  effect  of  surface 
finish  in  subcooled,  low-velocity  flow  boiling  at  moderate  heat  flux  (Fig.  14).  Berenson  (Ref. 
41)  and  later  Ramilison  and  Lienhard  (Ref.  42)  found  a  modest  surface  finish  dependence 
of  boiling  burnout  (up  to  approximately  20-percent  decrease  from  smooth  to  rough  surface) 
in  pool  boiling.  Other  than  Weatherhead  (Ref.  43),  who  concluded  that  CHF  is  decreased 
for  very  smooth  surfaces,  other  investigators  such  as  Leung,  et  al.  (Ref.  44),  Aladyev,  et 
al.  (Ref.  45),  and  Bergles  and  Morton  (Ref.  46)  found  little  or  no  effect  of  surface  micro¬ 
roughness  on  CHF  in  flow  boiling.  Bergles  and  Morton  further  determined  that  any 
improvements  of  the  burnout  heat  flux  caused  by  roughness  appear  to  diminish  with  increasing 
velocity  and  subcooling. 

Macro-roughness  has  been  found  to  have  a  more  pronounced  effect  on  forced  flow  cooling 
efficiency  than  the  micro-roughness  previously  discussed.  Types  of  macro-roughness  used 
in  forced  flow  include  knurling  (Durant,  et  al.,  Ref.  47),  fins  (Kovalev,  et  al.,  Ref.  48,  and 
Shim,  et  ah,  Ref.  49),  and  ribs  (Ravigururajan  and  Bergles,  Ref.  50,  and  Akhanda  and  James, 
Ref.  51).  Durant  and  his  team  determined  that  liquid  film  heat-transfer  coefficients  of  knurled 
surfaces  were  up  to  75-percent  higher  than  those  of  smooth  surfaces.  Burnout  heat  flux  was 
shown  to  increase  by  as  much  as  80  percent  (Fig.  15).  The  types  of  roughness  investigated 
by  Durant’s  team  were  comparatively  large  (greater  than  130  gm)  over  those  discussed  in 
the  previous  paragraph.  Macro-roughness  is  considered  a  heat-transfer  enhancement  technique 
rather  than  a  basic  parameter  affecting  heat  transfer  that  would  be  inherent  in  a  given  system. 
An  excellent  review  of  this  enhancement  technique  is  presented  by  Webb  (Ref.  52);  however, 
most  of  the  surface  geometries  have  been  developed  for  pool  boiling  only. 
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Coolant  Properties 

The  understanding  of  the  effects  of  coolant  properties  on  the  heat-transfer  processes  in 
boiling  remains  limited,  even  though  selection  of  the  fluid  is  a  most  important  consideration 
in  boiling  applications  (Ref.  1).  The  physical  properties  of  a  given  fluid  such  as  liquid/vapor 
density  ratio,  latent  heat  of  vaporization,  surface  tension,  heat  capacity,  and  transport 
properties  (i.e.,  viscosity  and  thermal  conductivity)  may  play  a  role  in  the  cooling  process; 
however,  distinguishing  their  individual  contribution  is  difficult,  and  little  or  no  data  exist 
to  ascertain  their  effect.  The  properties  are  related  to  the  pressure  and/or  the  temperature 
of  the  fluid,  and  the  problem  of  parametric  distortion  discussed  earlier  becomes  evident. 

The  limited  studies  of  the  microscale  mechanisms  during  boiling  have  produced  conflicting 
information  as  to  the  effect  of  latent  heat.  A  number  of  investigations  (Refs.  53-56)  have 
revealed  that  latent  heat  transport  accounts  for  as  much  as  50  percent  of  the  measured  heat 
flux  in  subcooled  boiling;  however,  Forster  and  Grief  (Ref.  57),  Brown  (Ref.  40),  and  Plesset 
and  Prosperetti  (Ref.  58)  found  very  little  latent  heat  contribution. 

In  their  work  with  a  pool  boiling  apparatus,  Pike,  et  al.  (Ref.  59)  found  that  increasing 
viscosity  hindered  the  onset  of  boiling  (decreased  heat-transfer  coefficient,  increased  wall 
superheat)  and  lowered  the  CHF.  As  will  be  discussed  under  transient  operation  effects,  it 
has  been  proposed  that  transient  operation  is  viscosity-dependent,  thereby  contributing  to 
an  observed  difference  in  CHF  for  the  transient  and  steady-state  heating  of  certain  fluids. 

Surface  tension  may  or  may  not  play  a  beneficial  role  in  the  cooling  process  (Ref.  11). 
Fiori  and  Bergles  (Ref.  20)  identified  several  investigations  that  clearly  showed  that  CHF 
decreases  with  decreasing  surface  tension.  The  addition  of  various  amounts  of  surface-active 
agent  has  been  found  to  increase  the  heat  transfer  quite  significantly  at  the  expense  of  reducing 
the  CHF.  For  example,  Frost  and  Kippenhan  (Ref.  54)  showed  that  for  a  given  wall 
temperature,  the  heat  flux  increased  by  approximately  50  percent  when  a  surface-active  agent 
was  added  to  water,  thereby  reducing  the  surface  tension  (Fig.  16).  The  reduction  of  surface 
tension  (or  increased  wettability)  had  a  detrimental  (decreasing)  effect  on  the  CHF,  which 
agrees  with  previous  results,  and  is  evident  in  Fig.  16.  However,  a  recent  study  by  Sadasivan 
and  Lienhard  (Ref.  60)  showed  an  increased  CHF  with  the  addition  of  the  surfactant  sodium 
dodecyl  sulfate  to  water.  They  also  noted  that  bubble  diameters  were  smaller,  and  the  number 
of  bubbles  was  much  greater  with  the  surfactant  present. 

Improvements  to  the  CHF  have  been  demonstrated  in  a  limited  number  of  investigations 
using  binary  mixture  components.  Papaioannou  and  Koumoutsos  (Ref.  61)  and  Wei  and 
Maa  (Ref.  62)  found  that  polymeric  additives  to  water  increased  heat-transfer  coefficients 
and  broadened  the  peak  (CHF)  of  the  boiling  curve.  Tolubinskiy  and  Matorin  (Ref.  63)  found 
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that  for  ethanol-water  and  acetone- water  mixtures,  the  maximum  CHF  (occurring  at  20-percent 
ethanol/80-percent  water  and  10-percent  acetone/90-percent  water)  exceeded  the  CHF  for 
water  by  30  and  40  percent,  respectively  (Fig.  17). 

A  number  of  coolants  such  as  refrigerants,  cryogenics,  organics,  binary  mixtures,  liquid 
metals,  and  propellants,  in  addition  to  water,  have  been  investigated  for  various  cooling 
applications.  A  propellant  is  probably  not  attractive  for  use  as  a  coolant  in  a  high-enthalpy 
facility,  except  in  specific  test  article  configurations.  The  use  of  organics  is  unlikely,  except 
possibly  as  a  trace  additive  or  in  a  binary  mixture,  because  of  their  lower  performance 
compared  to  that  of  water  (Ref.  64)  and  their  tendency  to  decompose  at  surface  temperatures 
experienced  at  nucleate  boiling  heat  flux  levels  (Ref.  8).  Binary  mixtures,  as  discussed 
previously,  may  be  an  attractive  approach  to  improve  the  CHF  in  a  water-cooled  system. 
Liquid  metals  may  be  considered  in  the  future  as  an  alternative  or  complement  to  water  cooling. 
Because  of  their  high  thermal  conductivity,  liquid  metals  have  proved  useful  in  convective 
cooling  applications.  Boiling  in  liquid  metals,  however,  is  somewhat  confused.  Lyon,  et  al. 
(Ref.  65)  showed  that  mercury,  for  instance,  appears  to  enter  film  boiling  at  surprisingly 
low  wall  superheats.  And  because  of  the  high  wettability  of  liquid  metals,  high  superheats 
are  normally  required  to  initiate  boiling,  in  some  cases  having  an  explosive  transition  that 
can  cause  structural  damage  to  hardware.  A  thorough  understanding  of  boiling  metals  is 
needed  before  they  can  be  applied  in  high  heat  flux/cooling  situations.  Refrigerants  such 
as  Freon®  have  been  used  in  previous  experiments  (Refs.  66-69)  to  simulate  water  cooling 
at  high  heat  flux  conditions  since  the  CHF  can  be  achieved  at  much  lower  heat  flux.  Scaling 
methods  are  applied  to  the  simulation  to  achieve  high  heat  flux/water  cooling  predictions 
(see,  for  example,  Ahmad,  Ref.  70).  A  good  example  of  this  “fluid-to-fluid”  scaling  is 
presented  by  Bergles  in  Ref.  8.  Use  of  refrigerants  for  actual  cooling  in  high-enthalpy  facilities 
is  not  recommended  because  of  their  low  thermal  conductivity  and  CHF  levels.  In  addition, 
refrigerants  suffer  from  a  high  wettability  problem  not  unlike  that  of  liquid  metals.  Similarly, 
Powell  (Ref.  71)  found  that  low  heat -transfer  coefficients  and  CHF  levels  existed  for  cryogenics 
such  as  hydrogen,  nitrogen,  and  oxygen  which  would  limit  the  use  of  liquified  gases  as  a 
heat -transfer  fluid.  Refrigerants  and  cryogenics  perform  well  in  systems  that  have  relatively 
low  heat  flux  requirements  and  limited  space  availability. 

Wall  Material 

Although  it  is  considered  to  be  a  secondary  effect,  surface  material  properties  appear 
to  contribute  to  cooling  efficiency.  Early  experiments  in  pool  boiling  demonstrated  a  difference 
between  stainless  steel  and  copper  in  nonboiling  and  boiling  water  (Ref.  72);  copper,  gold, 
and  chromium  in  boiling  ethanol  (Ref.  73);  nickel,  tungsten,  and  Chromel®  in  boiling  water 
(Ref.  74);  and  aluminum,  copper,  iron,  and  chrome-plated  copper  in  water  and  various  organic 
liquids  (Ref.  75),  among  others.  Later,  Magrini  and  Nannei  (Ref.  76)  found  that  the  heat- 
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transfer  coefficient  appears  to  be  dependent  on  the  thermal  activity,  Vgck,  for  heaters  of 
different  metals  above  a  limiting  value  of  thickness  (Fig.  18).  Kovalev,  et  al.  (Ref.  77)  showed 
that  copper  with  coatings  of  low  thermal  conductivity  tended  to  increase  the  degree  of  wall 
superheat  over  that  of  clean  copper  boiling  in  Freon  113.  Cheng  and  Ragheb  (Ref.  78)  found 
that  a  greater  wall  superheat  was  required  for  an  Inconel®  surface  than  a  copper  surface 
in  flow  boiling  of  water  (Fig.  19). 

Numerous  pool  boiling  experiments  have  demonstrated  that  the  CHF  has  a  significant 
material  dependence  (Refs.  73-75  and  77)  which  has  not  been  observed  in  flow  boiling.  The 
data  of  Cheng  and  Ragheb  presented  in  Fig.  19  show  the  CHF  to  be  equivalent  for  Inconel 
and  copper.  Jacket,  et  al.  (Ref.  79)  found  no  effect  on  CHF  due  to  changing  wall  material 
from  nickel  to  Zircaloy®  or,  later,  stainless  steel  (Ref.  31)  in  flow  boiling  of  water  at  high 
pressure.  Van  Huff  and  Rousar  (Ref.  29)  identified  studies  in  which  various  metals  had  no 
effect  on  CHF,  except  where  the  material  and  coolant  were  incompatible  (additional  reactions 
occur).  Similarly,  Vandervort,  et  al.  (Ref.  80)  found  no  difference  in  CHF  for  subcooled 
forced-convection  boiling  of  water  on  stainless  steel  304,  stainless  steel  316,  nickel  200,  brass, 
and  Inconel  600.  However,  Vandervort’s  team  suggested  that  materials  with  low  values  of 
thermal  conductivity,  such  as  stainless  steel  and  Inconel,  may  exhibit  hot  spots  due  to  changes 
in  metal  crystalline  structure  when  high  wall  temperature  gradients  are  present  (e.g.,  high 
heat  flux  in  thin  wall  tubes).  Fiori  and  Bergles  (Ref.  20)  suggest  that  materials  with  high 
thermal  diffusivity  (e.g.,  aluminum  and  molybdenum)  are  expected  to  have  higher  CHF. 

Wall  Thickness 

The  thickness  of  the  heated  wall  has  been  shown  to  have  an  effect  on  the  boiling  heat 
flux,  but  only  below  a  critical  or  limiting  value  of  wall  thickness.  Magrini  and  Nannei  (Ref. 
76),  for  instance,  found  for  saturated  pool  boiling  of  water  that  the  limiting  value  was 
approximately  15  /am  (0.0006  in.)  for  tin  and  nickel,  and  negligible  for  copper  and  silver. 
They  also  determined  the  limiting  value  for  zinc  to  be  approximately  70  /xm  (0.0028  in.), 
although  it  is  not  obvious  from  the  data  presented  in  Ref.  76.  Del  Valle  and  Kenning  (Ref. 
81),  while  studying  subcooled  water  flow  at  high  heat  flux,  found  that  thin  heated  walls  (0.08- 
to  0.2-mm  thickness)  had  no  effect  in  the  nonboiling  region,  but  the  rate  of  nucleate  boiling 
heat  transfer  at  a  given  wall  superheat  increased  with  increasing  wall  thickness  (Fig.  20).  Tippets 
(Ref.  28)  found  that  CHF  measured  for  0.15-mm-thick  heater  ribbons  was  approximately 
20  percent  lower  than  for  0.25-mm-thick  ribbons  in  saturated  water  flow.  They  suggest  that 
thicker  walls  can  maintain  a  higher  internal  heat  generation  rate  per  unit  volume  (q"/V)  before 
the  critical  condition  is  reached  than  thin  walls,  since  the  latter  will  have  a  larger  temperature 
rise  at  a  given  heat  flux  (i.e.,  q"/V  will  be  greater  for  a  thin-walled  tube  since  the  volume 
is  smaller  and  the  heat  flux  constant).  Fiori  and  Bergles  (Ref.  20),  using  subcooled  water 
in  a  flow  boiling  apparatus,  showed  that  thicker-walled  tubes  (2  mm  versus  0.3  mm)  had 
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a  CHF  up  to  58  percent  higher  than  thin-walled  tubes.  Aladyev,  et  al.  (Ref.  44),  however, 
found  no  effect  of  wall  thickness  (0.4  mm  to  2  mm)  on  the  CHF  for  water  flow  boiling. 

Geometry 

Geometric  considerations  that  have  been  experimentally  shown  to  have  an  effect  on  the 
cooling  process  include  the  diameter  (or  hydraulic  diameter)  of  the  heater  (or  channel),  heated 
length  (or  L/D),  the  channel  cross  section  and  longitudinal  configuration,  and  the  flow 
orientation  (Ref.  1).  Because  of  various  applications,  past  experiments  associated  with  studying 
cooling  problems  have  encompassed  a  wide  range  of  configurations,  and  one  must  be  cautious 
when  comparing  data  from  one  experiment  to  another  (Ref.  11). 

The  channel  diameter  appears  to  have  a  varied  effect  on  boiling  heat  transfer.  Schweppe 
and  Foust  (Ref.  82)  found  that  for  flow  boiling  with  low-velocity  water,  the  wall  superheat 
increased  (or  heat-transfer  coefficient  decreased)  as  channel  diameter  increased  from  11.1 
mm  to  26.6  mm.  Collier  (Ref.  7)  presented  subcooled,  flow  boiling  data  of  Lee  and  Obertelli 
(Ref.  83)  and  Matzner  (Ref.  84)  that  indicated  that  CHF  increased  as  channel  diameter 
increased  from  5.6  mm  to  37.5  mm  at  constant  subcooling  (Fig.  21).  Bergles  (Ref.  85),  however, 
found  for  subcooled  flow  boiling  that  CHF  increased  as  channel  diameter  decreased  from 
8.4  mm  to  0.6  mm  (see  Fig.  22),  although  the  percentage  increase  in  CHF  is  somewhat  less 
when  parametric  distortion  is  accounted  for  (Ref.  20).  Similarly,  Glushchenko  (Ref.  86)  found 
that  CHF  increased  as  the  tube  diameter  was  reduced  below  2  mm.  Van  Huff  and  Rousar 
(Ref.  29)  found  no  effect  of  channel  diameter  on  CHF  for  subcooled  flow  boiling  with 
diameters  ranging  from  2.7  to  14.9  mm. 

Schweppe  and  Foust  (Ref.  82),  by  varying  channel  diameter  (as  previously  discussed)  and 
holding  the  heated  length  constant,  showed  for  flow  boiling  with  low-velocity  water  that  the 
wall  superheat  increased  as  L/D  decreased  below  15  to  20.  Bergles  (Ref.  85)  found  for  flow 
boiling  that  CHF  increased  as  L/D  was  decreased  below  35  (Fig.  23).  Similarly,  Ornatskiy 
(Ref.  87),  using  subcooled  water  flowing  in  circular  ducts,  found  this  limiting  L/D  to  be 
between  20  to  24,  and  Jens  and  Lottes  (Ref.  88)  showed  that,  for  high-pressure,  subcooled 
flow  boiling,  changing  L/D  from  110  to  21  had  an  increasing  effect  on  CHF.  Lee  and  Obertelli 
(Ref.  83)  found  a  considerable  effect  on  CHF  in  subcooled,  low-velocity  flow  for  L/D  values 
less  than  300  (Fig.  24).  Van  Huff  and  Rousar  (Ref.  29),  however,  noted  that  data  acquired 
by  Aerojet  General  Corporation  for  high-velocity,  highly  subcooled  water  showed  very  little 
effect  of  L/D  (ranging  from  13  to  100)  on  CHF.  As  Boyd  (Ref.  1)  points  out,  since  the  L/D 
limit  is  related  directly  to  the  flow  development,  this  limit  is  not  a  constant  and  is  related 
to  the  flow  parameters  and  fluid  properties. 

The  geometric  considerations  thus  far  discussed  have  been  primarily  for  circular  tube 
configurations.  External  flow  over  rods  (or  wires)  and  rod  bundles,  jets,  internal  flow  in 
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multiple  tube  (or  tube  bundle)  configurations,  annuli,  and,  to  a  lesser  extent,  rectangular 
channels  have  been  studied  in  the  past.  Internal  tube  flow  and  annular  flow  closely  simulate 
the  cooling  configurations  encountered  in  a  high-enthalpy  arc  facility.  Becker  and  Hernborg 
(Ref.  89)  found  for  water  flowing  in  an  annulus  with  inlet  subcooling  and  an  inner  heated 
wall  (similar  to  an  arc  facility  nozzle)  that  CHF  data  were  considerably  below  those  for 
comparable  flow  in  a  circular  tube,  and  Zenkevich,  et  al.  (Ref.  90)  also  found  a  difference 
in  CHF  between  annuli  and  tubes  (at  the  same  length  and  diameter).  Barnett  (Ref.  91), 
however,  points  out  significant  difficulties  which  arise  when  comparing  tube  and  annulus 
data  and  the  definition  of  an  appropriate  equivalent  diameter.  In  addition,  Alekseev,  et  al. 
(Ref.  92)  found  significant  heated  length  effects  for  annuli  at  much  greater  lengths  than  for 
circular  tubes.  Collier  (Ref.  7)  discusses  other  characteristics  of  annular  flow  of  interest. 
Specifically,  eccentricity  of  the  internal  heated  wall  will  tend  to  reduce  the  CHF  (see  also 
Ref.  91).  In  addition,  Collier  points  out  that  the  outer  wall  in  a  concentric  annulus  system 
has  little  or  no  effect  on  the  CHF  for  the  inner  surface. 

The  channel  longitudinal  configuration  is  important  to  the  boiling  heat  transfer  when 
significant  curves  or  twists  are  introduced.  Such  longitudinal  changes  have  an  effect  on  the 
acceleration  and  turbulence  of  the  coolant  flow.  A  discussion  of  these  types  of  configurations 
is  included  in  the  section  entitled  “Acceleration  and  Turbulence.” 

Coolant  flow  orientation  (horizontal,  inclined,  or  vertical)  has  been  shown  to  have  an 
effect  on  cooling  efficiency  due  primarily  to  buoyancy  (gravity)  causing  stratification  of  the 
liquid/vapor  flow  (Refs.  93  and  94).  The  effect  of  flow  orientation,  though,  diminishes  with 
increasing  subcooling  and/or  mass  velocity,  and  there  appears  to  be  a  threshold  above  which 
no  effect  of  orientation  has  been  observed  (Refs.  1,  7,  94,  and  95).  Zeigarnik,  et  al.  (Ref. 
96),  for  instance,  showed  that  the  orientation  of  the  heater  plate  in  a  rectangular  flow  channel 
had  little  effect  on  the  CHF  at  mass  velocities  above  4,500  kg/m2  sec  (922  lbm/ft2  sec)  in 
subcooled  water  flow  (Fig.  25).  Boyd,  et  al.  (Ref.  97)  suggests  that  the  mass  velocity  limit 
is  between  3,000  and  4,500  kg/m2  sec  (600  to  900  lbm/ft2  sec),  which  agrees  with  the  velocity 
limit  of  4,000  kg/m2  sec  (820  lbm/ft2  sec)  found  by  Merilo  (Ref.  94). 

Gas  Content 

A  number  of  experiments  in  pool  boiling  (Refs.  59  and  98)  and  flow  boiling  (Refs.  21, 
24,  and  99)  have  demonstrated  the  effect  of  dissolved  gases  on  pure  convection  and  nucleate 
boiling  heat  transfer.  McAdams,  et  al.  (Ref.  24),  for  example,  found  little  effect  of  dissolved 
air  on  forced-convection  heating,  but  the  presence  of  air  tended  to  reduce  the  wall  superheat 
necessary  to  initiate  and  sustain  nucleate  boiling  (Fig.  26)  which  is  typical*  of  data  obtained 


*  Only  Jens  and  Lottes  (Ref.  88)  found  that  boiling  incipience  is  not  affected  by  gas  content. 
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in  boiling  experiments  where  the  effect  of  dissolved  gas  was  evaluated.  In  several  experiments, 
the  effect  of  dissolved  gases  was  found  to  diminish  in  fully  developed  nucleate  boiling  (Refs. 
24  and  100).  This  is  evident  in  Fig.  26  at  high  wall  superheat  where  the  data  for  boiling  with 
dissolved  air  match  that  for  degassed  boiling. 

Data  obtained  to  evaluate  the  effect  of  dissolved  gases  on  CHF  are  contradictory;  however, 
in  experiments  where  an  effect  was  noted  (Refs.  24,  99-101),  the  presence  of  dissolved  gases 
tended  to  reduce  the  CHF  from  that  of  degassed  boiling  at  the  CHF.  Fisenko,  et  al.  (Ref. 
99)  found  this  reduction  to  be  as  high  as  23  percent  at  higher  subcoolings.  A  number  of 
investigators  (e.g.,  Refs.  21,  31,  and  102)  identified  no  appreciable  effect  of  dissolved  gases 
on  CHF.  Kalayda,  et  al.  (Ref.  103)  found  that  for  water  flow  boiling,  the  gas  content  had 
little  effect  on  CHF  at  pressures  below  100  bar,  but  at  higher  pressures  the  CHF  was  reduced 
with  increased  concentrations  of  nitrogen. 

Other  experimental  results  have  been  obtained  with  gassed  and  degassed  coolants.  Jens 
and  Lottes  (Ref.  88)  and  Buchberg,  et  al.  (Ref.  100)  in  water  flow  boiling  found  that  the 
presence  of  dissolved  nitrogen  had  a  small  effect  on  the  pressure  loss  along  the  heated  wall. 
Baranenko,  et  al.  (Ref.  104)  found  for  atmospheric  pool  boiling  of  water  that  the  bubble 
diameters  are  larger  when  degassed  water  is  used  rather  than  gassed  water.  The  presence  of 
dissolved  oxygen  can  have  an  additional  effect  on  the  heat-transfer  process  due  to  oxidation 
of  surfaces  in  contact  with  the  coolant.  A  discussion  of  this  oxidation  effect  is  included  in 
the  following  section. 

Surface  Aging,  Deposits,  and  Coatings 

Surface  aging*  has  been  shown  in  the  past  to  increase  the  wall  superheat  in  pool  boiling 
(Refs.  11  and  13).  Bonilla  and  Perry  (Ref.  73),  for  instance,  found  this  to  be  the  case  on 
a  chrome-plated  copper  heater  with  boiling  ethanol.  Hughes  (Ref.  21)  found  for  subcooled 
flow  boiling  that  surface  aging  tended  to  decrease  the  CHF,  and  Akhanda  and  James  (Ref. 
51)  found  that  aging  increased  the  wall  superheat  in  forced-convection  boiling  of  water. 

The  effect  of  surface  deposits  (or  fouling)  on  cooling  heat  transfer  is  much  more  difficult 
to  ascertain.  Epstein  (Ref.  105)  identified  several  categories  in  which  fouling  (primarily  in 
forced-flow  heat  exchangers)  can  occur.  These  include  precipitation  (dissolved  substance 
deposits),  particulates  (suspended  particle  deposits),  chemical  reaction  (surface  is  not  a 


*  Nucleate  boiling  on  clean  surfaces  will  not  reach  a  steady-state  condition  for  a  sizeable  period  of  time  (in  some  cases,  an 
hour  or  more)  due  to  several  reasons  including  surface  outgassing,  chemical  reactions,  etc.  (Ref.  1 1).  This  process  of  reaching 
steady-state  conditions  is  known  as  aging. 
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reactant),  and  corrosion  (surface  reactions)  among  others.  Collier  (Ref.  7)  pointed  out  that 
porous  deposits  on  heat-transfer  surfaces  can  influence  the  CHF  both  favorably  and  adversely 
(see  also  Westwater,  Ref.  11).  Bui  and  Dhir  (Ref.  106),  found  that  the  presence  of  an  oxide 
increased  the  CHF  in  pool  boiling  of  water.  Mel’nikov,  et  al.  (Ref.  107),  however,  found 
that  CHF  decreased  nearly  36  percent  (at  a  quality  of  0.06)  with  the  presence  of  a  55-/xm 
(0. 002-in. )-thick  ferric-oxide  scale  as  compared  to  a  scale-free  surface  in  forced-convection 
boiling  of  water. 

The  effect  of  surface  coatings  (primarily  used  for  heat -transfer  enhancement)  on  cooling 
heat  transfer  is  similar  to  that  of  deposits.  For  example,  Zhukov,  et  al.  (Ref.  108)  found 
that  the  wall  superheat  and  CHF  primarily  decreased  for  an  aluminum  oxide  coating,  but 
both  increased  for  an  enameled  surface  in  Freon  pool  boiling.  As  discussed  previously  in 
the  wall  material  effects,  Kovalev,  et  al.  (Ref.  77)  found  that  low  thermal  conductive  coatings 
increased  the  wall  superheat.  Collier  (Ref.  7)  discusses  improvements  in  heat  transfer 
demonstrated  by  several  investigators  using  thin  Teflon®  films;  however,  Sadasivan  and 
Lienhard  (Ref.  60)  showed  that  Teflon  microcoatings  on  Nichrome®  surfaces  reduced  the 
wettability  and,  hence,  reduced  CHF  values. 

Acceleration  and  Turbulence 

High-enthalpy  facility  components  invariably  require  coolant  flow  direction  and  area 
changes  which,  in  turn,  cause  flow  acceleration/deceleration  and  turbulence.  Merte  and  Clark 
(Ref.  109)  found  for  pool  boiling  of  saturated  water  that  increasing  the  acceleration  (g-loading 
of  1  to  20)  decreased  the  wall  temperature  and  thereby  increased  the  heat -transfer  coefficient 
in  both  the  pure  convection  and  nucleate  boiling  regions,  although  the  effect  of  acceleration 
diminished  as  the  boiling  became  fully  developed.  Forced-convection  acceleration/turbulence 
effects  have  been  studied  with  curved  channels  and  coiled  tubes.  Gu,  et  al.  (Ref.  110),  studying 
a  fluorocarbon  liquid  flowing  over  a  concave-shaped  heater,  found  that  for  the  same  flow 
velocity  and  subcooling,  both  the  heat-transfer  coefficient  and  CHF  are  higher  in  a  curved 
channel  (concave  section)  than  in  a  straight  channel.  Gu’s  team  also  found  that  the  increase 
in  CHF  is  significant  (40-percent  higher  than  straight  channel  CHF)  at  higher  flow  velocities. 
Hughes  (Ref.  21)  found  similar  results  for  concave  heater  surfaces  in  subcooled  Freon  113 
flow  boiling.  Hughes  found  that  the  CHF  was  increased  up  to  50  percent  in  concave  channels 
over  straight  channels  (Fig.  27).  Hughes’  data,  however,  illustrate  the  disadvantages  of  such 
a  curved  channel  configuration.  Should  the  channel  be  heated  entirely,  the  convex  portion 
could  have  a  detrimental  effect  on  CHF,  as  shown  in  Fig.  27.  The  CHF  for  convex  surfaces 
was  as  much  as  22  percent  lower  than  that  for  a  straight  channel.  Jensen  and  Bergles  (Ref. 
95)  found  low  CHF  levels  in  subcooled  R-113  flow  over  inside  surfaces  (convex  portion)  of 
helically  coiled  tubes,  and  Bergles  (Ref.  8)  noted  similar  results  found  by  other  investigators. 
Winovich  and  Carlson  (Ref.  Ill)  made  use  of  the  advantages  of  a  concave  surface  in  the 
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design  of  an  undulating  flow  channel  used  in  various  high-enthalpy  facility  components.  They 
found  that  the  undulating  channel  significantly  increased  (as  much  as  100  percent)  the  CHF 
over  that  of  a  straight  channel  (Fig.  28). 

Swirl  flow  devices  and  other  inserts  have  been  used  in  the  past  to  induce  turbulence,  thereby 
enhancing  the  cooling  heat  transfer  in  applications  such  as  heat  exchangers.  Kudryavtsev, 
et  al.  (Ref.  112),  used  a  twisted  tape  in  tubes  with  saturated  water  flow  boiling  to  show  that 
the  turbulence  generated  by  density  gradients  and  acceleration  from  the  twisted  tapes  decreased 
the  rate  of  nucleate  boiling  (increased  wall  superheat).  However,  Lopina  and  Bergles  (Ref. 
113)  found  little  effect  of  swirl  on  the  wall  superheat  for  subcooled  water  flow  boiling. 
Ornatskiy,  et  al.  (Ref.  114)  demonstrated  that  induced  swirl  at  the  inlet  of  a  concentric  annulus 
increased  the  CHF  by  10  to  60  percent  over  unswirled  flow  for  low  qualities.  Gambill,  et 
al.  (Ref.  115)  found  similar  results  for  subcooled  water  flow  boiling.  Ornatskiy’s  team  also 
showed  that  the  effect  of  swirl  on  CHF  decreases  with  increasing  pressure  (Fig.  29),  but 
increases  with  increasing  mass  velocity  (Fig.  30). 

Megerlin,  et  al.  (Ref.  116)  augmented  the  heat  transfer  in  subcooled  water  flow  in  tubes 
by  using  brush  and  mesh  inserts,  showing  that  heat-transfer  coefficients  up  to  nine  times 
the  coefficients  of  empty  tubes  can  be  obtained.  Megerlin  and  his  team  found  that  the  inserts 
produced  very  large  pressure  drops  and  very  large  wall  superheats.  Bergles  (Ref.  8)  recommends 
that  the  inserts  are  more  suitable  for  use  in  single-phase  flow  heat-transfer  augmentation. 
As  mentioned  previously  in  the  discussion  of  surface  roughness  effects,  macro-roughness 
can  be  used  to  enhance  the  heat  transfer  in  forced  flow,  which  probably  involves  some  type 
of  turbulence  mechanism. 

Heat  Flux  Distribution 

Variations  in  heat  flux  distribution,  both  axial  and  circumferential,  have  been  shown  to 
have  an  effect  on  boiling  heat  transfer,  primarily  on  the  CHF.  Stein  and  Begell  (Ref.  117), 
using  subcooled  water  flowing  in  internally  heated  annuli,  found  no  significant  effect  of  an 
axial  heat  flux  distribution  (cosine  shape)  on  the  heat-transfer  coefficients.  Various  researchers, 
however,  have  found  significant  effects  when  CHF  is  approached.  Styrikovich,  et  al.  (Ref. 
118),  while  studying  subcooled  and  saturated  water  flow  in  pipes  with  uniform  and  linearly 
varying  heat  flux,  found  that  the  CHF  for  a  linearly  increasing  heat  flux  (along  the  tube 
length)  was  nearly  twice  that  of  a  tube  with  uniform  heat  flux  for  low-quality  water  at  100 
atm  (mass  velocity  of  410  lbm/ft2  sec)  and  40  to  50  percent  higher  at  180  atm.  For  these 
cases,  the  crisis  always  occurred  at  the  exit  of  the  heated  tube.  Styrikovich’s  team  also  found 
that  the  CHF  for  a  linearly  decreasing  heat  flux  occurred  near  the  inlet  of  the  heated  section 
and  was  larger  than  that  with  a  uniform  heat  flux,  although  data  errors  prevented  the  increase 
from  being  quantified  for  the  subcooled  case.  Ornatskiy,  et  al.  (Ref.  1 19)  found  similar  results 
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for  subcooled  water  flow  in  annuli  with  linearly  or  parabolically  increasing  heat  flux;  however, 
they  found  quite  different  results  for  other  nonuniform  heat  flux  applications.  Earlier  (Ref. 
114),  Ornatskiy’s  team  had  shown  that  for  a  cosinusoidal  heat  flux  distribution  on  annuli 
with  subcooled  water  flow,  the  CHF  can  be  as  much  as  80  percent  lower  than  that  of  a 
uniformly  heated  annulus  at  similar  conditions,  and  the  difference  between  the  nonuniform 
heating  CHF  and  uniform  heating  CHF  increased  with  a  decrease  in  pressure  and  mass  flow 
rate.  In  addition,  data  presented  in  Ref.  119  show  that,  for  subcooled  water  flow  in  annuli, 
the  CHF  with  sinusoidal,  exponential,  or  linearly /parabolically  decreasing  heat  flux  is  generally 
lower  than  the  CHF  with  uniform  heat  flux.  Similar  results  were  presented  by  Zenkevich, 
et  al.  (Ref.  90)  for  subcooled  and  saturated  water  flowing  in  tubes  with  uniform  and 
cosinusoidal  heat  flux  distribution,  although  at  higher  pressure  and  subcooled  conditions, 
the  difference  appears  to  diminish  or  even  reverse  (Fig.  31).  They  point  out  that  for  complex 
nonuniform  heat  flux  applications,  the  location  of  the  burnout  point  along  the  length  of 
tube  or  annulus  is  difficult  to  predict.  Swenson,  et  al.  (Ref.  120)  presented  some  interesting 
results  concerning  saturated  water  flow  in  tubes  with  nonuniform  heating  along  their  length. 
They  showed  that  the  effect  of  nonuniform  heating  on  DNB  (CHF)  was  similar  to  that  of 
Ornatskiy’s  team  (CHF  for  nonuniform  heating  is  less  than  CHF  for  uniform  heating); 
however,  as  the  region  over  which  the  nonuniform  heating  occurred  became  very  small  (i.e., 
a  spike  heat  flux),  the  CHF  for  the  nonuniform  heating  case  approached  and  surpassed  that 
for  a  uniformly  heated  tube. 

Aladyev,  et  al.  (Ref.  45)  demonstrated  that  nonuniform  heating  around  the  perimeter 
of  a  tube  with  subcooled  or  saturated  water  flow  can  affect  the  level  of  CHF.  They  showed 
that  the  CHF  generally  occurs  where  the  heat  flux  is  the  highest  and  increases  with  increasing 
unevenness  of  the  heating  (characterized  by  the  ratio  of  the  maximum  to  the  average  heat 
flux).  It  was  also  shown  that  the  effect  of  nonuniform  heating  around  the  perimeter  decreased 
with  increasing  diameter.  Leontiev,  et  al.  (Ref.  121)  found,  for  subcooled  water  flowing  in 
horizontal  tubes  with  nonuniform  circumferential  heating  (the  ratio  of  maximum  heat  flux 
to  average  heat  flux  was  1.5),  that  the  local  values  of  CHF  were  somewhat  higher  than  those 
of  uniformly  heated  tubes.  However,  when  they  averaged  the  heat  flux  at  the  station  where 
the  critical  point  occurred,  the  averaged  CHF  for  the  nonuniform  heating  case  was  much 
lower  than  that  for  the  uniform  heating  case  at  subcooled  conditions,  and  the  two  coincided 
at  high-quality  conditions. 

Transient  Operation 

The  effects  of  transient  power  (heat  flux),  flow  rate,  and  pressure  operation  on  heat  transfer 
have  been  studied  in  previous  experiments  primarily  associated  with  petrochemical  and  nuclear 
power  plant  safety  (e.g.,  Ref.  122).  A  major  portion  of  these  experiments  involved  the  study 
of  transient  heat  flux  effects.  Ragheb,  et  al.  (Ref.  123),  showed  that  the  wall  superheat  was 


20 


AEDC-TR-93-3 


larger  for  transient  boiling  (found  by  quenching  a  high  thermal  capacity  tube  with  subcooled, 
low-pressure  water  flow)  than  steady-state  boiling  at  comparable  conditions  (Fig.  32).  Cheng, 
et  al.  (Ref.  124)  had  earlier  found  little  difference  in  wall  superheat  for  transient  or  steady- 
state  heating  with  the  same  apparatus,  although  a  limitation  in  the  steady-state  heating 
capability  prevented  data  comparison  at  heat  flux  levels  approaching  the  CHF.  Kataoka, 
et  al.  (Ref.  125)  found  an  effect  on  wall  superheat  similar  to  the  one  found  by  Ragheb’s 
team  for  an  exponentially  increasing  heat  input  to  a  platinum  wire  in  subcooled  forced  water 
flow. 

Several  interesting  results  have  been  found  during  studies  of  the  effect  of  transient  heat 
flux  on  CHF.  Borishanskiy  and  Fokin  (Ref.  64)  and  later  Tolubinskiy,  et  al.  (Ref.  126)  found 
that  the  CHF  was  equal  for  transient  and  steady-state  heating  in  pool  and  forced-convection 
boiling  of  water  with  and  without  subcooling.  However,  when  organic  coolants*  were  used, 
it  was  found  that  the  CHF  for  steady-state  heating  was  much  larger  than  the  CHF  for  transient 
heating  (Tolubinskiy’ s  team  found  this  difference  to  be  as  large  as  a  factor  of  3).  Tolubinskiy’s 
team  theorized  that  the  transient  operation  is  viscosity-dependent,  and  pointed  out  that  the 
organics  have  higher  viscosities  than  that  of  water,  thereby  contributing  to  the  difference 
in  CHF  for  the  transient  and  steady-state  heating  of  organic  fluids.  Roemer  (Ref.  127)  found 
no  effect  of  transient  heating  (time  constant  of  the  applied  power  was  approximately  0.1 
sec)  on  CHF  for  water  flow  normal  to  cylindrical  test  elements  as  long  as  the  test  surfaces 
were  aged.  This  result  was  found  to  be  independent  of  the  test  element  wall  thickness  and 
the  speed  of  the  power  transient.  The  CHF  for  unaged  specimens  was  as  much  as  25  percent 
lower  than  the  CHF  at  comparable  steady-state  conditions.  Other  researchers  have  found 
results  conflicting  with  those  observed  by  Borishanskiy’ s  team  and  Tolubinskiy’s  team  for 
water.  As  seen  in  Fig.  32,  the  data  of  Ragheb,  et  al.  (Ref.  123)  show  the  CHF  for  transient 
conditions  to  be  higher  than  that  for  steady-state  conditions.  Kataoka,  et  al.  (Ref.  125)  also 
found  the  CHF  for  transient  heating  to  be  greater  than  that  for  steady-state  heating,  and 
the  difference  between  the  two  CHF  varied  approximately  as  the  -0.6  power  of  the  period 
of  the  exponentially  increasing  heat  flux  (Fig.  33). 

Fewer  experiments  have  been  performed  to  evaluate  the  effects  of  transient  flow  rate  or 
pressure  operation  on  boiling  heat  transfer.  Celata,  et  al.  (Ref.  128),  using  subcooled 
Refrigerant- 12  as  a  coolant,  found  that  the  difference  between  the  CHF  for  a  flow  rate  transient 
and  the  CHF  for  a  steady-state  flow  rate  became  significant  when  the  time  to  reach  half 
of  the  initial  flow  became  less  than  approximately  twice  the  flow  rate  transit  time.  The  CHF 
for  the  transient  flow  was  typically  less  than  the  comparable  steady-state  CHF  in  these  cases. 
Celata,  et  al.  (Ref.  129)  later  evaluated  the  combined  effect  of  transient  flow  rate  and  thermal 
power  on  CHF  using  flowing  subcooled  R-12.  His  team  found  that  the  CHF  for  the  transient 


*  Borishanskiy’s  team  used  ethyl  alcohol;  Tolubinskiy’s  team  used  ethanol  and  acetone. 
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cases  was  approximately  equal  to  that  for  steady-state  conditions  when  the  time  to  reach 
half  flow  was  greater  than  the  transit  time.  When  the  time  to  reach  half  flow  dropped  below 
the  transit  time,  the  CHF  for  the  transient  case  became  greater  than  the  corresponding  steady- 
state  CHF.  These  trends  appeared  to  be  insensitive  to  the  level  of  pressure  or  mass  flux  at 
which  the  runs  were  made. 

Weisman,  et  al.  (Ref.  130),  using  depressurization  of  saturated  water  to  produce  boiling 
in  a  nonflowing  apparatus,  found  that  the  wall  superheat  for  the  pressure  transients  was 
significantly  higher  (as  high  as  a  factor  of  two)  than  a  slow  pressure  change,  although  the 
data  are  probably  applicable  only  near  boiling  incipience.  Celata,  et  al.  (Ref.  131)  used  flowing, 
subcooled  R-12  to  demonstrate  that  pressure  transients  tended  to  reduce  the  CHF  (higher 
depressurization  rates  led  to  more  marked  crisis  conditions)  from  that  exhibited  during 
comparable  steady-state  conditions.  Finally,  Celata,  et  al.  (Ref.  122)  evaluated  the  combined 
effect  of  transient  pressure,  power,  and/or  flow  rate  on  CHF,  although  no  steady-state  data 
were  obtained  for  comparison.  Celata’ s  team  did  show  the  inadequacy  of  using  the  available 
steady-state  CHF  correlations  in  predicting  transient  power/flow  rate/pressure  situations. 
The  data  acquired  by  Celata’s  team  for  one-,  two-,  and  three-parameter  transients  were  recently 
tabulated  in  Ref.  132. 

Flow  Instabilities 

Thermal-hydraulic  flow  instabilities  are  directly  related  to  the  flow  system  design,  and 
can  seriously  affect  the  value  of  CHF  and  the  integrity  of  the  flow  system  (Ref.  17).  Podowski 
(Ref.  133)  identified  several  instabilities  that  can  appear  in  two-phase  flow  and  pointed  out 
that  these  instabilities  can  be  compounded  in  some  cases.  The  classes  of  instabilities  include 
excursive  (Ledinegg),  flow  regime  relaxation,  nucleation,  density-wave  oscillation,  pressure 
drop  oscillation,  acoustic,  and  condensation-induced  instabilities.  Numerous  authors  (Refs. 
17,  133-136)  have  provided  excellent  reviews  of  these  instabilities  and  their  effects;  therefore, 
a  rigorous  treatment  is  excluded  here.  In  general,  flow  instabilities  have  a  detrimental  effect 
on  boiling  heat  transfer,  primarily  a  reduction  in  the  CHF  (Ref.  14). 

According  to  Podowski  (Ref.  133),  one  of  the  most  predominant  of  the  instabilities  is 
that  classified  under  density-wave  oscillation.  It  includes  flow  loop,  parallel-channel,  and 
channel-to-channel  instabilities.  Flow  loops  consisting  of  sections  with  liquid  flow  combined 
with  sections  having  liquid/vapor  flow  will  have  areas  with  different  speeds  of  propagation 
of  perturbations.  In  such  a  system,  a  perturbation  can  induce  velocity  and  pressure  oscillations 
which  may  diverge  or  reach  a  self-sustained  periodic  mode.  In  the  latter  case,  the  oscillations 
can  exceed  the  thermal  limits  of  the  system,  thereby  reaching  a  premature  CHF.  Premature 
CHF  has  also  been  demonstrated  in  parallel  channels  (see  Veziroglu  and  Lee,  Ref.  137)  that 
share  a  common  inlet  manifold  or  have  a  bypass  channel.  In  configurations  where  the  channels 
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have  similar  operating  conditions,  the  phase  of  oscillations  may  be  opposite,  depending  on 
the  number  of  channels  (Ref.  133).  Such  oscillations  are  called  channel-to-channel  instabilities. 
In  many  cases,  the  manifolded  cooling  approach  is  used  in  high-enthalpy  arc  facilities  and 
may  be  prone  to  these  types  of  instabilities  when  two-phase  flow  is  present.  Concepts  for 
future  arc  facility  nozzle  liners  include  thin  walls  with  structural  webbing  or  ribs  on  the  cooling 
channel  surface.  Such  designs  may  introduce  channel-to-channel  instabilities  between  the 
cooling  passages  on  the  nozzle  liners. 

Density-wave  oscillations  frequently  combine  with  other  phenomena  such  as  pressure- 
drop  oscillations  to  create  additional  instabilities.  For  instance,  Lowdermilk,  et  al.  (Ref.  101) 
and  later  Aladyev,  et  al.  (Ref.  45)  showed  that  pulsations  in  flowing  subcooled  and  saturated 
water  caused  by  an  upstream  compressible  volume  reduced  the  CHF  (from  that  of  a  stable 
system  with  no  pulsations)  as  much  as  80  percent.  Lowdermilk’s  team  found  that  the  flow 
instability  could  be  removed  by  throttling  the  flow  just  upstream  of  the  test  section  (Fig. 
34).  In  addition,  Gambill,  et  al.  (Ref.  115)  suggested  that  stability  is  improved  as  the  test 
section  pressure  drop  becomes  a  smaller  fraction  of  the  total  system  pressure  drop,  and 
restricted  flow  upstream  of  the  test  section  for  most  of  their  tests.*  It  appears  that  the  stability 
of  a  cooling  apparatus  which  involves  flow  boiling  can  be  improved  by  avoiding  compressible 
volumes  in  the  flow  loop  (or  throttling  the  flow  between  the  compressible  volume  and  the 
test  section)  and  forcing  the  test  section  pressure  drop  to  be  a  small  portion  of  the  total  pressure 
drop  of  the  system.  However,  as  pointed  out  by  Boyd  (Ref.  17),  in  cases  where  a  compressible 
volume  is  in  the  test  section,  or  in  very  long  test  sections,  no  amount  of  throttling  can  eliminate 
or  reduce  the  instability. 

Heat  Application/Orientation 

Although  not  typically  associated  with  an  actual  cooling  configuration,  the  method  of 
heat  application  when  resistive  (Joule)  heating  is  used  in  an  experimental  simulation  can  be 
important.  Ellion  (Ref.  138)  was  one  of  the  first  to  note  an  effect  of  alternating  current  (a-c) 
heating  on  an  experimental  boiling  apparatus.  Ellion  found  that  a  60-cycle/sec  a-c  power 
supply  caused  a  120-cycle/sec  growth  and  collapse  of  the  bubbles  in  flowing  subcooled  water 
in  an  annulus.  The  problem  was  rectified  by  switching  to  direct  current  (d-c)  power  with 
the  amplitude  of  the  voltage  ripple  less  than  1  percent  over  the  full  operational  range.  Gambill, 
et  al.  (Ref.  115)  suspected  a  similar  problem  when  they  used  a-c  power  to  electrically  heat 
a  tube  with  internally  flowing  subcooled  water.  Gambill’s  team  calculated  that  the  a-c  heating 
reduced  the  CHF  as  much  as  16  percent  over  that  for  d-c  heating.  Tippets  (Ref,  28),  who 
used  a-c  heating  of  stainless  steel  ribbons  to  study  subcooled  and  saturated  water  flow  boiling, 


*  Gambill’s  team  restricted  the  ratio  of  test  section  pressure  drop  to  the  total  system  pressure  drop  to  values  between  0.045  and  0.40. 
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found  that  the  CHF  for  0.006-in. -thick  heater  ribbons  averaged  20  percent  lower  than  those 
for  0.0 10-in. -thick  ribbons,  and  in  a  few  cases,  the  difference  was  substantially  larger.  Tippets 
theorized  that  the  reduction  was  due  in  part  to  the  lesser  thermal  time  constant  of  the  thinner 
ribbons  coupled  with  the  a-c  heating.  Although  Sadasivan  and  Lienhard  (Ref.  60)  showed 
a  definite  difference  between  a-c  and  d-c  heating  on  bubble  character,  they  found  little  effect 
on  CHF. 

Leung,  et  al.  (Ref.  44)  found  little  effect  (less  than  6-percent  difference)  of  direct  and 
indirect  heating  of  the  interior  wall  of  an  annulus  with  flowing  water  at  subcooled  inlet 
conditions.  As  mentioned  previously,  Collier  (Ref.  7)  points  out  that  the  outer  wall  in  a 
concentric  annulus  system  has  little  or  no  effect  on  the  CHF  for  the  inner  surface.  However, 
Tong  (Ref.  10)  suggests  that  the  presence  of  an  unheated  wall  in  the  proximity  of  a  critical 
point  can  adversely  affect  the  cooling  effectiveness  (termed  the  “unheated  wall  effect”). 
Tolubinsky,  et  al.  (Ref.  139)  found,  for  subcooled  and  saturated  water  flowing  in  annuli, 
that  higher  CHF  (as  much  as  30  percent)  existed  on  a  heated  inner  wall  than  on  a  heated 
outer  wall  at  comparable  conditions.  Yucel  and  Kakac  (Ref.  140)  found  that  the  heater 
orientation  had  an  effect  on  the  CHF  in  a  rectangular  channel  with  flowing  subeooled  water. 
Their  data  indicated  that  higher  CHF  is  achieved  with  the  heated  surface  (along  one  side 
of  the  rectangular  channel)  facing  up  rather  than  facing  down.  However,  they  attributed 
the  difference  not  to  the  presence  of  an  unheated  wall,  but  to  buoyancy  forces  moving  the 
vapor  away  from  the  lower  surface,  resulting  in  a  lower  surface  temperature  and  higher  CHF 
(as  opposed  to  a  higher  surface  temperature  and  lower  CHF  on  the  downward  facing  heater). 

Miscellaneous  Effects 

Several  other  parameters  have  been  shown  to  have  an  effect  on  the  cooling  process,  albeit 
only  a  small  number  of  experiments  have  been  performed  that  demonstrate  their  effect.  Tong 
(Ref.  10)  suggested  that  the  local  enthalpy  can  impair  the  CHF,  primarily  when  there  is 
significant  vapor  voidage  near  the  heated  wall.  Such  conditions  could  exist  for  subcooled* 
or  saturated  flow,  even  at  high  velocity.  The  local  enthalpy  near  the  wall,  even  at  upstream 
locations,  can  be  excessively  high  and  can  cause  burnout  when  significant  vapor  exists  in 
the  bubble  layer  near  the  wall. 

Johnson  (Ref.  141)  pointed  out  that  the  presence  of  an  electric  field  has  some  effect  in 
the  pure  convection  region,  little  effect  in  the  nucleate  boiling  region,  and  a  significant  effect 


*  Fiori  and  Bergles  (Ref.  20)  found,  contrary  to  prior  beliefs,  that  large  vapor  voidages  can  exist  even  in  subcooled  flow.  They 
noted  that  instantaneous  void  fractions  larger  than  50  percent  were  possible  at  high  heat  flux  (3.1  x  106  to  17.3  x  106  W/m2), 
high  exit  subcooling  (28°  to  56°C),  and  low  pressure  (less  than  6.2  bar). 
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at  CHF  and  beyond.  Markets  and  Durfee  (Ref.  142)  applied  a  3,000-v  a-c  voltage  across 
the  gap  of  an  annulus  with  the  inner  wall  heated  and  flowing  subcooled  water.  They  found 
that  the  CHF  was  augmented  as  much  as  40  percent  with  the  applied  electric  field  as  opposed 
to  no  field  present.  No  effect  of  the  applied  voltage  was  found  where  flow  instability  burnouts 
occurred. 

Another  heat  flux  augmentation  method  which  has  been  demonstrated  involves  the  use 
of  flow  and  ultrasonic  vibration.  Although  the  use  of  vibration  has  not  resulted  in  an  increase 
in  CHF  for  forced-flow  configurations  (Refs.  7  and  143),  an  enhancement  of  the  heat  transfer 
(increase  in  the  heat-transfer  coefficient)  in  the  single-phase  forced-convection  region  and 
boiling  onset  has  been  demonstrated  (see  Bergles,  Ref.  143). 

Boyd  (Ref.  1)  identified  two  other  heat  flux  enhancement  techniques  which  should  be 
included  here  for  completeness.  As  Boyd  noted,  Inoue  and  Bankoff  (Ref.  144)  found  that 
the  destabilization  of  film  boiling  (initially  subcooled)  by  the  introduction  of  a  pressure  wave 
increased  the  transient  heat  transfer  up  to  20  times  the  steady-state  values.  Boyd  also  noted 
that  Weede  and  Dhir  (Ref.  145)  showed  that  the  CHF  can  be  enhanced  and  its  location 
controlled  using  tangential  injection  of  a  fluid  into  the  mainstream.  This  technique  probably 
takes  advantage  of  the  induced  swirl  and  turbulence  whose  benefits  were  discussed  previously. 

3.2  SUMMARY  OF  PARAMETRIC  TRENDS  AND  IMPLICATIONS 

As  suggested  previously,  velocity  and  subcooling  have  a  primary  contribution  on  cooling 
effectiveness,  and  it  appears  that  a  highly  subcooled,  forced-convection  cooling  system  is 
most  attractive  for  high-enthalpy  arc  facility  component  cooling.  Figure  35  summarizes  the 
effect  of  a  majority  of  the  parameters  and  enhancement  techniques  previously  discussed  for 
a  subcooled,  forced-convection  cooling  environment.  Those  parameters  and  techniques  that 
tend  to  decrease  heat  transfer  or  increase  wall  temperature  in  the  pure  convection  (nonboiling) 
and  nucleate  boiling  regions,  respectively,  are  not  necessarily  detrimental  to  the  cooling  process 
unless  the  wall  temperature  increase  approaches  the  melt  temperature  of  the  material  or 
contributes  to  plasticity  to  the  point  of  structural  failure.  Probably  the  most  important 
contributions  to  cooling  effectiveness  and,  hence,  component  survival,  are  those  parameters 
and  enhancement  techniques  which  tend  to  increase  the  CHF.  A  number  of  parameters,  such 
as  pressure,  can  have  either  a  beneficial  or  a  detrimental  effect  on  CHF  and  must  be  optimized 
for  the  conditions  of  interest.  However,  system  limitations  and  configuration  requirements 
may  limit  the  degree  of  optimization  that  can  be  achieved. 

The  implications  of  requiring  a  subcooled,  forced-convection  cooling  configuration  are 
numerous  and  are  not  generally  conducive  to  experimental  evaluation.  As  previously  discussed, 
high  velocity,  high  subcooling,  and  high  pressure  tend  to  decrease  bubble  diameter  and  lifespan, 


25 


AEDC-TR-93-3 


where 

/rb/jtw  is  between  0  and  40 
Ren  is  between  104  and  5  X  106 
Pr  is  between  0.5  and  200  for  6-percent  accuracy 
Pr  is  between  200  and  2,000  for  10-percent  accuracy 
n  =  0.11  for  Tw  >  Tb 
n  =  0.25  for  Tw  <  Tb 
n  =  0  for  gases 
and  where 

f  =  1/(1. 82  logio  ReD  -  1.64)2  for  smooth  pipes 

or  from  Moody  chart  for  smooth  or  rough  pipes. 

Sleicher  and  Rouse  (Ref.  152)  recommended  a  simpler  expression  with  comparable  accuracy 
for  tube  flow, 

Nud  =  5  +  0.015  (ReD)a  (Pr)b  (6) 

where 

a  =  0.88  -  0.24  /  (4  +  Pr) 
b  =  (1/3)  4-  0.5  e(~~°-6p0 
and 

Pr  is  between  0.1  and  105 
ReD  is  between  104  and  106 

The  properties  in  ReD  are  evaluated  at  the  local  film  temperature,  Tf  =  (Tb  +  Tw)/2, 
and  the  properties  in  Pr  are  evaluated  at  the  local  wall  temperature. 
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Kays  and  Leung  (Ref.  153)  solved  the  turbulent-flow  energy  equation  for  constant  heat 
flux  in  annuli  over  wide  ranges  of  Reynolds  number,  Prandtl  number,  and  annulus  radius 
ratio.  Table  1  in  Ref.  153  or  Tables  13-3  to  13-5  in  Ref.  5  present  the  results  of  their 

computations. 

/ 

For  fully  developed  laminar  flow  in  a  circular  tube  with  uniform  heat  flux,  the  Nusselt 
number  and,  hence,  the  heat-transfer  coefficient  is  constant,  independent  of  Reynolds  number 
and  Prandtl  number, 

Nud  =  4.36  constant  heat  flux,  Pr  >  0.6  (7) 

Treatment  of  heat  transfer  in  the  entry  region  where  velocity  and  temperature  profiles 
are  not  fully  developed  is  discussed  in  Chapter  8  of  Ref.  154. 

Pressure  drop  is  also  of  interest  in  a  cooling  process,  and  the  pressure  gradient  has  even 
been  shown  to  have  an  effect  on  local  heat  transfer  (see  Chapter  6  of  Ref.  148).  Solution 
of  the  equations  of  motion  (using  semiempirical  shear  correlations  for  the  turbulent  case 
where  random  fluctuations  are  important)  provides  a  means  of  determining  the  pressure  drop 
in  a  flow  cooling  configuration  as  outlined  by  White  (Ref.  155).  Friction  factors  as  a  function 
of  Reynolds  number  for  flow  in  smooth  or  rough  circular  tubes,  annuli,  or  noncircular  ducts 
are  provided  by  the  Moody  chart  or  expressions  such  as  those  presented  in  Chapter  6  of  Ref. 
155. 

4.2  TRANSITION  TO  FLOW  BOILING,  FLOW  BOILING,  AND  CHF: 

THEORETICAL  MODELS 

The  theoretical  determination  of  heat  transfer,  pressure  drop,  and  other  flow  characteristics 
once  a  cooling  process  has  transitioned  into  boiling  is  one  of  the  greatest  challenges  facing 
researchers  in  two-phase  flow.  Three  currently  accepted  two-phase  flow  or  “mixture”  models 
have  been  developed  based  on  the  more  dominant  flow  regimes  in  two-phase  flow.  The  drift 
flux  model,  which  is  based  on  a  molecular  diffusion  model,  allows  for  motion  of  one  phase 
relative  to  the  other,  and  the  two  are  coupled  through  the  use  of  correlations  accounting 
for  interaction  effects.  The  correlations  are  difficult  to  obtain  because  they  require  a  knowledge 
of  the  average  effects  of  one  phase  moving  relative  to  the  other  (see,  for  example,  Ref.  156). 
Another  theoretical  model,  the  two-fluid  model  based  on  the  application  of  the  Navier-Stokes 
equations,  treats  the  phases  separately  and  requires  matching  conditions  at  both  solid 
boundaries  and  liquid-vapor  interfaces.  Difficulty  is  encountered  in  obtaining  local  and 
instantaneous  mass,  momentum,  and  energy  balance  equations  at  the  interfaces  (see,  for 
example,  Ref.  157).  The  drift  flux  model  is  typically  more  applicable  to  flows  where  the  one 
phase  is  dispersed  throughout  the  other.  The  two-fluid  model  is  generally  used  for  flows  where 
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the  two  phases  are  separated  into  layers,  but  has  been  accepted  for  use  with  dispersed  flows. 
If  the  motion  of  the  two  phases  can  be  assumed  to  be  equivalent  (i.e.,  no-slip  condition), 
a  simplified  model,  the  homogeneous  model  that  treats  the  mixture  as  a  single-phase  fluid, 
may  be  used.  The  transport  properties  required  for  the  model  are  difficult  to  determine  due 
to  the  lack  of  existing  relationships.  Although  extensive  progress  in  theoretical  modeling  has 
been  achieved  in  recent  years,  its  application  has  been  limited  because  of  the  difficulties 
described  above.  A  more  complete  treatment  of  these  models,  especially  the  two-fluid  model, 
is  included  in  Jones,  Ref.  158,  and  Drew,  Ref.  159. 

4.3  TRANSITION  TO  FLOW  BOILING:  CORRELATIONAL  APPROACHES 

Early  studies  in  flow  boiling  suggested  that  boiling  incipience  occurred  at  the  intersection 
of  the  pure  convection  and  fully  developed  nucleate  boiling  curves,  although  experimental 
data  indicated  a  transition  region  (McAdams,  et  al.,  Ref.  24).  Several  investigators  later  derived 
heat-transfer  expressions  for  the  transition  from  pure  convection  to  fully  developed  nucleate 
boiling  (termed  partial  nucleate  boiling).  Kutateladze  (Ref.  160)  and,  later,  Forster  and  Greif 
(Ref.  57)  suggested  correlations  based  on  the  pure  convection  heat  flux  at  the  incipience  of 
boiling  and  the  corresponding  value  of  heat  flux  in  pool  boiling.  Bergles  and  Rohsenow  (Ref. 
161),  however,  later  showed  that  the  transition  region  between  pure  convection  and  nucleate 
boiling  cannot  be  based  on  results  from  saturated  pool  boiling,  and  recommended  the  following 
relation  based  on  an  extension  of  the  fully  developed  nucleate  boiling  curve, 

q"/q"conv  =  [1  +  {(q'VcT'conv)  (1  -  q"Bi/q"B)}2]1/2  (8) 

where  q"B  is  calculated  from  a  fully  developed  nucleate  boiling  correlation  at  various  wall 
temperatures,  and  q"Bi  is  determined  from  an  extension  of  the  fully  developed  nucleate 
boiling  curve  to  the  temperature  corresponding  to  the  boiling  incipient  point  (see  Fig.  37). 
In  addition,  Bergles  and  Rohsenow  developed  a  relation  for  heat  flux  at  subcooled  boiling 
incipience  by  solving  a  bubble  growth  equation  graphically  to  obtain 

q"LB  =  15.60  pi-i5«  (Tw  -  Tsat)2-3/P0-0234  (9) 

4.4  FLOW  BOILING:  CORRELATIONAL  APPROACHES 

The  team  of  McAdams,  et  al.  (Ref.  24)  was  one  of  the  earliest  boiling  research  teams 
to  propose  an  empirical  correlation  for  fully  developed  nucleate  boiling  in  subcooled  forced 
flow.  They  found  for  flow  in  an  annulus  (heated  stainless  steel  inner  wall)  that  the  heat  flux 
during  boiling  was  independent  of  coolant  velocity  (1  to  36  ft/sec),  degree  of  inlet  subcooling 
(20°  to  150°F),  pressure  (30  to  90  psia),  and  equivalent  diameter  (0.17  to  0.48  in.).  They 
correlated  the  heat  flux  data  with  the  expression, 

q"  =  C '  ATsat3-86  (10) 


30 


AEDC-TR-93-3 


where  the  value  of  C'  decreased  from  0.19  to  0.074  as  the  concentration  of  dissolved  gas 
decreased  from  0.30  to  0.06  ml  of  air  at  standard  conditions  per  liter  of  water.  Jens  and 
Lottes  (Ref.  88)  analyzed  high-pressure,  high  heat  flux,  flow  boiling  data  obtained  by 
researchers  at  UCLA,  Purdue,  and  MIT,  and  derived  the  following  empirical  correlation, 

q"  =  (0.527  eU900  ATsat)4  (11) 

where  p  is  in  pounds  per  square  foot.  Jens  and  Lottes  found  that  the  wall  temperature  during 
boiling  was  independent  of  the  flow  velocity,  gas  content,  and  degree  of  subcooling  and  was 
dependent  only  on  heat  flux  and  water  pressure.  All  of  the  data  were  obtained  for  flow  inside 
round  or  square  tubes  made  of  stainless  steel  or  nickel. 

Giimour  (Ref.  162)  derived  a  correlation  based  on  dimensionless  groups  rather  than 
temperature  difference: 

(h/cf  G)  (c^/k/-6  (Qfa/p2)0.425  =  0.001  /(DG/fx^  (12) 


where 


G  =  (mv/A)  (er/fiv) 

and  where  p  is  in  pounds  per  square  foot.  Although  Giimour  substantiated  the  correlation 
with  data  from  various  experiments,,  including  several  in  which  water  was  used,  it  was  noted 
that  surface  material  had  a  significant  effect  that  was  not  accounted  for  in  the  correlation 
and  that  data  from  experiments  involving  boiling  from  wires  did  not  agree  well  with  the 
correlation. 

Later,  Levy  (Ref.  163)  derived  a  semiempirical  correlation  based  on  the  Forster  and  Zuber 
(Ref.  164)  bubble  growth  model  and  applicable  to  nucleate  boiling  with  flow  as  well  as  pool 
boiling,  arguing  that  the  fluid  velocity  normally  used  exerts  a  minor  effect  in  comparison 
with  the  local  turbulence  produced  by  the  bubble  motions.  For  subcooled  liquids,  Levy 
proposed  the  correlation, 

q"  =  [kfc(,e2/aTs  (e,  -  Gv)]  [1/BL]  [{ifg  +  cf  (Ts  -  Tb)}  /  ifg]  (ATsat)3  (13) 

where  Ts  and  Tb  are  in  °R,  a  in  Btu/ft2,  and  BL  is  an  empirical  constant  which  is  a  function 
of  qv  ifg  as  shown  in  Fig.  38.  Levy  found  that  the  above  correlation  agreed  reasonably  well 
with  data  from  McAdams,  et  al.,  and  the  UCLA  data  used  by  Jens  and  Lottes  discussed 
previously. 
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About  the  same  time  Levy  derived  his  correlation,  Forster  and  Greif  (Ref.  57)  theorized 
that  the  primary  boiling  heat-transfer  mechanism  is  the  exchange  between  the  vapor  and  liquid, 
and  that  microconvection  in  the  liquid  sublayer  beneath  the  bubbles,  bubbles  acting  as  artificial 
roughness,  and  latent  heat  transport  by  the  bubbles  contribute  very  little.  Based  on  this 
argument,  they  used  a  Reynold’s  analogy  to  derive  the  following  expression  for  general  boiling 
of  various  liquids  including  water: 

q"  =  Ksf  [aCfPfTs/ JifgQ v(j1/2]  [c(Tsa^/i  (ifgev)2]1/4  [ef/W]5/8  [^/k^  Ap2  (14) 

where  Ap  is  the  pressure  difference  in  pounds  per  square  foot  corresponding  to  the  superheat, 
ATsat.  Reasonable  agreement  was  achieved  when  the  correlation  was  compared  to 
experimental  data  at  pressures  up  to  50  atm  with  Ksf  equal  to  0.0012.  In  the  discussion  in 
Ref.  57,  the  authors  noted  that  the  constant  does  change  for  other  data. 

Shah  (Ref.  165)  presented  a  correlation  for  the  generalized  prediction  of  heat  transfer 
during  subcooled  boiling  in  annuli  based  on  earlier  work  involving  tubes  (Ref.  166).  His 
correlation  is 


q"  =  hf  (Tw  -  Tb)  +  hf  tyo  -  1)  (Tw  -  Tsat)  (15) 


where 

i/'o  =  1  +  46  B£5  for  B0  <  0.3  X  10~4 

t c  =  230  B°-5  for  B0  >  0.3  X  10~4 

B0  =  q"/Gifg 

hf  =  nonboiling  heat-transfer  coefficient 

Shah  found  that  the  expression  correlated  450  representative  data  points  with  a  mean 
deviation  of  8.9  percent.  The  range  of  test  conditions  was:  reduced  pressure  from  0.009  to 
0.89,  inlet  subcooling  from  2°  to  109°C  (3.6°  to  196°F),  mass  flux  from  278  to  7,774  kg/m2 
sec  (57  to  1,592  lbm/ft2  sec),  inner  tube  diameter  from  4.5  to  42.2  mm  (0.177  to  1.66  in.), 
and  annular  gap  from  1  to  6.6  mm  (0.039  to  0.26  in.).  Various  fluids,  including  water,  were 
evaluated  and  heating  was  from  the  inner,  outer,  or  both  walls  of  the  annuli. 

Probably  the  most  widely  accepted  and  verified  flow  boiling  correlation  is  based  on 
Rohsenow’s  (Ref.  167)  pool  boiling  correlation  (a  Nusselt  number  derived  from  the  fluid 
Prandtl  number  and  a  bubble  Reynolds  number)  and  proposed  by  Rohsenow  (Ref.  168)  in 
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1952.  The  heat-transfer  rate  associated  with  forced-convection  boiling  was  found  to  be 
accurately  predicted  by  addition  (or  superposition)  of  the  pool  boiling  effect  with  the  pure 
forced-convection  effect  discussed  above,  i.e., 

Q  =  Q  forced  convection  +  Q  pool  boiling  (16) 

where  q"forced  convection  is  determined  from  the  expressions  presented  in  the  previous  discussion 
of  pure  forced-convection  correlations  and  q"pooi  boiling  is  determined  from 

Cf  (Tw  -  Tsat)  /  ifgPrs  -  Csf  [{q"//ififg}  {gc  a/g  ( q{  -  pv)}1/2F  (17) 

where  the  properties  are  evaluated  at  Tsat  (except  as  noted),  r  equals  0.33,  s  equals  1  for  water 
cooling  and  1.7  for  other  fluids,  and  CSf  equals  0.006  for  nickel  and  brass  surfaces  and  0.013 
for  copper  and  platinum  surfaces.*  Vachon,  et  al.  (Ref.  169)  later  presented  new  values  of 
r  and  Csf  which  account  for  surface  preparation  technique  as  well  as  liquid-surface 
combination.  However,  their  analysis  applies  only  to  the  linear  portion  of  the  nucleate  boiling 
curve,  i.e.,  that  portion  which  excludes  partial  nucleate  boiling  and  approaching  CHF.  In 
addition,  the  effects  of  subcooling  and  pressure  have  not  been  fully  investigated  and  may 
affect  the  constants  presented  by  Vachon’s  team.  Tong  (Ref.  10)  pointed  out  that  the 
superposition  approach  may  be  questionable  in  light  of  the  experiments  of  Bergles  and 
Rohsenow  (Ref.  161)  which  demonstrated  the  differences  in  the  fluid  mechanics  of  flow  boiling 
and  pool  boiling.  Bergles  and  Rohsenow  recommended  a  less  restrictive  construction  approach 
for  predicting  forced-convection  boiling  when  the  net  vapor  generation  is  relatively  small, 
which  was  presented  previously  in  Eq.  (8).  Bjorge,  et  al.  (Ref.  170)  and  Stephan  and  Auracher 
(Ref.  171)  later  presented  variations  of  the  superposition  approach.  For  additional  approaches, 
Guglielmini,  et  al.  (Ref.  172)  provided  a  survey  of  flow  boiling  correlations  as  of  1980,  including 
boiling  incipience  and  partial  nucleate  boiling  expressions. 

Large  vapor  void  fractions  have  been  observed  in  low-pressure,  subcooled  flow,  possibly 
greater  than  20  percent  on  a  time-averaged  basis  and  exceeding  50  percent  on  a  local 
instantaneous  basis  (Ref.  20).  Such  void  fractions  contribute  to  the  accelerative  and  static 
head  components  of  the  pressure  gradient  along  the  heated  channel.  Flow  boiling  experiments 
where  pressure  drop  characteristics  were  evaluated  (e.g.,  Refs.  22  and  173)  generally  have 
shown  that  as  heat  flux  increases,  the  pressure  drop  decreases  until  boiling  begins,  after  which 
the  pressure  drop  increases  (Fig.  39).  Boyd  (Ref.  174),  however,  found  pressure  drops  for 
high-pressure,  moderately  subcooled  water  boiling  to  be  slightly  less  than  single-phase  pressure 
drops.  Boyd  noted  that  CHF  occurred  before  substantial  pressure  drops  occurred.  Bergles 


*  The  Rohsenow  correlation  is  valid  only  for  clean,  relatively  smooth  surfaces.  Rohsenow  (Ref.  167)  found  that  s  varied  erratically 
between  0.8  and  2.0  when  the  surface  was  dirty. 
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and  Dormer  (Ref.  173),  who  performed  experiments  in  subcooled  water  flow  boiling  at  pressure 
below  100  psia,  noted  that  little  success  had  been  reported  in  the  general  correlation  of  pressure 
drop,  and  presented  graphical  correlations  for  the  pressure  drops  encountered  in  their 
experiments.  About  the  same  time,  Staub  and  Walmet  (Ref.  175)  identified  the  two  regions 
before  and  after  the  point  of  significant  vapor  generation  (SNVG)  where  the  effects  of  void 
fraction  on  pressure  drop  are  different  in  subcooled  flow  boling.  Recently,  Hoffman  and 
Kline  (Ref.  176)  and  Lu  and  Jia  (Ref.  177)  have  presented  correlations  for  the  prediction 
of  void  fraction  and  the  subsequent  gravity,  acceleration,  and  friction  pressure  drop  before 
and  after  SNVG. 

4.5  CHF  CORRELATIONAL  APPROACHES 

No  theory  of  subcooled  burnout  has  yet  been  created  for  flow  boiling  (Ref.  33).  Complexity 
of  the  boiling  mechanisms  at  the  critical  point,  their  unsteady  nature,  and  the  small  scale 
of  the  heat-transfer  processes  have  hindered  the  development  of  theoretical  models  and 
understanding  of  the  phenomena  from  experimental  studies.  Moreover,  the  cooling  flow 
structure  and,  hence,  the  heat-transfer  processes  vary  depending  on  the  heating  condition, 
fluid  properties,  and  flow/surface  geometries  as  indicated  in  the  previous  discussion  of 
parametric  effects.  This  has  led  to  an  excessive  number  of  correlations  derived  for  the  prediciton 
of  CHF  in  flow  boiling.  Most  of  the  CHF  correlations,  not  unlike  those  for  fully  developed 
nucleate  boiling,  were  developed  either  by  the  analytical  (mechanism-based)  approach, 
dimensional  analysis/similitude-based  approach,  or  optimization  of  test  data  (empirical-based). 

4.5.1  Mechanism-Based  CHF  Correlations 

In  an  attempt  to  derive  analytical  relationships  among  the  macroscopic  heat-transfer  and 
flow-field  characteristics  in  a  boiling  environment,  several  models  describing  the  mechanisms 
on  a  microscopic  scale  have  been  proposed.  The  mechanisms  appear  to  vary  with  the  degree 
of  subcooling,  and  several  current  approaches  divide  the  subcooled  boiling  regions  into  three 
areas:  (1)  high  subcooling  and  mass  velocity,  (2)  moderate  subcooling,  and  (3)  low  subcooling, 
velocity,  and  pressure  (Ref.  8).  At  high  subcooling  and  mass  velocity,  it  has  been  proposed 
that  local  wall  overheat  occurs  when  high  heat  flux  leads  to  dryout  and  a  resulting  hot  spot 
beneath  a  bubble;  however,  there  is  little  experimental  support,  and  no  correlations  have 
been  developed  based  on  this  mechanism  (Ref.  8).  At  moderate  subcoolings,  it  has  been 
proposed  that  the  liquid  flow  stagnates  near  the  heated  wall  when  the  boundary  layer  separates. 
At  high  heat  flux,  the  stagnant  liquid  evaporates,  resulting  in  vapor  blanketing  near  the  wall 
and  subsequent  burnout.  Tong  (Ref.  178)  presented  a  correlation  based  on  this  mechanism 
which  agreed  reasonably  well  with  data  from  subcooled,  high-pressure  water  flow  in  tubes, 
annuli,  and  rod  bundles.  Purcupile  and  Gouse  (Ref.  179)  presented  a  similarly  based  correlation 
in  nondimensional  form.  Other  studies  associated  with  this  proposed  mechanism  include  the 
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work  of  Kutateladze  and  Leont’ev  (Ref.  180)  and  Hancox  and  Nicoll  (Ref.  181).  At  low 
subcooling,  velocity,  and  pressure,  Fiori  and  Bergles  (Ref.  20)  suggested  that  dry-out  of  the 
microlayer  follows  when  inadequate  quenching  of  the  surface  by  adjacent  liquid  occurs  after 
the  passage  of  a  bubble.  Subsequent  burnout  occurs  when  high  surface  temperatures  prevent 
rewetting  of  the  surface.  Because  of  difficulties  in  defining  the  flow  structure,  no  correlation 
could  be  obtained  by  the  authors.  Lee  and  Mudawwar  (Refs.  182  and  183)  proposed  a  similar 
mechanism  by  suggesting  that  CHF  occurs  as  a  result  of  the  Helmholtz  instability  at  the 
microlayer-vapor  interface,  leading  to  dry-out  of  the  microlayer.  They  demonstrated  that 
their  correlation  agreed  with  numerous  subcooled  (over  a  wide  range  of  inlet  subcoolings) 
water  data  points  with  a  mean  deviation  of  approximately  11  percent.  Katto  (Refs.  184  and 
185)  presented  a  correlation  based  on  this  same  sublayer  dry-out  mechanism  and  suggested 
that  CHF  of  subcooled  flow  boiling  has  strong  similarities  to  pool  boiling  CHF.  Earlier, 
Katto  first  presented  generalized  correlations  for  CHF  in  subcooled  and  saturated  water  flow 
in  tubes  (Refs.  186  and  187)  and  annuli  (Ref.  188)  based  on  the  existence  of  four  characteristic 
regimes  of  CHF  and  the  following  dimensionless  groups: 

q"cHF  /  {Gifg  [1  +  K  (Aij/ifg)]}  =  f  (ev/Qt,oQt/&UL/D)  (18) 

where  K  is  determined  analytically  to  account  for  inlet  subcooling.  A  reasonable 
correspondence  between  the  correlations  and  experimental  data  was  found;  however,  later 
studies  (e.g.,  Ref.  189)  identified  problems  with  applying  the  correlations  to  the  various  regimes. 
In  addition,  difficulty  was  encountered  in  deriving  an  expression  for  CHF  at  very  high  mass 
velocity  and  pressure  conditions.  Nishikawa,  et  al.  (Ref.  190)  later  correlated  high-pressure 
tube  data  at  various  Freon  mass  flows  with  Katto’ s  proposed  dimensionless  groups  using 
the  boiling  length  as  the  characteristic  length.  Their  proposed  expressions  correlated 
approximately  150  data  points  to  within  ±10  percent. 

Other  proposed  mechanisms  of  general  subcooled  flow-boiling  CHF  include  liquid  flow 
blockage  models  and  vapor  removal  limit/bubble  crowding  models  (Ref.  184).  The  liquid 
flow  blockage  models  assume  that  the  flow  of  liquid  normal  to  the  heated  surface  is  blocked 
by  the  flow  of  vapor  along  the  surface.  The  mechanisms  proposed  by  Bergel’son  (Ref.  191) 
and  Smogalev  (Ref.  192)  involve  some  form  of  liquid  flow  blockage.  Griffith  (Ref.  193) 
presented  a  dimensionless  correlation  based  on  a  critical  packing  of  bubbles  on  the  surface. 
Griffith  proposed  that  a  saturated  pool  boiling  bubble  growth  term  be  multiplied  by  a 
correction  factor  which  included  velocity  and  subcooling  effects.  Approximately  94  percent 
of  his  data  were  correlated  by  the  expression  to  within  ±  33  percent.  Hughes  (Ref.  21)  suggested 
that  CHF  occurred  when  a  maximum  number  of  nucleation  sites  available  for  a  surface  was 
exceeded  (i.e.,  bubble  crowding).  For  Freon  flow  over  a  straight  surface,  he  proposed  the 
correlation: 
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Q"chf  —  [2dsGv  (cfATsub  +  2ifg)]  /  3K  {4/9ira  (evATbuikds/ GfATsat)2  + 

W3  (evifgds/kfATsat)2  +  (1/B)  to  [(U0  +  Bds)  /  U0J}  (19) 

where 

ds  =  6.12  W  Dh  /  Re7/8  (W  =  30  in  Ref.  21) 

B  -  (12/11)  (CDW/e(d2) 

U0  =  (3/4t to)  (k2  ATsat2  /  G2  ifg2d) 

and  q"cHF  is  in  Btu/in.2  sec,  q  is  in  lbm/in.3,  k  is  in  Btu/in.  sec  °F,  fi  is  in  lbm/in.  sec, 
Dh  and  d  are  in  in.,  and  a  is  in  in.Vsec.  The  term  K  was  determined  experimentally  as 

K  =  jSfs  Vm1/2/ATsub  (/3fs  is  approximately  equal  to  7.94) 

Although  not  given  by  Hughes,  the  drag  coefficient  for  a  bubble  is  approximated  by 

CDSphere  *  +  +  ^  ^  ^  ChaP'  ^ 

and,  for  a  straight  test  section,  the  bubble  diameter  for  the  preceding  equations  is  approximated 
by 

d  -  ds 

Hughes  also  presented  somewhat  cumbersome  algebraic  expressions  for  CHF  on  convex 
and  concave  surfaces  (see  Ref.  21).  Hebei,  et  al.  (Ref.  194)  similarly  suggested  that  CHF 
occurs  when  a  limiting  vapor  removal  rate  is  reached.  Other  proposed  mechanisms  included 
in  this  group  are  those  suggested  by  Yagov  and  Puzin  (Ref.  195)  and  later,  Weisman  and 
Ileslamlou  (Ref.  196).  Weisman  and  Ileslamlou  (see  also  Yang  and  Weisman,  Ref.  197),  by 
extending  the  earlier  work  reported  in  Ref.  198  based  on  limited  turbulent  interchange  at 
the  outer  edge  of  a  bubbly  layer  near  a  heated  surface,  proposed  the  following  CHF  correlation 
for  use  in  highly  subcooled  flow: 

q"cHF  =  Gi/4  [if  (1  -  x2)  +  iv  (x2)  -  i]  (20) 
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where 


r P  =  [1/(2)1/27t]  exp  [-  1/2  (vn/a;)2]  -  1/2  {vn/o'v)  erfc  (vn  /  (2)1/2  a'v) 
ib  =  0.462  (K)°-6  (Re)-o.i  (Db/D)°-6  [1  +  a  (Qe  -  Gv)  /  Gv] 
i  is  the  average  fluid  enthalpy 

x2  is  the  average  quality  in  the  boundary  layer  corresponding  to  a  void  fraction  of  0.82 
K  =  2.4  (experimentally  determined) 

Db  =  [0.015  (uDh/ rw),/2]  {1  +  0.1  (g/gc)  [(q,~  Qv)  /  rw]  Dh}~^ 

buoyancy  correction  term  »  1  for  high  mass  velocity 
vn  =  q"ifg  (x2)  /  [{if  (1  -  x2)  +  iv  (x2)  -  i]  (Gvifg)} 

Uv  =  (G/ed  ib 

a  -  0.87  (1.36  -  VSL)  +  0.135  0.5  <  VSL  <  1.36  m/sec 

a  =  0.135  G  <  9.7  x  106  kg/m2  hr,  Vsl  2:  1.36  m/sec 

a  =  0.135  (G/9.7  x  106)-°-3  G  >  9.7  x  106  kg/m2hr 

where  VSl  is  the  superficial  velocity  of  two-phase  mixture  in  m/sec.  The  correlation  requires 
iteration  on  q"cHF  since  Vy  is  dependent  on  q".  Andreyev,  et  al.  (Ref.  199)  suggested  that 
CHF  that  occurs  in  the  transition  from  bubbly  to  dispersed-annular  flow  can  be  based  on 
a  transition  void  fraction  and  recommended  an  equation  for  this  specific  void  fraction  that 
agreed  reasonably  well  with  experimental  data  for  subcooled  water  flow  boiling.  Povarnin 
(Ref.  200)  assumed  a  concept  of  corresponding  states  and  derived  a  correlation  that  agreed 
reasonably  well  with  subcooled  CHF  data  for  various  fluids  including  water.  No  further  work 
has  been  performed  in  recent  years  to  substantiate  Povarnin’s  approach.  Bergles  (Ref.  8) 
identified  several  other  proposed  mechanism-based  correlations  for  subcooled  flow  boiling 
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CHF.  Specifically,  Chang  (Ref.  37)  developed  a  correlation  for  CHF  based  on  the  idea  that 
large  bubbles  are  broken  up  in  a  hydrodynamic  limit,  thereby  reducing  the  heat-transfer 
coefficient.  Thorgerson,  et  al.,  (Ref.  201)  proposed  a  correlation  which  allowed  a  burnout 
condition  as  a  result  of  a  critical  friction  factor.  His  correlation  falls  somewhat  in  the  class 
of  expressions  based  on  boundary-layer  separation.  Both  of  these  correlations  have  met  with 
little  success  in  experimental  verification  or  general  acceptance. 

Probably  the  most  widely  verified  mechanism-based  correlation  uses  the  superposition 
approach  originally  proposed  by  Rohsenow  for  fully  developed  nucleate  flow  boiling  and 
extended  to  CHF  prediction  by  Gambill  (Ref.  202): 

q"cHF  =  [q"cHF]pool  boiling  +  [q"cHFlforced  convection  (21) 


where 

[q  CHFlpool  boiling  =  Klfg  Q\  [°Sc§  (@P  ~  Qv)/Qv\^^4  [1  F  (Qf/Qv)^'^  (CpATsub/25ifg)] 
in  which 


K  =  0.12  —  0.17,  cp  is  evaluated  at  T  =  Tsat  -  (ATsub/2) 

[q"cHF]  forced  convection  —  hconv  1(Tw)chf  -  Tb] 

where  (TW)CHF  is  evaluated  with  Bernath’s  (Ref.  203)  generalized  plot  of  wall  superheat  at 
burnout  (Fig.  40).  Gambill  found  that  the  correlation  predicted  burnout  for  a  large  number 
of  experimental  data  points  within  a  maximum  deviation  of  40  percent  for  various  fluids 
and  17.8  percent  for  water.  The  data  compared  were  for  numerous  fluids  and  flow 
configurations  with  a  range  of  conditions  of  velocity:  0  to  174  ft/sec,  pressure:  4  to  3,000 
psia,  subcooling  (various  locations):  0°  to  506°F,  acceleration:  1  to  57,000  g,  CHF:  0.1  x 
106  to  37.4  x  106  Btu/hr  ft2.  Levy  (Ref.  204)  presented  a  similar  superposition  correlation 
for  CHF: 


qwcHF  -  q"p  +  q"c  +  qwF  (22) 

where 

q"P  =  0.131  ifg  pv  [ag2c  (q(  -  qv)/qW/4 
q"c  =  0.696  (kfe<cf)1/2  [(fir  -  Q,V^l/4  [o&Qt  ~  6v)/e?]1/8  ATsub 
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and 


Q"f  =  hf  (ATsat)  +  hpATsub 

where  hf  is  determined  from  a  pure  forced-convection  Nusseit  number, 

hf  =  Nu  k/D 

and  (ATsat)  is  determined  by  trial  and  error  from  a  fully  developed  nucleate  boiling  relation 
such  as  Jens  and  Lottes  [Eq.  (11)]  at  q"cHF- 

The  list  of  mechanism-based  CHF  correlations  is  by  no  means  limited  to  those  presented 
here.  Those  for  which  the  equations  are  presented  appear  to  be  generally  applicable  to  highly 
subcooled,  water  flow  boiling  which  will  be  encountered  in  the  cooling  of  high-enthalpy  arc 
facility  components.  A  statement  made  by  Bergles  in  Ref.  8,  which  is  applicable  here,  is  that 
there  are  enough  “adjustable”  constants  in  any  predictive  equation  to  permit  an  acceptable 
correlation  of  data. 

4.5.2  Dimensional  Analysis/Similitude-Based  CHF  Correlations 

Although  dimensional  analyses  have  been  performed  by  Barnett  (Ref.  205),  Zenkevich 
(Ref.  206),  and  Kampfenkel  (Ref.  207),  Glushchenko  (Ref.  86)  used  partial  modeling  to  obtain 
a  correlation  for  highly  subcooled  (greater  than  25 °C)  water  flow  in  tubes  and  annuli  based 
on  the  following  dimensionless  groups: 


Kl  -  q"cHF/ifg  6v  Vavg 


K-2  —  Cp  ATsub  Qf/ifg  Qv 


K-3  —  k|>/  VavgDecp  q  p 
K4  =  ifg/ cpATsat 


Glushchenko  proposed  the  correlation 


Kj  =  18.25  K®'35  K3'5  K.42  (23) 

which  agreed  reasonably  well  with  experimental  data  (80  percent  of  nearly  200  selected  data 
points  correlated  within  25  percent).  The  data  covered  the  ranges:  mass  velocity  from  500 
to  40,000  kg/m2  sec  (100  to  8,200  lbm/ft2  sec),  pressure  from  4.9  x  105  to  197  x  105  N/m2 
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(72  to  2,900  psia),  subcooling  (location  not  discussed)  from  25°  to  250°C  (45°  to  450°F), 
diameter  from  2  to  12  mm  (0.08  to  0.5  in.),  and  L/D  from  10-20  to  60-120.  Ornatskii,  et 
al.  (Ref.  208),  later  extended  the  approach  to  low  subcooling  and  moderate  quality  water  flow. 

Probably  the  most  successful  approach  to  similitude  analysis  of  the  CHF  phenomena 
was  performed  by  Ahmad  (Ref.  70).  The  primary  impetus  for  such  analyses  was  the  need 
for  fluid-to-fluid  modeling  capability.  Refrigerants  such  as  Freon  have  been  used  in  previous 
experiments  (Refs.  66-69)  to  simulate  water  cooling  at  high  heat  flux  conditions,  since  the 
CHF  for  Freon  can  be  achieved  at  much  lower  heat  flux.  Scaling  methods  were  necessary 
to  provide  a  means  of  relating  the  Freon  test  results  to  water  for  high  heat  flux/water  cooling 
predictions.  Ahmad  used  classical  dimensional  analysis  to  identify  12  independent  groups, 
of  which  6  were  eliminated  through  inductive  arguments.  Of  the  remaining  six  terms,  3  terms 
(subcooling  number  Ai/ifg,  liquid/vapor  density  ratio  q/qv,  and  L/D)  can  be  independently 
satisfied  in  a  controlled  experiment  (controlled  inlet  temperature,  pressure,  and  geometry). 
The  remaining  3  terms  (Reynolds  number,  Weber-Reynolds  number,  and  liquid/vapor  viscosity 
ratio)  are  distorted  through  the  liquid  and  vapor  viscosities  when  the  above  conditions  are 
fixed.  Ahmad  solved  this  problem  of  multiple  distortion  by  expressing  the  remaining  3  terms 
as  a  CHF  modeling  parameter  \Pchf,  and  using  a  compensated  distortion  model  to  relate 
the  terms  through  two  empirical  exponents.  The  final  form  of  the  modeling  parameter  is 
written  as 


'/'CHF  -  KGD/fid  (w  /  oDQ()2/i  0c.//c)1/5]  (24) 

From  the  Buckingham  Pi  theorem  the  dependent  parameter,  boiling  number 
(q"cHF/Gifg),  is  written  as  follows: 


Q"cHF/Gifg  —  f  (i/chf>  Ai/ifg,  Q(/qv,  L/D)  (25) 

In  addition  to  the  independent  parameters  identified  above,  Ahmad  included  the  ratio 
of  the  heated  equivalent  diameter  (4  x  flow  area/heated  perimeter)  to  the  hydraulic  diameter 
(4  x  flow  area/ wetted  perimeter),  Dhe/D,  for  more  complex  geometries.  Ahmad 
demonstrated  reasonable  scaling  of  subcooled  Freon  CHF  data  from  Coffield,  et  al.  (Ref. 

209)  to  water  CHF  using  the  above  approach.  Analysis  of  complex  geometries  by  this  method 
is  limited  by  the  modeling  parameter  in  which  the  empirical  exponents  are  highly  sensitive 
to  geometric  parameters.  Use  of  the  correlation  for  situations  other  than  extrapolation  of 
refrigerant  CHF  data  to  water  CHF  has  been  very  limited  (see,  for  example,  Katto,  Ref. 

210) . 
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4.5.3  Empirical-Based  CHF  Correlations 

Probably  the  most  common  approach  to  obtaining  a  CHF  correlation  is  the  curvefitting 
of  experimental  test  data,  although  some  confusion  exists  as  to  whether  or  not  a  correlation 
is  considered  empirical  when  some  or  all  of  the  parameters  used  in  the  correlation  are  obtained 
from  dimensional  analysis  or  limited  boiling  mechanism  analyses.  Because  most  correlations 
rely  on  some  form  of  empiricism,  this  is  probably  a  moot  point,  and  some  overlap  of  the 
mechanism-based,  dimensional  analysis/similitude-based,  and  empirical-based  approaches 
exists.  A  word  of  caution  about  the  use  of  empirical-based  correlations:  the  expressions 
typically  have  been  derived  for  a  specific  range  of  experimental  conditions,  and  extrapolation 
of  a  correlation  outside  the  specific  range  can  lead  to  very  large  errors. 

The  earliest  empirical  CHF  correlations  appeared  in  the  late  1940’s  into  the  early  1950’s, 
most  notably  those  of  McAdams,  et  al.  (Ref.  24),  Jens  and  Lottes  (Ref.  88),  Gunther  (Ref. 
27),  Buchberg,  et  al.  (Ref.  100),  and  McGill  and  Sibbitt  (Ref.  211).  The  correlations  proposed 
by  McAdams’  team  and  Gunther  are  limited  to  low  pressure  and  only  McGill  and  Sibbitt 
evaluated  subcoolings  greater  than  250°F.  All  of  the  correlations  are  limited  to  low  coolant 
velocity  (less  than  40  ft/sec  or  a  mass  velocity  less  than  2,500  lbm/ft  sec).  Only  the  Jens 
and  Lottes  correlation  included  a  pressure  effect,  and  none  included  geometry  effects.  The 
test  configurations  were  typically  tubes  except  those  used  by  McAdams’  team  (annulus)  and 
Gunther  (metal  strip  in  a  rectangular  channel). 

Weatherhead  (Ref.  43)  later  modified  the  Jens  and  Lottes  correlation  to  include  the  effects 
of  geometry  and  latent  heat  of  vaporization.  Weatherhead  subdivided  subcooled  flow  boiling 
into  four  classifications  and  accounted  for  variations  in  boiling  characteristics  at  low  and 
high  pressure,  low  and  high  mass  velocity,  small  and  large  geometries,  and  low  and  high 
inlet  subcooling.  The  correlation  at  high  subcooling,  however,  could  only  be  validated  at 
low  L/D,  and  Tong  (Ref.  10)  later  found  that  Weatherhead’s  correlation  did  not  improve 
the  accuracy  of  prediction  when  other  data  were  evaluated. 

A  number  of  empirical  correlations  appeared  in  the  late  1950’s  into  the  mid-  1960’s,  many 
of  which  are  summarized  in  Ref.  10  (see  also  Refs.  212  and  213).  Most  notable  are  correlations 
derived  by  Bernath  (Ref.  203),  Janssen,  et  al.  (Refs.  214  and  215),  Labuntsov  (Ref.  216), 
Van  Huff  and  Rousar  (Ref.  29),  and  the  United  Kingdom  Atomic  Energy  Establishment  at 
Winfrith,  AEEW  (Refs.  83,  91,  217,  and  218).  Bernath  developed  an  empirical  correlation 
based  on  the  premise  that  at  CHF  the  two-phase  flow  near  the  heated  surface  is  highly  turbulent 
and,  therefore,  well  mixed.  Assuming  the  convective  heat  transfer  through  this  mixture, 
Bernath  empirically  formulated  expressions  for  the  wall  superheat  and  the  heat-transfer 
coefficient  at  the  critical  point  (Ref.  219).  He  later  (Ref.  203)  refined  the  correlation  for  an 
extended  range  of  geometric  variables,  and  for  water  is  given  as 
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Q"CHF  =  hCHF  (Tw.CHF  “  Tb) 


(26) 


where 


Tw,chf  =  57  fti  p  -  54  [p/(p  +  15)]  -  V/4 

=  19,602  [De/(De  4-  Dj)]  +  (slope)  V 

slope  =  86.4/D°-6  for  De  <  0.1  ft 

slope  —  162  +  (18/De)  for  De  >  0.1  ft 

where  hCHF*  is  in  Btu/ft2hr°C,  TVV)chf  and  Tb  are  in  degrees  Centigrade,  and  Dj  is  the  heated 
perimeter  (inner  diameter)  in  feet  divided  by  pi.  Bernath  found  that  the  correlation  predicted 
the  CHF  for  nearly  250  forced  water  flow  experiment  data  points  within  16  percent.  The 
range  of  parameters  for  the  correlation  is:  pressure  from  23  to  3,000  psia,  velocity  from  4 
to  54  ft/sec,  subcooling  (various  locations)  from  0°  to  615°F,  and  hydraulic  diameter  from 
0.143  to  0.66  in.  Multiple  configurations  including  tubes,  annuli,  ducts,  and  ribbons  were 
used  in  the  experiments. 

Janssen  and  others  (Refs.  214  and  215)  developed  rather  cumbersome  design  correlations 
(presented  in  Ref.  10)  for  CHF  in  subcooled  and  saturated  water  flow  in  tubes  and  annuli 
(inner  wall  heated).  Their  correlations  are  valid  over  the  parameter  ranges:  pressure  from 
600  to  1,450  psia,  mass  velocity  from  0.2  to  6.2  lbm/ft2  hr,  and  hydraulic  diameters  from 
0.25  to  1.25  in.  Labuntsov  (Ref.  216)  assumed  (in  the  first  approximation)  that  the  CHF 
is  dependent  on  pressure,  velocity,  and  subcooling,  and  developed  a  CHF  correlation  for 
water  flow  at  pressures  up  to  200  atm.  Van  Huff  and  Rousar  (Ref.  29)  compiled  CHF  data 
on  23  different  fluids,  including  water,  during  the  study  of  heat  flux  limits  of  storable 
propellants.  They  correlated  the  higher  VATSUb  water  data  as  shown  in  Fig.  7  with  the 
expression 


Qwchf  -  5.1  +  0.000860  VATsub  (27) 

where  q"cHF  is  in  Btu/in.2sec.  Considerable  scatter  was  noted  at  low  VATsub  (less  than  10,000 
ft  °F/sec),  and  a  significant  pressure  effect  was  observed  for  data  obtained  in  swirl 


*  Bernath’s  original  work  computed  hCHF  in  units  of  P.c.u./hr  ft2oC,  where  P.c.u.  is  a  pound-centigrade  heat  unit.  His 
equations  presented  here  have  been  converted  to  Btu  rather  than  P.c.u. 
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flow.  The  correlation  is  valid  over  the  range  of  parameters:  pressure  from  10  to  2,000  psia, 
velocity  from  7.5  to  205  ft/sec,  and  bulk  fluid  temperature  (various  locations)  from  76°  to 
470°F. 

Lee  and  Obertelli  (Ref.  83)  at  AEEW  in  the  UK  satisfactorily  correlated  early  AEEW 
subcooled  flow  boiling  CHF  data  obtained  at  1 ,000  psia  using  a  modified  version  of  an  widely 
used  correlation  derived  by  DeBortoli,  et  al.  (Ref.  220)  for  higher  pressure.  Lee  and  Obertelli 
also  found  that  a  more  general  correlation  (based  on  a  local  conditions  hypothesis)  under 
consideration  by  AEEW  colleagues  at  about  the  same  period  of  time  (see  Barnett,  Ref.  221) 
provided  excellent  accuracy  in  the  prediction  of  over  600  CHF  data  points  obtained  at  AEEW 
over  a  range  of  pressure  from  560  to  1,600  psia.  Macbeth  (Ref.  217)  and,  later,  Thompson 
and  Macbeth  (Ref.  218)  correlated  nearly  4,400  CHF  data  points  (with  a  root-mean-square 
error  of  7.5  percent)  obtained  for  water  flow  boiling  in  round  tubes  with  expressions  based 
on  Barnett’s  local  condition  hypothesis.  The  hypothesis  assumes  that  the  CHF  is  only  a  function 
of  the  mass  quality  at  the  point  of  overheating.  Macbeth  showed  that  a  nearly  linear 
relationship  exists  at  low  mass  velocity;  therefore,  he  divided  the  boiling  map  into  a  low- 
velocity  regime  and  a  high-velocity  regime.  Figure  41  presents  the  boundary  limits  of  the 
two  regimes  as  identified  by  Macbeth.  The  CHF  correlations  for  the  two  regimes  as  given 
by  Thompson  and  Macbeth  are 

High-Velocity  Regime 

Q"chf  =  {[A'  4-  0.25D  (G  X  10-6)  Aid  /  (C'  +  L)}  X  106  (28) 

where  A'  and  C'  are  given  in  Fig.  42, 

Low-Velocity  Regime 

Q"chf  =  {[(G  X  10-«)  (ifg  +  Ai;)]  /[1 58  D01  (G  x  10-6)0-49  +  4L/D]}  X  106  (29) 

The  terms  L  and  D  are  expressed  in  inches.  Barnett  (Ref.  91)  later  modified  the  basic  form 
of  these  equations  to  obtain  a  correlation  more  applicable  to  annulus  data 

q"cHF  =  [(A  +  BAiO  /  (C  +  L)]  x  10«  (30) 


where 

A  =  67.45  Dhe0'68  (G  X  lO”6)0-192  {1  -  0.744  exp  [-6.512  DHY  (G  x  10~6)]} 
B  =  0.2587  Dhel261  (G  X  IQ-6)0-817 
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C  =  185.0  DHyL415  (G  X  10-6)0-212 
Dhe  =  (Ph  -  Di)/D[ 

Dry  -  Do  -  Di 

The  terms  L,  D0,  and  Dj  are  expressed  in  inches.  With  the  correlation  above,  Barnett 
correlated  724  annulus  CHF  data  points  with  a  root-mean-square  error  of  5.9  percent.  Limits 
of  the  correlation  are  Dj  from  0.375  to  3.798  in.,  Do  from  0.551  to  4.006  in.,  L  from  24 
to  108  in.,  G  x  10~6  from  0.14  to  6.2  lbm/hr  ft2,  and  Aii  from  0  to  412  Btu/lbm. 

A  number  of  empirical  CHF  correlations  developed  for  subcooled  flow  boiling  have 
appeared  since  the  1960’s.  Tolubinsky,  et  al.  (Ref.  139),  at  the  Institute  of  Engineering 
Thermophysics,  USSR,  proposed  a  design  equation  for  determining  CHF  in  annular  channels 
with  the  inner  wall  heated,  noting  that  because  the  CHF  is  affected  by  the  heated  surface 
shape  and  dimensions,  the  CHF  data  obtained  from  tube  flow  cannot  be  applied  to  annular 
flow.  Their  expression  correlated  more  than  90  percent  of  approximately  400  experimental 
data  points  to  within  ±  25  percent.  The  correlation,  however,  requires  knowledge  of  the 
CHF  during  pool  boiling  conditions  for  the  annular  configuration  of  interest,  which  may 
hinder  its  usefulness.  More  recent  work  performed  in  the  USSR  regarding  CHF  in  flow  boiling 
correlations  includes  the  studies  by  Levitan  and  others  (Refs.  222  and  223)  and  the  Scientific 
Council  of  the  Academy  of  Sciences  of  the  USSR  (Ref.  224).  Levitan  and  his  team  proposed 
the  following  correlation  for  subcooled  and  low-quality  water  flow  in  tubes: 

Q"chf  =  [10.3  -  7.8  (p/98)  +  1.6  (p/98)2]  e"1-5*  (G/l.OOO)1-2!0-25^98)  -  «  -  x>(31) 

where  q"cHF  is  in  MW/m2,  p  is  in  bars,  G  is  in  kg/m2  sec,  and  x  is  the  relative  quality  at 
CHF.  The  correlation  is  good  for  water  flow  in  circular  tubes  with  a  diameter  within  the 
range  of  4  to  16  mm  (0.16  to  0.63  in.).  Levitan  and  his  team,  along  with  several  other 
researchers,  later  assisted  the  Scientific  Council  of  the  Academy  of  Sciences  of  the  USSR 
(Ref.  224)  in  compiling  tabular  data  for  predicting  CHF  in  uniformly  heated,  8-mm  (0.31-in.)- 
diam  tubes.  An  approximation  formula  is  provided  for  calculating  CHF  for  tubes  whose 
diameter  is  different  than  8  mm: 


Q"chf  =  q"cHF,8  (8/D)0'5  (32) 

where  q"CHF  and  q^CHF.s  are  in  MW/m2,  q"cHF,8  is  the  value  of  q"cHF  in  the  8-mm-diam 
tube,  and  D  is  in  millimeters.  The  tabular  data  are  quoted  to  have  a  root-mean-square  error 
of  10  percent.  The  tabular  correlation  is  limited  to  the  following  range  of  conditions:  pressure 
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from  29.5  to  196  bar  (435  to  2,880  psia),  mass  velocity  from  750  to  5,000  kg/m2  sec  (154 
to  1,024  lbm/ft2  sec),  subcooling  (location  not  discussed)  from  0°  to  75°C  (0°  to  135°F), 
diameter  from  4  to  16  mm  (0.16  to  0.63  in.),  and  length-to-diameter  ratio  of  20  or  greater. 
Doroschuk,  et  al.  (Ref.  222)  modified  Eq.  (31)  to  include  a  variable  critical  pressure  such 
that  the  correlation  could  be  used  for  other  fluids. 

Bowring  (Ref.  225)  built  upon  the  work  of  Thompson  and  Macbeth  (Ref.  218)  and 
developed  a  CHF  correlation  using  four  basic  variables  as  functions  of  pressure.  He  showed 
that  the  correlation  predicted  approximately  3,800  data  points  with  a  root-mean-square  error 
of  6.96  percent.  At  about  the  same  time,  Becker,  et  al.  (Ref.  226)  proposed  a  CHF  correlation 
for  round  tubes  (their  data  were  for  primarily  10-mm-diam  tubes).  Their  correlation  is  given  as 

q"cHF  =  [G(450  +  Ai;)/ (40L/D  +  156G0’45)]  [1.02  -  (pr  -  0.54)2]  (33) 

where  q"cHF  is  m  W/cm2,  D  and  L  are  expressed  in  meters,  Ai;  in  kJ/kg,  G  in  kg/m2  sec, 
and  pr  in  bars.  The  expression  correlated  over  500  data  points  to  a  root-mean-square  error 
of  5.7  percent,  and  is  good  for  the  following  range  of  conditions:  pressure  from  120  to  200 
bar  (1,765  to  2,940  psia),  mass  velocity  from  G(p),  as  given  in  Fig.  43,  to  7,000  kg/m2  sec 
(up  to  1,435  lbm/ft2  sec),  inlet  subcooling  from  8°  to  272°C  (14°  to  490°F),  heated  length 
from  2,000  to  5,000  mm  (79  to  197  in.),  and  steam  quality  from  -0.3  to  0.6,  where  negative 
quality  represents  a  subcooled  liquid. 

Knoebel,  et  al.  (Ref.  227)  developed  a  CHF  correlation  for  subcooled  water  flow  based 
on  previous  work  performed  at  the  Savannah  River  Laboratory: 

Q'chf  =  K  (1  +  0.0515V)  (1  +  0.069ATsub)  (34) 

where  K  equals  153,600  for  stainless  steel  heaters.  The  correlation  is  limited  to  low-pressure 
water  flow  (less  than  45  psia)  and  higher  subcooling  (greater  than  45°F).  Green  and  Lawther 
(Ref.  228)  developed  a  CHF  correlation  for  high-pressure  water  as  well  as  Freon, 

q"CHF  -  Re?  (£?f/£>v)m  Pry  <  f  (Lsat/D)  (1  +  5)  (1  -  Sf)  (35) 


where 

n  =  1  —  exp  [-0.0067  (Lsat/D)] 
m  =  0.1  -I-  exp  (  —  0.007  (Lsat/D)] 
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p  -  -0.5  {0.15  +  exp  [-0.007  (Lsat/D)]} 
w  =  -{0.21  +  0.55  exp  [-0.007  (Lsat/D)]} 

f  (Lsat/D)  =  9  X  10“5  exp  {-0.00055  (Lsat/D)  +  3.83  exp  [-0.00396  (Lsat/D)]} 

8  =  exp  {-[0.14  x  108  (JN  +  0.02  (L/D)  Prv]} 

<tn  =  ff/(DAivQv) 

8(  =  0.75  exp  {-B(L/D)  (Vma/pf/if  Aiv)} 

B  =  130.5  exp  {5.0  exp  (-0.02  L/D)} 

The  expression  was  quoted  to  have  correlated  over  4,250  data  points  with  a  root-mean- 
square  error  of  less  than  10  percent,  although  the  correlation  proposed  by  Thompson  and 
Macbeth  (Ref.  218)  predicted  the  higher  pressure  data  better.  The  correlation  is  limited  to 
reduced  pressure  less  than  0.7,  mass  fluxes  greater  than  200  kg/m2  sec  (41  lbm/ft2  sec),  and 
exit  quality  greater  than  0.1. 

Shah  presented  graphical  correlations  for  CHF  in  subcooled  and  saturated  flows  in  tubes 
(Ref.  229)  and  annuli  (Ref.  230).  The  tube  correlation  in  functional  form  is 

Bo  =  ^(x^L/D)  for  Y  <  104 

Bo  =  f2(xcr,Y,pr)  for  Y  >  105  (36) 


where 

Y  -  (GDCp/kf)  (G2/ pfgD)0-4  W/O0'6 

The  graphical  expressions  correlated  90  percent  of  1,271  tube  data  points  within  ±30 
percent.  The  data  include  various  fluids  including  water  at  reduced  pressure  from  0.0012 
to  0.94,  mass  flux  from  6  to  24,300  kg/m2  sec  (1.2  to  4,980  lbm/ft2  sec),  and  inlet  quality 
from  -3.0  to  positive  values.  The  annulus  correlation  in  functional  form  is 


Bo  —  f3(xin,L/Dhp) 


(37) 
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where 


DhP  -  4  x  Flow  area/Heated  perimeter 

The  graphical  expression  correlated  88  percent  of  825  annulus  data  points  and  one 
noncircular  geometry  data  point  within  ±  30  percent.  The  data  include  various  fluids  including 
water  at  reduced  pressure  from  0.017  to  0.9,  mass  flux  from  100  to  15,780  kg/m2  sec  (20 
to  3,230  lbm/ft2  sec),  heated  perimeter  diameter  from  5.3  to  96.3  mm  (0.2  to  3.8  in.),  gap 
width  from  0.5  to  11.1  mm  (0.02  to  0.44  in.),  L/Dhp  from  3.7  to  335,  and  inlet  quality  from 
-3.1  to  0.0. 

Finally,  Vandervort,  et  al.  (Ref.  80),  recently  presented  an  optimized  CHF  correlation 
for  subcooled  forced  convection  in  tubes  based  on  the  five  parameters  of  mass  flux,  subcooling, 
pressure,  diameter,  and  length-to-diameter  ratio.  The  correlation  accounts  for  cross-coupling 
of  the  parameters,  and  is  given  as 

Q " CHF  =  17.05  [(G')°-0732  +  0.2390(D')J 

X  [(T ' )°.3°60(g')  +  0.001730(T')  -  0.0353  (D')] 

X  [(P') -0.1289] 

x  [1  +  0.01213(D')-2-946  +  °-7821(G')  +  0.009299(T')] 

X  [1.540  -  1.280  (L/D)'] 


where 


G'  =  0.005  +  G/105 
T'  =  5  +  ATsub 

P'  -  (0.0333  +  p)/3.0 

D'  =  D/0.003 

(L/D)'  =  (L/D)/40.0 

and  q"cHF  is  in  W/m2,  G  is  in  kg/m2-sec,  ATsub  in  degrees  Centigrade,  p  is  in  MPa,  and 
D  is  in  meters.  The  expression  correlated  85  percent  of  over  800  data  points  within  ±  25  percent. 
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Q"cHF,max  =  (GDifg/4L)  [1  +  cf  (ATsub)i/ifg]  (41) 

Gambill  and  Lienhard  (Ref.  231)  presented  an  interesting  semiempirical  approach  to 
determining  a  practical  upper  limit  to  CHF.  From  kinetic  theory  they  derived  an  expression 
for  Qniax,max>  the  highest  heat  flux  that  can  conceivably  be  achieved  in  a  phase-transition 
process: 


max, max  =  Ggifg  (RT/^Tr)"2  ^ 

A  practical  upper  limit  to  heat  transfer  by  phase  change  was  shown  to  be  0.1qmaXjmax; 
however,  the  method  is  limited  to  pressures  less  than  about  one-tenth  of  the  critical  pressure. 
As  a  point  of  interest,  the  highest  CHF  achieved  for  any  configuration  and  test  condition 
was  nearly  3.4  x  108  W/m2  (30,000  Btu/ft2  sec)  measured  by  Ornatskii  and  Vinyarskii  (Ref. 
232)  with  a  nonuniformly-heated  (circumferential),  small-bore  tube. 

4.6  COMPARISON  OF  FLOW  BOILING  AND  CHF  CORRELATIONS 

Guglielmini,  et  al.  (Ref.  172)  presented  an  interesting  comparison  of  flow  boiling 
correlations  including  boiling  incipience  and  partial  nucleate  boiling.  Most  of  the  boiling 
incipience  correlations,  when  compared  to  representative  data,  had  an  accuracy  within  30 
percent.  A  quantitative  comparison  of  the  partial  nucleate  boiling  expressions  could  not  be 
performed  because  of  lack  of  data.  The  correlations  for  fully  developed  nucleate  boiling 
generally  predicted  representative  data  to  within  ±  30  percent,  although  in  some  cases  the 
data  scatter  caused  error  in  the  prediction  to  be  as  large  as  70  percent. 

Gambill  (Ref.  233)  compared  ten  CHF  correlations  for  water  flow  in  a  tube  (diameter 
of  0.1  in.  and  length  of  20  in.)  at  a  pressure  of  600  psia  and  subcooling  of  100°F  over  a 
range  of  velocity  from  10  to  80  ft/sec.  Prediction  of  CHF  at  the  higher  velocity  using  the 
correlations  varied  by  nearly  a  factor  of  four  and  by  a  factor  of  two  at  the  lower  velocity. 
Factors  of  two  or  greater  between  correlation  predictions  at  higher  mass  velocity  have  also 
been  noted  by  Zeigarnik,  et  al.  (Ref.  96)  and  Boyd  (Ref.  17),  and  the  disagreement  tends 
to  worsen  as  mass  velocity  increases.  It  should  be  reiterated  that  empirical  CHF  correlations 
typically  have  been  derived  for  a  specific  range  of  experimental  conditions,  and  extrapolation 
of  a  correlation  outside  the  specific  range  can  lead  to  very  large  errors.  In  addition,  a  given 
analytical  CHF  correlation  typically  has  enough  “adjustable”  constants  in  the  predictive 
equations  to  permit  an  acceptable  correlation  of  a  specific  data  set;  however,  general 
application  of  the  correlation  to  various  configurations  or  test  conditions  may  result  in 
significant  errors. 
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5.0  APPLICABLE  EXPERIMENTAL  APPROACHES  AND  MEASUREMENT 
TECHNIQUES  FOR  BOILING  HEAT  TRANSFER 

A  representative  list  of  flow  boiling  experiments  that  have  been  performed  within  the 
last  50  years  is  presented  in  Appendix  B,  and  a  few  comments  are  in  order.  The  heating  method 
of  choice  is  electrical  resistance  heating  (Joule  heating).  A  fewer  number  of  experiments  made 
use  of  cartridge  heaters,  and  very  few  incorporated  some  other  means  of  heating.  The  test 
sections  were  primarily  tubes,  although  a  significant  number  of  experiments  had  annular 
test  sections.  Of  the  more  than  175  listed  experiments  where  water  was  used,  46  had  a  maximum 
heat  flux  greater  than  1  x  107  W/m2  (1,000  Btu/ft2  sec),  and  only  10  had  a  heat  flux  greater 
than  5.5  x  107  W/m2  (5,000  Btu/ft2  sec),  with  most  of  these  performed  at  a  pressure  less 
than  35  bar  (500  psia),  subcooling  less  than  170°C  (300°F),  or  mass  velocity  less  than  30 
Mg/m2  sec  (6,200  lbm/ft2  sec).  Clearly,  data  obtained  at  high  heat  flux,  pressure,  subcooling, 
and  velocity  are  severely  limited. 

A  number  of  measurement  techniques  used  in  these  and  other  experiments  are  worth 
reviewing.  One  of  the  most  important  measurements  in  a  boiling  experiment  is  surface 
temperature.  Placement  of  a  discrete,  intrusive  temperature  measuring  device  such  as  a 
thermocouple  on  the  surface  at  the  fluid  interface  may  disrupt  the  flow  field,  causing  large 
errors  in  the  temperature  measurement.  For  this  reason,  backside  wall  temperature  is  typically 
measured  and  used  to  estimate  surface/fluid  interface  temperature  through  the  use  of  an 
analytical  temperature  distribution.  On  an  electrically  heated  surface,  the  temperature 
measurement  with  a  thermocouple  (which  provides  temperature  from  millivolt  changes)  is 
further  complicated  by  voltage  drop  in  the  test  specimen  due  to  current  flow.  Hughes  (Ref. 
21)  incorporated  a  single  wire  junction  at  the  surface  and  the  two-wire  thermocouple  junction 
formed  0.8  mm  (1/32  in.)  from  the  surface,  thereby  preventing  any  extraneous  voltages  due 
to  voltage  drop  between  the  thermocouple  leads.  Numerous  researchers  have  used  thin  sheets 
of  mica  or  other  electrical  insulating  material  to  isolate  the  thermocouple  from  the  current- 
carrying  surface.  One  of  the  most  attractive  approaches  as  suggested  by  Dutton  and  Lee  (Ref. 
234)  involves  the  use  of  a  three-wire  thermocouple.  The  thermocouple  circuit  is  “balanced” 
by  applying  and  subsequently  reversing  the  current  and  nullifying  the  voltage  drop  with  a 
voltage-dividing  potentiometer,  the  net  effect  being  the  same  as  welding  all  three  wires  to 
the  location  of  the  middle  wire. 

Burnout  detectors  have  been  used  by  a  number  of  investigators  (e.g.,  Refs.  21, 235,  and 
236)  to  protect  the  test  sections  from  destruction  as  the  CHF  is  approached.  However,  as 
pointed  out  previously,  Hughes  (Ref.  21)  noted  that  the  speed  and  intensity  of  the  transition 
to  film  boiling  (i.e.,  burnout)  was  observed  to  increase  with  increasing  velocity  and  subcooling, 
and  in  some  cases  burnout  has  occurred  instantaneously  (Ref.  147).  Therefore,  the  effectiveness 
of  a  burnout  detector  to  protect  the  test  section  in  highly  subcooled,  high-velocity,  high- 
pressure  flow  is  questionable. 
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Gunther  (Ref.  27)  was  one  of  the  earliest  to  perform  a  detailed  photographic  study  of 
subcooled  flow  boiling  at  higher  heat  flux,  although  pressure  was  limited  to  below  11.3  bar 
(164  psia)  and  velocity  limited  to  below  12  m/sec  (40  ft/sec).  His  photographic  system  was 
capable  of  20,000  frames/sec  and  gave  a  resolution  of  bubbles  as  small  as  0.1  mm  (0.004 
in.)  diameter.  Fiori  and  Bergles  (Ref.  20)  used  several  photographic  techniques,  including 
high-speed  photographs,  movies,  and  video.  They  successfully  identified  cyclic  variations 
in  surface  temperature  from  variations  in  the  flow  structure  by  simultaneously  recording  surface 
temperature  and  high-speed  photographs  of  a  specific  surface  location.  Ram,  et  al.  (Ref. 
237),  although  studying  bubble  generation  in  pool  boiling  of  water,  analyzed  intensity 
fluctuations  registered  by  a  photomultiplier  tube  to  determine  bubble  size  and  growth  cycle. 
Lineberger  (Ref.  238)  used  a  dual-frequency  sound  field  to  determine  bubble  size  in  pool 
boiling  of  water.  Brown  (Ref.  40)  used  reflected  light  from  bubbles  detected  by  a  light- 
dependent  resistor  to  identify  the  presence  of  first  vapor  when  boiling  begins. 

Several  approaches  have  been  used  to  determine  the  vapor  thickness  or  void  fraction  near 
the  heated  surface  in  local  boiling.  Costello  (Ref.  239)  measured  the  number  of  beta  particles 
emitted  from  a  vial  containing  radioactive  strontium  90  to  determine  the  vapor  thickness 
in  a  low-velocity,  water  flow  boiling  apparatus  with  an  annular  test  section.  Rogers,  et  al. 
(Ref.  240)  used  a  traversing  gamma-ray  densitometer  (cobalt-57  source)  to  measure  void 
fractions  in  water  flow  boiling  in  an  annulus.  Similar  densitometers  have  been  used  by 
Buchberg,  et  al.  (Ref.  100)  and  Edelman,  et  al.  (Ref.  241)  to  measure  the  void  fraction  of 
subcooled  water  flow  in  tubes.  Jiji  and  Clark  (Ref.  242)  used  a  specially  developed  traversing 
thermocouple  probe  to  measure  bubble  boundary-layer  thickness  and  temperature  profiles 
in  low-velocity,  subcooled  water  flow.  Later,  Stefanovic,  et  al.  (Ref.  243)  used  a  similar  probe 
to  measure  temperature  fluctuations  within  a  superheated  boundary  layer  adjacent  to  a  heated 
annular  wall  in  forced  water  flow. 

The  measurement  of  acoustic  noise  due  to  boiling  has  permitted  various  researchers  to 
study  hydrodynamic  instabilities  in  heated  channels  and  to  diagnose  the  different  regimes 
of  boiling.  Pressure  transducers  have  typically  been  incorporated  in  an  experimental  apparatus 
where  flow  oscillations  are  to  be  studied  (e.g.,  Refs.  244  and  245),  although  Romberg  and 
Harris  (Ref.  246)  used  piezoelectric  accelerometers  to  measure  vibrations  associated  with 
parallel  channel,  density  wave,  and  acoustic  oscillations  in  a  flow  boiling  loop.  Several 
investigations  have  been  performed  to  show  that  information  on  the  boiling  mode  can  be 
ascertained  from  boiling  noise  (Refs.  247-251).  Hydrophones  (i.e.,  transducers)  with 
piezoceramic  sensing  elements  have  typically  been  used  to  measure  the  acoustic  emissions 
from  boiling. 

A  number  of  important  considerations  can  be  learned  from  the  success  or  failure  of 
previous  experiments  and  associated  analyses.  Westwater  (Ref.  11)  pointed  out  that  one  basic 
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method  of  proving  the  correctness  of  instrumentation  in  these  types  of  experiments  which 
is  widely  overlooked  is  a  simple  energy  balance  of  the  system.  Accuracy  of  the  all-important 
wall  temperature  can  be  severely  affected  by  overlooked  items  such  as  conduction  effects 
or  temperature-dependent  material  properties.  Flow  loop  layout  can  affect  the  hydrodynamic 
stability  of  the  system,  which  can  seriously  degrade  the  heat-transfer  and  CHF  measurements. 
Gas  content  and  cooling  fluid  impurities  have  been  shown  to  have  an  effect  on  boiling  heat 
transfer  and  CHF,  yet  in  many  experiments  these  have  been  ignored  or  have  not  been 
adequately  quantified.  As  shown  in  Appendix  B,  many  researchers  fail  to  present  basic  test 
condition  values  such  that  their  data  may  be  used  in  other  analyses.  In  many  cases,  the  location 
in  which  the  subcooling  and  pressure  are  measured  or  determined  is  not  given.  Finally,  the 
geometric  scale  of  a  heater  surface  has  been  shown  to  be  very  important.  Specifically,  Bakhru 
and  Lienhard  (Ref.  252)  demonstrated  that  the  hydrodynamic  processes  that  give  rise  to 
burnout  cease  to  occur  when  the  Laplace  number*,  R' ,  is  reduced  below  a  value  on  the  order 
of  0.1  (see  also  Refs.  33  and  253).  Generalized  correlations  based  on  data  acquired  from 
a  configuration  with  a  Laplace  number  below  this  value  may  have  very  large  errors  when 
applied  to  scaled-up  or  different  configurations. 

6.0  SUMMARY 

An  extensive  review  of  backside  water  cooling  processes  that  are  applicable  to  high-enthalpy 
facility  components  has  been  conducted.  The  processes  can  be  identified  on  the  basic  boiling 
curve  which,  for  the  configurations  of  interest,  includes  the  pure  convection  regime,  the 
nucleate  boiling  regime,  and  the  boiling  crisis  point  where  the  CHF  occurs.  More  than  20 
parameters  have  been  shown  to  affect  various  portions  of  the  boiling  curve,  and  optimum 
cooling  (and  greater  CHF)  occurs  at  high  subcooling  and  high  mass  velocity.  The  primary 
parameters  which  affect  steady-state  boiling  heat  transfer  are  pressure,  mass  velocity, 
subcooling,  heater  diameter,  and  heater  length  (Ref.  7);  however,  other  parameters  have  been 
experimentally  shown  to  have  significant  effects.  Several  heat-transfer  enhancement  techniques 
have  also  been  reported. 

Because  of  the  complexity,  the  unsteady  nature,  and  the  small  scale  of  the  flow  boiling 
processes,  no  general  theory  has  been  developed.  Reasonable  confidence  in  correlating  heat- 
transfer  data  in  the  pure  convection  and  nucleate  boiling  regimes  is  shown.  The  more  recent 
pure  forced-convection  correlations  of  Petukhov  [Eq.  (5)]  for  turbulent  flow  in  a  tube  and 
Kays  and  Leung  (Ref.  153)  for  turbulent  flow  in  an  annulus  are  probably  the  most  useful 
for  the  nonboiling  regime.  Bergles  and  Rohsenow’s  partial  nucleate  boiling  expression  [Eq. 
(8)]  best  correlates  the  transition  from  pure  convection  to  fully  developed  nucleate  boiling, 


*  The  Laplace  number  is  defined  as  R'  =  r[g(gr  -  ev)/tf]1/2,  where  r  is  a  characteristic  radius  dimension  of  the  heater  (see 
Ref.  252). 
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and  the  widely  accepted  Rohsenow  relation  [Eqs.  (16)  and  (17)]  ranks  among  the  best 
correlations  for  fully  developed  nucleate  boiling  in  forced  flow.  Considerable  disagreement 
between  the  large  number  of  analytical  and  empirical  CHF  correlations  currently  exists. 
Moreover,  because  the  CHF  correlations  are  sensitive  to  the  range  of  conditions  for  the  data 
used,  no  single  correlation  can  be  recommended.  Those  CHF  correlations  developed  with 
data  from  high  mass  velocity,  pressure,  subcooling,  and  heat  flux  experiments  using  water 
as  the  coolant  are  more  appropriate  for  use  in  the  design  of  high-enthalpy  arc  heater 
components.  Recommended  correlations  in  this  group  include  Bernath  (Ref.  203),  Van  Huff 
and  Rousar  (Ref.  29),  Yagov  and  Puzin  (Ref.  195),  Levy  (Ref.  204),  Labuntsov  (Ref.  216), 
and  Vandervort,  et  al.  (Ref.  80).  Additional  correlation  development  will  be  required  to  account 
for  factors  such  as  acceleration  or  transient  effects  at  these  extreme  conditions. 

Very  few  data  exist  that  were  obtained  at  high  heat  flux,  high  subcooling,  high  mass 
velocity,  and  high  pressure,  which  are  optimum  for  cooling  high-enthalpy  facility  components. 
Finally,  numerous  shortcomings  of  previous  experiments  have  been  identified,  and  care  must 
be  taken  in  any  experimental  program  where  boiling  heat-transfer  data  are  acquired  for  future 
analysis  and  correlation. 
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Figure  2.  Effect  of  pressure  and  velocity  on  subcooled  boiling. 
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Figure  5.  Effect  of  liquid  velocity  on  bubble  growth  and  collapse. 
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Figure  7.  Subcritical  pressure  water  burnout  heat  flux. 
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Figure  13.  Effect  of  surface  roughness  on  heat  transfer. 
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Figure  15.  Burnout  heat  flux  of  knurled  surfaces. 
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Figure  16.  Effect  of  a  surface  active  agent  on  heat  transfer. 


80 


10 


'5 


AEDC-TR-93-3 


P  =  6.6  bars 
AT  Sat  =  50°C 

AC  =  Excess  Concentration  of 

the  Volatile  Component 
in  the  Vapor 


Concentration  of  Ethanol,  percent 


a.  Ethanol-water 


from  Ref.  63 


Concentration  of  Acetone,  percent 
b.  Acetone-water 

Figure  17.  Effect  of  binary  mixture  concentrations  on  critical  heat  flux  (CHF). 
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Figure  19.  Effect  of  wall  material  on  CHF. 
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Figure  20.  Effect  of  wall  thickness  on  heat  transfer. 
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Figure  21.  Influence  of  tube  diameter  on  CHF  at  various  subcoolings. 
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Figure  22.  Effect  of  tube  diameter  on  CHF. 
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Figure  24.  Effect  of  tube  length  on  CHF  at  various  subcoolings. 
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Figure  25.  Effect  of  position  of  the  heating  surface  on  CHF. 
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Figure  27.  Effect  of  curved  test  section  on  CHF. 
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Figure  28.  Effect  of  undulating  flow  on  heat  transfer. 
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Figure  29.  Pressure  effect  with  flow  swirl. 
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c.  G  =  1,700  kg/m2-sec 
Figure  30.  Mass  velocity  effect  with  flow  swirl. 
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ATsat,  K 


a.  D  =  1.2  mm  and  Tsub  =  10  K,  inlet  b.  D  =  0.8  mm  and  ATsub  =  30  K,  inlet 

Figure  33.  Effect  of  transient  heating  on  boiling  curve. 
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from  Ref.  101 


Figure  34.  Effect  of  flow  restriction  on  flow  stability  and  burnout. 
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Figure  35.  Parametric  trends  for  subcooled,  forced-convection  cooling. 
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Figure  36.  Effect  of  heat-transfer  rate  on  bubble  characteristics. 
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Stainless  Steel  0.094  in.  I.D. 
p  =  22  lbf/in.2  abs 
v  V  =  1.7  ft/sec 


Figure  37.  Construction  of  the  partial  boiling  curve. 
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O  Addoms,  Water,  at  14.7,  283,  770, 1,205, 1 ,985,  2,465  psia  (Ref.  1 5) 
+  Cichel ll-Bon i lla.  Benzene,  at  14.7,  55, 115,  265,  51 5  psia  (Ref.  36) 
v  Cichell i-Bonilla,  n-Heptane,  at  6.6, 14.7,  50,  1 15,  21 5  psia  (Ref.  36) 
9  Piret-lsbin,  Water,  Carbon  Tetrachloride,  Iso-propyl  and 
n-Butyl  Alcohol,  at  14.7  psia  (Ref.  260) 
x  Cichel  I  i-Bon  ilia,  n-Pentane,  at  22,60, 1 1 5,  21 5,  31 5  psia  (Ref.  36) 

A  Cichelli-Bonilla,  Ethanol,  at  14.7,  55,  115,  265,  515  psia  (Ref.  36) 
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Figure  38.  Determination  of  the  coefficient  BL  in  the  Levy  correlation. 
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Figure  39.  Dependence  of  overall  pressure  drop  on  operating  conditions. 
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Upper  Limit  of  G  x  10-6  Ibm/hrft2 

Figure  41.  Approximate  limits  of  velocity  regimes  for  round  tubes. 
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b.  C'  for  p  =  560  psia 
Figure  42.  Continued. 
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c„  A'  for  p  =  1,000  psia 
Figure  42.  Continued. 
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d.  C'  for  p  =  1,000  psia 
Figure  42.  Continued. 
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G  x  10-6  lbm/hr  ft2 


e.  A'  for  p  =  1,550  psia 
Figure  42.  Continued. 


012  3456789 

G  x  10-6  Ibm/hrft2 

f.  C'  for  p  =  1,550  psia 
Figure  42.  Continued. 
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Figure  42.  Continued. 
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from  Ref.  218 


h.  C '  for  p  =  2,000  psia 
Figure  42.  Concluded. 
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Figure  43.  Mass  velocity  range  for  the  Becker  burnout  correlation. 
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APPENDIX  A 
Table  of  Conversion  Factors 


Multiply 

By 

To  Obtain 

Parameter 

ft2 

9.290E-02 

m2 

A 

ft/sec 

3.048E-01 

m/sec 

V 

ft 

3.048E  — 01 

m 

D,L 

Btu/lbm°F 

4.187E  +  00 

kJ/kg°C 

cp 

Btu/ft2hr 

3.155E  +  00 

W/m2 

q",CHF,q, 

in. 

2.540E  +  01 

mm 

d 

lbm/ft2hr 

1.356E-03 

kg/m2sec 

G 

ft/hr2 

2.352E  — 08 

m/sec2 

g,a 

Btu/ft2hr°F 

5.679E  +  00 

W/m2°C 

h 

Btu/lbm 

2.326E  +  00 

kJ/kg 

Btu/ft  hr°F 

1.731E  +  00 

W/m°C 

k 

psia 

6.895E-02 

bar 

P 

psia 

6.895E  +  03 

Pa 

P 

°F  (differential) 

5.556E  — 01 

°C,  K 

ATsub>  AT, 

ft2/hr 

2.581E  — 05 

t 

m2/sec 

a 

Ibm/ft  hr 

4.134E  — 04 

kg/m  sec 

lbm/ft3 

1.602E  +  01 

kg/m3 

<? 

lbf/ft 

1.459E  +  01 

N/m 

a 

lbf/ft2 

4.788E  +  01 

N/m2 

T 

°C  =  5/9  (°F  —  32) 

K  =  °C  +  273.15 
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Forced-Convection  Boiling  Experiments 
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!  55.2-138 

.16-2 

.27-2.2 

.48-. 56 

Bankoff  k  Mason 

274 

62  1  Northwestern 

X 

X 

brass 

water 

X 

X 

horiz 

X 

X 

i  X 

1  1 

1 

X 

x  i  18-73 

L 

• 

!  1 

1  207 

? 

1 

bubble  growth;  jet  iepingaent 

Swenson,  et  al. 

275 

62  1  Babcock  k  Mil 

X 

X 

SS 

water 

X 

X 

vert 

ac 

:  x 

X  1 

1 

1 

x  1  0 

L 

1  176 

.28-. 57 

1  .95,  1.36 

straight  k  swirl  flow 

Swenson,  et  al. 

120 

62  ! Babcock  k  Mil 

X 

X 

SS 

water 

X 

X 

vert 

ac 

1  X 

X  1 

1 

1 

! 

X 

i  0 

L 

1  164 

i  130 

.22-1.0 

!  .68-1.76 

Styrikovich,  et  al. 

118 

62  1  USSR 

X 

X 

X 

SS 

water 

1 

vert 

ac 

1 

X  1  X 

X  *  X 

I 

! 

• 

X 

:  .4-2 

276 

62  !  MSA 

Inconel 

water 

X 

vert 

dc 

1 

X  1  X 

X  1  X 

! 

x  : 

X 

!  x  !  x 

1 

I  105-191 

A 

1  28-61.5 

!  4.8-21.2 

0-133 

22-62 

!  22-62 

Papell 

277 

63  !  NASA 

X 

X 

Inconel  X 

water 

X 

horn? 

ac 

X 

X 

X 

! 

! 

!  : 

! 

!  :  100-146 

L 

1  21 

1  2.6-12.4 

.6-2.6 

1.2-3.B 

Judd  k  Hade 

278 

63  !  McMaster  U 

0 

0 

X 

SS 

water 

X 

vert 

X 

1 

X 

1  x  1 

1 

x  !  ? 

!  16-96 

1  ? 

? 

eccentricity  effects 

Becker 

279 

63  !  Sweden 

X 

X 

SS 

water 

X 

vert 

dc 

X 

X 

X 

: 

x  !  0 

L 

:  64 

1  2.5-10.0 

.5-1.2 

.032-.  33 

tee  t  Obertelli 

03,280 

63  {  UK  AEEN 

x 

X 

SS 

water 

X 

vert 

dc 

X 

X 

X 

X 

X 

1  X 

X 

X 

I  2.8-139 

I 

!  20-360 

!  38.6-110 

0-7.7 

!  ,4-4 

Matzner 

04 

63  !  Columbia  li 

X 

X 

? 

water 

X 

vert 

dc 

X 

X 

X 

X 

X 

X 

!  9-210 

I 

I  26-151 

TW 

T^Ti 

!  .68-18.6 

Styrikovich,  et  al. 

281 

63  1  USSR 

X 

X 

X 

SS 

water 

X 

vert 

ac 

X 

X 

X 

X 

X 

i  0-42 

L 

1  20-118 

:  101 

0-9.9 

I  .4-2 

Janssen,  et  al. 

215 

63  i  6£ 

X 

X 

X 

? 

water 

X 

vert 

X? 

X 

X 

X 

X 

X 

X 

X 

X 

x  I  0-33 

E 

!  33-600 

!  41.4-100 

0-5.5 

1  .2-0.4 

eccentricity 

Bergles 

05 

63  !  MIT 

X 

X 

SS 

water 

X 

X 

X 

horiz 

dc 

X 

X 

X 

X 

X 

X 

X 

!  0-50 

E 

!  5-35 

i  1. 4-5.9 

!  1.5-24.3 

Tippets 

28 

64  :  GE 

X 

X 

X 

SS 

water 

X 

vert 

X 

X 

X 

X 

X 

X 

1  8-144 

1 

!  17.6 

!  60.7-72.0 

1  .24-1.9 

Jiji  k  Clark 

242 

64  iNYU,  U  of  MI  !  x 

1 

X 

SS 

water 

X 

X 

vert 

dc 

X 

X  I  X 

X 

x  !  28-167 

I 

I  20 

I  13.8-69 

.44-4.7 

.3-1.8 

flow  teap  probe  aeasureaent 

Becker  k  Hernborg 

89 

64  1  Sweden  !  x  ! 

X 

SS 

water 

; 

X 

vert 

dc 

X 

X 

X 

X 

1 

X 

X 

!  60-205 

I 

1  81 

f  8.3-36.8 

0-3 

1  ,07-. 96 

Alekseev,  et  al. 

92 

64  1  USSR 

X 

X 

X 

X 

X 

•5 

water 

1 

X 

various 

X 

X 

X 

x  1  X 

X 

x  i  ? 

I  40,67 

!  29.3-207 

0-11 

i  0-5.6 

aultiple  flow  configurations 

145 

64  1  MIT 

X 

X 

SS 

water 

X 

X 

vert 

i  X 

X  i  X 

X 

:  103 

I 

!  24 

i  1.5 

0-12.6 

2. 1-8,5 

Bergles  k  Rohsenow 

161 

64  !  HIT 

X 

X 

X 

X 

SS 

water 

X 

X 

X 

1 

honz 

dc 

1 

X  1 

X  1 

X 

I  0-100 

I 

!  20-58 

1  17.9 

0-1B.V 

0-17.4 

1 

Jeglic  k  Yang 

245 

65  !  NASA 

X 

X 

X 

Inconel  X 

water 

X 

X 

vert 

ac 

x  !  x  1  x 

X 

X 

X 

I  0-100 

T 

!  100 

!  .2-4 .1 

0-2.65 

.7-2.5 

Hosier 

282 

65  1  Hestinghouse 

X 

X 

x 

SS 

water 

X 

horiz 

dc 

! 

x  !  x  1  x 

! 

1 

i 

I 

i  X 

!  0-28 

I 

1  102 

1  10,4-41.4 

0-3.15 

.14-1.4  [primarily  photo  data 

Durant,  et  al. 

47 

65  idu  Pont/SRL 

X 

X 

SS 

water 

X 

vert 

dc 

X 

X 

X 

X 

!  x 

X  1 

i 

1 

!  38-88 

? 

!  40-96 

1  2. 1-4.1 

7-11.3 

3-8.5 

[knurled  surfaces 

Feldaanis 

287- 

65  1  AFFDL 

X 

X 

copper 

water 

X 

honz 

X 

X 

1 

1 

* 

i 

! 

1  0 

1 

!  118 

1  .03-. 2 

? 

.03-22 

•001-.0014 

ilow  gravity  flow 

Stenning  k  Veziroglu 

244 

65  !  U  of  Miami 

X 

X 

mchroee 

Freon  11 

honz'’ 

dc 

X 

X 

:  x 

I 

!  32 

1 

!  254 

!  5.2 

? 

.003-. 010 

Ornatskii  k  Vinyarskii 

232 

65  1  USSR 

X 

X 

Cu  alloy, brass/Ag 

water 

X 

vert? 

dc 

X 

X 

X 

X 

X 

X 

X 

!  50-190 

1  7-35 

1  11.2-61.9 

0-343 

I 

.01— .09 

highest  heat  flux  ever  achieved 

Waters,  et  al. 

284 

65  !  GE 

X 

X 

inconel 

water 

X 

honz,  vert 

dc 

x 

X 

X 

X 

X 

!  0-270 

I 

i  327 

!  70-105 

2-5.4 

6.45-9.6 

upstreaa  burnout 

Rousar  k  Van  Huff 

285 

66  !  Aerojet/JPL 

X 

X 

Inconel  718, SS 

H202, water 

X 

vert 

dc 

X 

X 

X 

X 

X 

X 

X 

X 

X 

x  !  151-331 

A 

!  23,31 

[  69-207 

0-45.7 

11.6-45.7 

286 

66  !  0t:  State  U 

X 

X 

X 

Inconel  600 

water 

X 

X 

vert 

dc 

X 

X 

X 

x  !  0-51 

1 

!  245 

!  1 

.06-. 25 

.009-.  039 

coiled  flow 

Todreas  !  Robsenow 

287 

66  1  AEC/MIT 

X 

X 

X 

aluiunua 

water 

X 

not  defined 

dc 

!  x 

X 

X 

X 

X 

X  1 

:  0-189 

I 

i  140-224 

1  2.6-15.9 

? 

.68-2.7 

253 

66  !Can/Colu».  U 

X 

x 

X 

water 

X 

vert 

dc 

X 

X 

X 

X 

1 

1 

!  356-579 

1  69 

0-2,7 

.69-6.8 

stripping  rings,eccentricity 

Haultoetsch  k  Griffith 

289 

66  1  HIT 

X 

x 

55  <§> 
Hastell  oy-C^ 

water 

X 

X 

horiz? 

dc 

X 

X 

X 

X 

X 

X 

X 

X 

!  108-139 

1 

i  25-200 

1  2. 1-5.8 

0-22 

3-6.8 

Larsen,  et  al. 

290 

66  !  Hestinghouse 

X 

X 

x,  sou  are 

water 

X 

vert 

dc 

X 

X 

X 

X 

[  7-70 

1 

1  118 

[  105-160 

0-4.9 

.66-4,1 

291 

66  1  MIT  lx 

X 

X 

1 

SS 

water 

X 

X 

x 

vert 

dc 

X 

X 

X 

X 

X 

!  0-139 

I 

1  60-240 

1  34.5-69 

.54-5.4 

40 

67  1  MIT 

X 

SS 

water 

X 

X 

horiz 

dc 

X 

X 

X 

X 

X 

!  6-100 

7 

!  20 

i  2. 1-6.2 

0-1.6 

0-1.7 

bubble  detector 

Frost  k  Kippenhan 

54 

67  1  UT/UNA 

X 

X 

SS 

water 

X? 

vert  1 

dc 

X 

X 

X 

X 

x  !  44-03 

i 

1  23 

1  2.1 

0-6.3 

.24-1.2 

varied  surface  tension 

Glickstein  k  Whitesides 

292 

67  iPrattkttutney 

X 

X 

X 

Inconel 

var 

vert  1 

dc 

X 

X 

X 

>  0-211 

L 

I  98.5 

!  10.4-55.2  !  .27-2.1 

.2-4.5 

nucleate  k  file  boiling 

Heisaan,  et  al. 

293 

68  [Hestinghs/PM. 

X 

1  X 

Inconel 

water 

x 

:  x 

1 

!  vert 

! 

x?  ! 

X 

X 

X 

X 

i  22-128 

I 

!  20-59 

!  125-160 

0-2 

.68-4 

Fieri  k  Bergles 

20 

68  1  HIT 

X 

X 

X 

SS 

watEr 

X 

X 

X 

vert 

1 

dc 

X 

X 

X 

X 

X 

X 

X 

X 

!  20-61 • 

E 

!  15-54 

!  1. 7-6.2 

3.15-17.3 

.68-10.2 

Markels  k  EXirfee 

142 

&8  !  ARC,  VA 

X 

X 

SS 

water 

X 

vert 

1 

dc 

X 

X 

X 

X 

X 

X 

x  i  0-258 

I 

!  31-80 

1  1.4-104 

0-5.4 

.2-2.7 

applied  a-t  field 

Lopina  k  Bergles 

2?4 

69  1  KIT.USAF  Ac. 

X 

X 

nickel 

water 

X 

x 

x 

horiz? 

dc 

X 

X 

X 

X 

X 

X 

x  :  ? 

i  0-58 

1  69-76 

"mrs?! 

? 

0-2.5 

.(llj-.lB 

nonboiling,  swirl  flow 

Borishanskjy  k  Fokin 

64 

69  1  USSR 

X 

X 

X 

X 

X 

SS, mchroee 

var, water 

X 

horiz, vert 

X? 

X  1  X 

X 

X 

X  1 

X 

? 

!  0.4-33 

!  1-24.6 

0-4.9 

1 

unsteady  heat  application 

Cof Field,  et  al. 

209 

69  !U  Pittsburgh 

X 

X 

SS 

Freon-113 

vert 

dc 

X  1  X 

X 

X 

1 

1 

1  2.8-24 

I 

I  74.6 

1  9.1-21.9 

0-.57 

I  1.3-5. 6 

DN3  detector 

Roeaer 

127 

69  1  Stanford  U 

X 

X 

SS 

water 

X 

X 

vert 

dc 

1 

X 

1 

1 

1 

i 

x  !  X 

x  ! 1.8-2.4,40  ? 

1  5.3-16 

0-2.3 

1 

1 

power  transients,  aging  studies 

Miropol'skiy  k  Pikus 

295 

69  !  USSR 

X 

X 

X 

SS 

water 

X 

not  defined 

dc 

X 

X 

X 

X  1 

1 

!  x  !  0-106 

? 

!  ? 

!  90-215 

0-4 

.1-2 

curved  channels 

Tolubinskiy,  et  al. 

296 

69  i  USSR 

X 

X 

X 

SS 

water 

X 

X 

vert 

X? 

X 

X 

X 

X 

X 

X 

1  :  0-143 

1 

1  7.7-750 

!  80-220 

0-7 

.05-.  8 

outer  wall  heated 

Bergles  k  Doraer 

173 

69  1  MIT 

X 

X 

SS .nickel 

water 

X 

X 

X 

honz 

dc 

X 

X 

X 

X 

X 

X 

X 

X 

!  !  58-141 

I 

!  24-195 

1  2. 1-5.5 

0-17.3 

1.5-18.3 

Peskov,  et  al. 

297 

69  1  USSR 

X 

X 

X 

SS 

water 

X 

vert 

ac 

X 

X 

X 

X 

X 

!  !  161-309 

I 

!  31-375 

:  101-203 

.58-4.9 

.82-5 

Zenkevich,  et  al. 

90 

70  •  USSR 

X 

X 

X 

X 

SS 

water 

X 

not  defined 

ac 

X 

X 

X 

X 

X 

X 

X 

!  !  0-60 

I 

1  93-730 

!  39,3-196 

0-4.5 

.3-2.8 

143 


APPENDIX  B 
Continued 


AEDC-TR-93-3 


: 

TYPE 

SURFACE 

COOLANT 

HEATING 

EFFECTS  STUDIED 

PARAMETER  RANGE 

Investigator 

Ref 

Yr  1  fr-ganixation 

Subcool 

Sat 

Tube 

Annulus 

Rod  / 

Tube 

Plate, 

Ribbon 

Katl 

Type 

Treatment  Ill 

Flow 

Coal, Gas 

Elect 

Elect  !  Liquid 

Steam, 

dTsub 

Press 

Vel 

Burnout 

Gas 

Accel 

L/D 

Hall  !  Dh 

Hall  i  RtXJOh- 

Press 

Uni /Nome* 

Stability 

Transient 

Optical 

Other 

dTsub 

1  L/D 

1  Press 

i  qdot  (6) 

Vel 

mdot  17) 

6  18) 

Comments 

1 

Isurfl 

(bulk) 

Hire 

Bundle 

Channel 

di 

dei 

dem 

deg 

gas 

untr 

nwtd 

Orientation 

Fired 

Htr 

Resis  I  Metal 

Fluid 

12) 

13) 

crit  pt 

Contnt 

£4) 

Matl  115) 

Thk 

ness 

Drop 

Ht  Flux 

deg-C  12) 

;  14) 

1  bar 

!  MH/sq-m 

m/s 

kg/s 

Mg/m2  s 

Tolubinskiy  l  Kostanchuk 

38 

70  :  USSR 

X 

X 

SS 

water 

X 

x 

vert 

dc  ! 

„ 

x 

x 

X 

X 

5-60 

? 

i  ? 

I  HO 

!  .05-1 

.08-.  2 

bubble  growth  measurements 

TolubinsViy ,  et  al. 

139 

70  J  USSR 

X 

X 

X 

SS 

water 

X 

vert 

x?  I 

X 

X 

X 

X 

X 

1  X 

0-313 

? 

!  67-1500?  1  50-200 

l  0-6.5 

.25-3 

varied  heating  modes 

Sterman  6  Nekrasow 

298 

70  ;  USSR 

X 

X 

steel 

water 

X 

vert 

ac  i 

X 

X 

X 

X 

1 

1  0 

? 

I  ? 

1  29.4-98 

1  .5-10.5 

.83-3.4 

5tefanovic,  et  al. 

213 

70  !  Yugoslavia 

X 

X 

X 

7 

water 

X 

vert 

dc  !  : 

X 

1  X 

x  !  0? 

? 

!  27,69 

:  l 

i  .07-. 27 

.26-. 34 

b.l.  probe  measurements 

Qsaachkin  i  Borisov 

299 

70  1  USSR 

X 

X 

X 

SS 

water 

X 

vert 

dc  1  i 

X 

x  1  ? 

i  519 

i  70 

:  .71 

6.9 

Saddis  i  Hall 

300 

70  1U  of  Mane, Eng 

X 

X 

copper,  Ni  plated 

water 

X 

horiz? 

t 

i  i 

X 

X 

i 

t  i 

X 

1 

!  x 

I  20-60 

L 

1  63 

1  ? 

:  o-.i5 

•  1-.5 

: 

drilled  cavity 

Johnson 

301 

70  !U  of  CA.Berkl 

X 

X 

nickel/iron 

water 

! 

X 

horiz 

x  1 

X 

X 

X  t  X 

x  :  i 

! 

• 

X  «  X 

!  5.6-62 

7 

!  7,4 

!  1-138 

:  o-.os 

.3,  4.3 

1 

Staub  l  Kaleet 

175 

70  i  re 

X 

X 

X 

7 

R-22, water 

! 

1  X 

vert 

?  i 

X 

X 

X  1 

:  ixi 

1  X 

1 

4-9! 

L 

1  25.8,150 

!  .1-69 

:  o-.oe? 

i  .18-2.8 

void  fraction  measurements 

Kawaaura,  et  al. 

302 

70  i  Japan 

X 

X 

X 

X 

SS,  nickel 

water 

1  X 

vert  ! 

x?  ! 

X 

X  1  X 

1 

:  i  i  x 

1 

1 

X  1  X 

0-70 

? 

!  2,2 

1  1 

i  0-15 

0-2 

1 

transient  power  application 

Ornatskii,  et  al. 

208 

71  1  USSR 

X 

X 

X 

!  1 

1  SS 

water 

i  X 

X 

1 

not  defined  1 

x?  !  I 

X 

X 

X  1  X 

! 

X  !  !  x  ! 

[ 

! 

0-07 

? 

!  30-150 

FKffSnK 

:  0-17,0 

1.5-8 

Veziroglu  k  Lee 

137 

71  IU  of  Miami 

X? 

X 

!  ni  throw 

Freon-11 

I 

vert  1 

dc  !  ! 

X 

1  X 

X 

? 

!  180,254 

5  ? 

'.WW-.OiS 

parallel  channels 

Motley,  et  al. 

67 

72  !Westingh/CMU 

x  1  ! 

X 

X 

1  ? 

Freon, water 

i 

1  1  X 

not  defined  ! 

dc  !  1 

X 

X 

X 

1 

1  lx! 

1 

X 

7 

!  40,227 

!  11.2-166 

!  0-4.5 

2-5.4 

Mattson 

303 

72  !  U  of  Ml 

x  ! 

X 

SS? 

Freon-113 

i 

horiz 

x?  :  i 

X 

X 

X 

X 

1  ! 

! 

x 

22-66 

E 

1  ? 

1  7.2-24.3 

!  .09-.96 

1.6-5. 6 

Becker,  et  al. 

226 

72  ! Royal  I,  Swe 

x  1 

X 

:  : 

: 

SS 

water 

i 

X 

vert 

I 

dc  :  :  s  x 

x  1  x  !  x 

i 

x  1  ! 

!  i 

1  1 

i  8-272 

I 

!  200-500 

WPICaPM 

i  .14-3.5 

.16-7,6 

Ahead  t  6roeneveld 

66 

72  :  AECL,  Can  1  x 

X 

! 

ni chrome  V 

Freon-12 

1  i 

not  defined 

:  x 

:  :  ix 

x  !  X  1 

1 

!  1 

!  1 

.4-10 

1 

1  155,222 

1  6.2-16.6 

i  0-3.6 

.14-3.3 

■  '  ii  mu  IIMMIS 

26 

72  :u  of  Toronto 

X 

! 

X 

SS 

water  i  x 

1 

i 

vert 

! 

x  :  : 

I 

X  i 

1 

:  i  i 

1 

X 

1.9 

? 

i  2.7? 

!  ? 

i  .19-. 46 

.9-2.3 

113 

73  !  MIT 

X 

X 

nickel 

water  i 

X 

X 

not  defined 

! 

dc  1 

! 

X  1  X 

1 

:  i  : 

x  1 

X 

X 

? 

!D=4.93mm 

!  2. 1-3. 5 

!  0-11 

2. 9-7. 6 

swirl  flow 

Plummer,  et  al. 

304 

73  !  re 

X 

X 

Inconel  X-750 

water 

x  1  vert? 

x?  1 

1  X 

X  1  X 

X 

1 

X 

0 

I 

!  445 

1  6.9-69 

I  0-2.  B 

.07-. 34 

fouling  evaluated,  fill  boiling 

Purcupile,  et  al. 

68 

73  IWestingh/DfJ 

X 

X 

? 

R-U 

vert 

?  ! 

X 

X 

X 

X 

18-28 

I 

!  90 

1  10.B-25 

i  0-1.2 

1.4-8. 8 

Knoebel,  et  al. 

227 

73  ICol  U/du  Pont 

X 

X 

X 

SS, aluminum 

water, D20 

X 

X 

X 

vert 

dc  ! 

X 

X 

X 

X 

X 

X  1  X 

X 

ESa 

L 

!  40-64 

1  2. 1-6.6 

i  0-20.5 

4,6-18.3 

bowing  of  heater  important 

■npmniHrziM 

63 

73  I  USSR 

X 

X 

? 

var  mixes 

not  defined 

-  ,  , 

X 

X 

X 

X 

X 

<  X 

? 

1  10-15 

1  3.3-13.2 

!  0-13.5 

2.5-10 

binary  mixes 

Megerlin,  et  al. 

116 

74  !W6r,Rayth,ISU 

X 

X 

SS 

water 

X 

X 

X 

not  defined 

ac  1 

X 

X 

X 

X 

X 

I  X 

82-203 

I 

!  9,5 

!  2-17.4 

I  0-31.5 

2.1-4 

tube  inserts 

Katto  l>  Monde 

305 

74  !  Japan  1 

X 

X 

Cr-SS 

water 

X 

horiz 

ac  i 

X 

X 

1 

X 

3 

L 

| 

:  l 

i  0-18 

0-60 

jet  impingement 

Hughes 

21 

74  1  Penn  State 

X 

X 

ni chrome 

Freon  113 

horiz 

dc  i 

X 

X 

X 

1  X 

X 

X 

28-61 

I 

i  11 

1  7.6-10.4 

i  0-1.95 

.3-1.3 

straight  and  curved  flow 

Levitan  1  Lantsaan 

223 

75  I  USSR 

X 

X 

X 

7 

water 

X 

not  defined 

?  1 

X 

X 

X 

X 

X 

? 

!  188-375 

!  49-186 

i  0-2.75 

1-3 

Qrnatskiy,  et  al. 

306 

75  !  USSR 

X 

X 

? 

water 

X 

not  defined 

x  ! 

X 

X 

X 

X 

X 

X 

38-337 

I 

1  75-162.5 

i  0-5.0 

,5-1.7 

Ornatskiy,  et  al. 

114 

75  i  USSR 

X 

X 

X 

•> 

water 

X 

vert 

dc  : 

X 

X 

X 

X 

X 

X 

0-159 

1 

5  53-150 

13151 

i  0-7.25 

.5-1.7 

with  l  without  swirl 

Iloeje,  et  al. 

307 

75  !  MIT 

X 

X 

Inconel  X-750 

water 

X 

vert 

x?  : 

X 

X 

X 

X 

0 

1 

i  8.1 

i  69 

i  0-2.75 

.07-.  34 

Kovalev,  et  al. 

43 

76  !  USSR 

X 

X 

copper 

water 

X 

vert? 

•note 

X 

X 

X 

X 

X 

15-70 

1 

l 

1  1 

:  0-10.4 

.02-. 14 

finned  geoms;  electron  bomb  htg 

flerilo 

94 

77  1 Chalk  R.,Can 

X 

X 

X 

1  Ingajel  600 

Freonl2, water 

X 

horiz, vert 

dc  ! 

X 

X 

X 

X 

X 

1 

0-169 

1 

i  97-387 

1  10.5-96.6 

i  0-2.3 

,7-5.7 

Spencer  6  Rousar 

308 

77  !  NASA 

X 

X 

MoneF¥-500,  Inconel 

liq.  oxygen 

vert 

dc  ; 

X 

X 

X 

X 

X 

? 

1  3.6-204 

1  173-345 

i  2-09.7 

!  24-122 

Ornatskiy  V  Sharayevskiv 

248 

77  1  USSR 

X 

X 

SS 

water  ! 

i  :  x 

not  defined 

dc  : 

X 

X 

X 

X 

! 

1  I 

1  X 

VTOKl 

I 

1  13.3 

I  1-4 

!  0-3.5 

1-2 

acoustic  emissions  measured 

Cu*o,  et  al. 

309 

77  !  Italy 

X 

:  x 

SS 

Freon 

vert 

dc  : 

X 

X 

X  i 

! 

1  1 

J 

1 

X 

5-55 

L 

1  20 

!  10.7-21.8 

i  .06-. 28 

.14-. 37 

vapor  thickness  measured 

Lung,  et  al. 

310 

77  !  Germany 

X 

X 

Inconel  600 

water 

X 

vert 

dc  1 

X 

X 

X  i 

1 

1 

1 

10-40 

A 

1  50 

!  0.5-7 

:  0-3 

0.5-8 

Eddington  It  Kenning 

311 

70  1  Oxford  U 

X 

X 

SS 

water 

X 

horiz 

dc  ! 

! 

X 

X 

3-6 

1 

|  ? 

1  1.2 

:  0-.21 

.15-.3 

Del  Valle/Kenning  rig 

Eddington,  et  al. 

312 

70  ! Oxford  U/USSR 

X 

X 

brass 

water 

X 

X 

X 

horiz 

dc  ! 

X 

X 

3 

I 

I  ? 

1  1.2 

!  o-.ll 

.15 

Del  Valle/Kenning  rig 

Oieng,  et  al. 

124 

78  1  U  of  Ottawa 

X 

X 

x 

copper 

water 

X 

vert 

X 

X 

X 

X 

0-28 

I 

!  8.5 

!  1 

i  0-2 

.068-.  21 

Kutateladze  It  Malenkov 

313 

79  !  USSR 

X 

X 

■icroporous  matl 

var .water 

X 

vert? 

X 

X 

X 

0 

L 

i  ? 

1  ? 

i  0 

? 

bubble  stability 

Monde  l>  Katto 

314 

78  !  Japan 

X 

X 

X 

copper 

Freon, water 

X? 

vert 

X 

X 

X 

X 

X 

0-30 

L 

! 

:  l 

:  0-19.8 

2-26 

impinging  jet 

Yucel  It  Kakac 

140 

70  !  Turkey 

X 

X 

SS 

water 

x 

horiz, vert 

dc  1 

X 

X 

X 

X 

!  x 

X 

3-42 

I 

1  14-28 

:  t 

i  0-6 

1.25-6.25 

varied  heater  orientation 

Drizius,  et  al. 

315 

73  1  USSR 

X 

X 

copper 

water 

X 

not  defined 

dc  i 

X 

X 

X 

X 

X 

X 

25-70 

L 

!  40-153 

!  3.9-13.4 

!  0-81.3 

“FS - 

swirl  flow 

Ooroschuk,  et  al. 

316 

70  1  USSR 

X 

X 

x 

X 

7 

var 

X 

vert 

?  i 

X 

X 

X 

X 

X 

X 

? 

!  225,? 

1  20-167 

!  0-1.2 

.75-3 

Ragheb,  et  al  . 

317 

70  1  U  of  Ottawa 

X 

x 

X 

copper 

water 

X 

not  defined 

X 

X 

X 

X 

X 

0-28 

I 

1  6.3 

;  i 

i  0-1 

.034-.  1 

transition  studied 

Mel 'nikov,  et  al. 

107 

79  !  USSR 

X 

X 

X 

SS 

water 

X 

vert 

?  ! 

X 

X  1  X 

:  i 

X 

0-86 

E 

I  ? 

1  30 

!  0-5.3 

2. 5-3.5 

fouling  studied 

Chen,  vt  al. 

319 

79  5  U  of  Ottawa 

X 

X 

SS, Inconel 

water 

X 

not  defined 

ac  ! 

X 

X 

X 

1 

X  !  x 

X 

X 

10-80 

I 

1  202.253 

:  i 

!  0-2.6 

.1-4 

Romberg  It  Harris 

246 

79  !  Australia 

X? 

X 

SS 

Freon  113 

not  defined 

i? 

1 

1 

X 

? 

1  ? 

:  ? 

!  ? 

? 

acoustic  meas. 

Ragheb  It  Cheng 

319 

79  !  U  of  Ottawa 

X 

X 

I 

copper 

water 

X 

vert? 

X 

X 

X 

X 

X 

0-28 

I 

1  8.3 

i  i 

i  0-1? 

.034-. 1 

Bennett  It  Chen 

320 

80  !  Lehigh  U 

X 

X 

Inconel 

glycol .water 

X 

vert 

x  ! 

X 

X 

0 

I 

1  3 

i  ? 

i  .7-3 

.16-1.6 

pure  and  binary  coolant  mixture 

Mori,  et  al. 

321 

80  1  Japan 

X 

X 

X 

nickel, glass 

water/oil 

:  i 

X 

horiz 

X 

dc  1 

X 

? 

!  75 

:  l 

i  0-1.5 

.005-.  007 

bubble  generation 

Petukhov,  et  al. 

322 

X 

X 

SS 

1 i o  heliui 

! 

vert 

dc  ! 

X 

X 

1  i 

0? 

I 

!  300 

!  1 

i  0-.004 

1  .015-.33 

Koba,  et  al. 

323 

SO  !  USSR  1  x  1  x 

X 

steel 

water 

X 

I 

horiz? 

dc  ! 

X 

X 

X 

X 

X  !  ! 

I3PM 

1 

i  15-200 

1  HI 

:  4-21 

.5-29.9 

| 

small  diameter  tubes 

Tolubinskiy,  et  al. 

250 

80  1  USSR  1  x  i 

X 

? 

water 

1  X 

not  defined 

dc  : 

X 

X 

X  1  X 

i 

I  1 

X 

? 

j  ? 

1  1-5 

:  0-3 

.5-2 

l 

acoustic  measurements 

Yilmaz  V  Hestwater 

324 

80  1HTRI/U  of  IL  !  J  x 

X 

: 

cooper 

Freon-113 

i 

; 

1 

vert 

! 

X 

i  X 

1  1 

I  1 

! 

1 

X 

I  0-5 

L 

1  20 

1  1-1.1 

!  0-.75 

0-6.7 

I 

Ornatskiy,  et  al. 

119 

80  i  USSR  i  x  1 

X 

! 

? 

water  ! 

1  1 

X 

vert 

! 

x?  5 

X 

X 

X  [  X  1 

1  X 

i  1 

1 

X 

l 

!  9-374 

I 

1  25-300 

!  73.6-201 

:  0-6.8 

i  .5-3.5 

Jain,  et  al. 

156 

80  !  U  of  IL  i  x  1 

X 

i 

SS 

R-113  1 

! 

vert 

f 

ac  1 

X  1  X 

X  1  1 

i  i 

1  X 

x  !  10-60 

I 

l  124 

!  3.5-4. 1 

i  .07-.  15 

J  .57-.  7 

void  fraction  measurements 

Ornatskiy  l  Sharayevskiy 

249 

80  1  USSR  !  x  1 

X 

l 

SS 

water  ! 

1 

X 

not  defined 

i 

dc  1 

X 

X 

X 

!  1 

1 

X 

10-120 

I 

!gap=.5M 

:  1-4 

!  ? 

1-2.5 

1 

acoustic  measurements 

Fisenko,  et  al. 

99 

80  :  USSR  lx 

X 

SS 

water  ! 

X 

x  ! 

vert 

i 

x?  1 

x 

X 

X 

X 

1 

10-100 

I 

1  125 

1  100-180 

i  0-2.9 

2 

Leontiev,  et  al. 

121 

81  !  USSR 

X 

X 

X 

SS 

water  ! 

X 

! 

horiz  ! 

ac  1 

X 

X 

X 

X 

X 

0-145 

L 

1  69 

i  68.7-137 

!  0-7.25 

.75-2 

Zeigarnik,  et  al. 

96 

81  \  U5SR 

X 

X 

copper 

water  ! 

X 

! 

horiz 

ac  1  1 

X 

X 

X 

X 

1  X 

20-160 

I 

1  54 

i  5-30 

I  8-62,1 

4-25 

Yao 

325 

81  ICarnegi e-Mel 1 

X 

X 

X 

l 

SS 

var, water  i  x 

X 

1 

1  vert 

dc  l 

X 

X 

0-4.5 

I 

I  7-136 

i  1-1.1 

i  O-.l 

.125-.2 

.19-.3 

gaps 

Styrikovich,  et  al. 

326 

81  1  USSR 

X 

X 

X 

SS 

water 

X1/ 

horiz 

ac  ! 

X 

X 

X 

X  1  ; 

X 

1 

0-116 

L 

1  75,152 

nisstij 

:  0-4.5 

1,2 

Ragheb,  et  al. 

123 

81  !  U  of  Ottawa 

X 

X 

X 

Cu,!nc,Al,Zircaloy 

water 

X 

vert? 

X 

X 

X 

X 

x  ! 

X 

0-28 

1 

!t*I2  h 

:  l 

!  0-2.5 

.07-. 2 

ss  k  transient  heating 

Ramamohan,  et  al. 

327 

81  1  India 

X 

X 

platinum 

gly,  ale /water 

X 

vert? 

x?  ! 

X 

X 

40-64 

A 

!  150 

j  7 

:  .08-2.4 

.  15-.55 

Jensen  L  Bergles 

95 

61  !U  of  WI,IA  St 

X 

X 

X 

SS 

R-113 

vert 

dc  ! 

X 

X 

X 

X 

X 

X 

0-110 

1 

1  83-174 

;  9.4 

i  .O5-.0 

.57-5.5 

straight  V  coiled  flow 

France,  et  al. 

328 

81  lArgonne  N  L 

X 

X 

chromoly  steel 

water 

X 

vert 

i  X 

X 

X 

t 

? 

1  1297 

i  70-153 

I  .5-1.6 

.7-3.2 

Agrawal,  et  al. 

329 

82  !  India 

X 

X 

SS 

R-12 

horiz 

ac  l 

X 

X 

0 

I 

!  210 

!  ? 

:  .008-.02 

.2-. 39 

plain  1  swirl  flow 

Aounallah,  et  al. 

330 

02  !  Oxford  U  ! 

X 

X 

cupro-nickel 

water 

X 

vert 

ac  ! 

1 

t 

X 

0 

I 

!  48 

!  1-1.4 

1  ? 

•  1,.2 

Bel  Valle  l>  Kenning 

81 

02  !  Oxford  U 

X 

X 

SS 

water 

X 

X 

horiz 

dc  ! 

X 

1  X 

X 

X 

EEa 

I 

1  21 

i  1 

i  .3-4.7 

.8-2 

Eddington  et  al  rig 

Fukuyama  It  Hirata 

331 

82  I  Japan 

X 

X 

SS 

Freon  113 

horiz 

ac  ! 

X 

X 

X 

X 

I 

i  1 

I 

I 

? 

1  51 

i  3.5-14.9 

1  0-9 

41.5-52.9 

questionable  mass  flux 

Katto  t  Ashida 

332 

82  !  Japan 

X 

X 

X 

SS 

R-12 

vert 

dc  I  1  I  x 

X 

X 

X  1 

I 

:  : 

1 

i 

!  0-10 

I 

!  50 

i  19.6-34.3  !  0-.83 

! 

.7-7 

Leung,  et  al. 

44 

k  :  *cl,cw/uca 

X 

X 

l 

Inconel  600 

water 

! 

X 

vert 

!  x 

x  1 

!  x 

X 

X 

| 

!  i 

1 

i  i 

x  :  0-106 

I 

!  1162 

!  97 

i  0-1.25 

2-5 

different  heater  configurations 

Nishima  It  Ishii 

333 

B2  !  Japan/ANL 

X 

X 

X  1 

l 

1  SS 

water 

X 

! 

vert 

1 

dc  1 

(  X 

X 

X 

! 

i  1  ! 

! 

[  ? 

!  108 

:  l 

I  ? 

.005-. 035 

Nishikawa,  et  a). 

190 

82  I  Japan 

X 

X 

X 

:  : 

SS 

Freon  22  i 

!  vert 

1 

ac  S 

X 

X  1  X 

X 

i  i 

j 

0-24 

I 

!  154 

!  34-48 

:  0-.09 

1  .2-1.3 

Papaioannou  t  Koumoutsos 

61 

82  I  Greece 

7 

X 

bronze 

water/pol ymer 

X 

horiz 

X 

X 

:  i 

! 

X 

X 

? 

! 

!  ? 

\  0-.02S 

O-.l  1 

bubble  generation 

Rogers,  et  al. 

240 

82  i  U  of  Ottawa 

X 

X 

X 

Inconel  710 

water 

X 

vert 

x?  i 

X 

X 

X 

x  !  1  x 

1 

X 

X 

33-93 

I 

!  28-54 

i  1.6 

!  .15-2.6 

.06-1.2 

void  fraction  measurements 

Katto  It  Yokoya 

32 

82  !  Japan 

X 

X 

X 

SS 

R-12 

vert 

dc  i 

X 

X 

X 

X 

x  !  !  x 

X 

0-16 

I 

!  200,333 

!  19.6-34.3 

i  0-.5 

1.1-9 

Sekoguchi,  et  al. 

334 

82  !  Japan 

X 

X 

X 

? 

water 

X 

X 

vert 

dc  ! 

X 

X 

X 

x  1  !  x 

0-84 

7 

!  91-154 

!  2-16 

1  .12-1.86 

.3-2 

Viannay  t  Kan  an 

335 

82  1  France 

X 

X 

X 

steel .copper 

var, water 

X 

not  defined 

X 

X 

X 

X 

t  i 

X 

7 

!  ? 

!  ? 

i  .HO 

.7-11 

Wei  k  Maa 

62 

82  !  Taiwan 

X 

SS, copper 

water/polymer 

X 

vert 

dc  1 

X 

x 

X 

x  ;  : 

X 

X 

X 

7 

1  10.5-32 

j  ? 

!  0-.9 

.03-. 3 

surface  tension  effect 

lleda,  et  al. 

336 

83  !  Japan 

X 

X 

X 

copper 

Freon-113 

vert 

I 

X 

X 

X 

x  1  1 

x«-.06to.6  I 

i  6.10 

:  2,3 

i  0-.55 

.36-1.5 

Sharon,  et  al. 

337 

83  ! Northwestern 

X 

X 

X 

SS 

water 

X 

vert 

X 

X 

X 

X 

1  !  x 

X 

X 

0-15 

I 

!  60,100 

:  1.3-7 

!  0-.5 

.33, .57 

Bilicki 

55 

83  !  Poland 

x 

X 

SS 

Freon  21 

horiz, vert 

x?  ; 

X 

X 

J  J 

1  0 

I 

!  38 

:  6 

1  O-.ll 

.62-1.2 

Weede  It  Dhir 

145 

83  !  UCLA 

X 

X 

Inconel  600 

Freon  113 

vert 

dc  ; 

X 

X 

X  1 

x  :  : 

1 

X 

x»-.2to-.3  1 

!  5-26.5 

:  2 

i  0-.B4 

2. 1-6.3 

with/without  tangntl  injection 

Klimeako  It  Sudarchikov 

338 

83  1  USSR 

X 

X 

X 

SS 

liq  nitrogen 

vert 

ac  : 

X 

X  1  X 

1  1 

i  ; 

1 

0.2-10 

1 

!  185 

!  1.8-7 

1  .01-. 04 

.17-. 75 

large  length 

Kataoka,  et  al. 

125 

83  !  Japan 

X 

X 

X 

platinum 

water 

X 

X 

vert 

dc  1 

x  !  x  !  x 

X  1 

1  X 

1 

I 

1  X 

0-70 

I 

!  26.2-126  !  1,4-15 

i  0-19.5 

.4-4 

transient  k  s.s.  heating 

Kalayda,  et  al. 

103 

B3  !  USSR 

X 

X 

SS 

water 

X 

X 

vert? 

ac  1 

!  x  !  x 

X  1  X 

:  i 

1 

0-189 

I 

1  .24? 

I  80-160 

!  t-2.55 

.5-2 

Baines,  et  al  . 

339 

84  1  MIT  lx 

X 

X 

1  copper 

water .Freon 

i 

1 

X 

inclined 

X 

1  1 

x  : 

X  1  X  1 

i  : 

1 

:  5-i8 

I 

| 

:  i 

!  0-1.62 

1-5 

gravity  &  jet  driven  films 

Katto  l<  Dhno 

189 

84  !  Japan 

X 

x  ; 

SS 

r-i: 

1 

vert 

dc  1 

x  !  x 

X  t  X 

i 

.1-10 

1 

1  100 

!  19.6-34.4 

1  0-.2 

.12-2.1 

Kozhelupenko,  et  al. 

340 

85  !  USSR 

X 

:  x  : 

5S 

water  i  x 

! 

horiz? 

I 

x  ! 

x  :  x 

X  1  X 

1 

1  X 

X  l 

X  1 

0-80 

I 

!  33-50 

!  HO 

!  .05-1.6 

.02-.  7 

Hina  It  lleda 

341 

B5  !  Japan 

x 

!  x  : 

SS 

Freon  113  ! 

i 

vert 

! 

ac  ! 

1  X 

i  X  i  X 

1 

1 

1 

X 

X  ! 

10-30 

I 

1  40 

!  1.5 

!  0-.2 

.51-1.24 

fluid  tesp  probe  survey 

145 
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97 

343 
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347 

Yr 

85 

85 

86 
86 
96 
86 
66 
86 
56 
86 

Organization 

TYJ 

Subcool 

(surf) 

>E 

Sat 

(bulk) 

Tube 

X 

X? 

Annulus 

Rod  / 
Hire 

SURF  AC 
Tube 
Bundle 

£ 

Plate, 

Channel 

Ribbon 

Matl 

Typ. 

dis 

X 

dei 

COC 

T 

dea 

1ANT 

eatat 

deg 

X 

jnt  I 
gas 

3 

untr 

nwtd 

X 

X 

Flow 

Orientation 

Coal,6as 

Fired 

Elect 

Htr 

X? 

HEATING 

Elect 

Resis 

dc 

dc 

X 

ac 

7 

X 

ac 

Liquid 

Metal 

Steaa, 

Fluid 

dTsub 

123 

X 

Press 

X 

X 

X 

X 

Vel 

133 

X 

X 

X 

X 

X 

X 

X 

Burnout 
crit  pt 

EFFECTS 
6as  !  Accel 
Contnt  ! 

TIDIED 

L/D  i  Hall 
[43  1  Matl 

i  X 

Dh 

153 

Hall  1  Rough- 
Thk  i  ness 

Press 

Drop 

Ihi/tawi 
Ht  Flux 

Stability 

Transient 

Optical 

Other 

x ' 

X 

X 

X 

dTsub  I  L/D 

deg-C  (23  !  143 

1  Press 
!  /  bar 

PARAMET 
!  qdot  C61 

1  MW/sq-a 

ER  RANGE 

Vel 

«/5 

adot  (73 
kg/s 

!  6  (91 

!.  Mg/a2  s 

Coaaents 

Evans,  et  al. 

Edelmam,  ft  al. 

Boyd,  et  al. 

Fujita,  et  al. 

Fukuyama,  ct  al. 

Grigoriev,  et  al. 
Muller-Steinhagen,  et  al. 
Tolubinskiy  b  Doeashev 
Celata,  et  ah 

Ueda  b  Kim 

Lehigh  U 
Technicn,  Isr 
Sandia  NL 
Japan 

Japan 

USSR 

Canada 

USSR 

Italy 

Japan  /Korea 

X 

X 

X 

X 

X 

X 

X 

X 

X? 

X 

X 

X 

Inconel  600 

SS 

cooper 

copper 

SS 

SS,  copper 

SS 

7 

SS 

SS 

water 

water 

water 

Freon  113 
Freon  113 
liq  nitrogen 
heptane 
water? 

R-12 

Freon  113 

vert 

vert 

horiz 

horiz 

horiz 

horiz, vert 
horiz,  vert 
vert? 
vert 

vert 

X 

X 

X 

X 

1 

X 

x  ! 

i  X 

i  X 

1 

X 

V 

( ,  •' 

. 

0 

78-83 

1B6 

59-140 

? 

92-139 

? 

23 

10-30 

1  !  88 
l  1  167,169 

I  1  116 

1  :  4.8 

L  :  50 

1  ? 

L  !  7,625 
t  100 

I  1  307 
i  :  40 

1  2. 4-5. 7 
!  1 

1  16 
!  1-5 

!  5.8-23.8 
!  1-20 
!  1-7.9 
!  98-147 
!  12.5-27.5 
:  1.5 

i  .02-.  06 
!  .03-.  24 

I  1.5-4.25 
i  0-.1 
!  0-10.1 
i  .006-. 15 
f  0-.4 
!  0-5.5 

I  .03-. 085 
!  0-.45 

S  .01 3-. 085 
!  .08-.25 

1  .63-3.5 
!  ? 

1  28-56 

!  .005-.7 

1  1 

:  .4-1.6 
!  .26-1.2 

vapor  probe  aeasureaents,^ 
void  fraction  measurements 

questionable  aass  flux 

heater  geometry  varied 

CHF  detector 
flow  transients 

Aounallah  l  Kenning 

348 

8/ 

Oxford  U 

X 

X 

cupronickel 

water 

X 

vert 

dc 

X 

X 

! 

i 

5-15 

?  :  52 

i  i.i 

1  0-.17  1 

.1-.2 

Yezhova,  et  al. 

349 

87 

X 

X 

1 

? 

water 

X 

x 

not  defined 

X 

! 

X 

X 

X 

X  I  X  i 

i 

0-85 

?  !D=8  aa 

!  98.1-167 

:  0-5.3  i 

1-3.7 

lee  b  Shen 

350 

87 

X 

X 

X 

SS 

water 

X 

not  defined 

1  ? 

X 

X 

1 

1  X 

I 

7 

!gap=12  aa 

j  ? 

!  0-1.5 

.05-2 

falling  files  also 

Inasaka  l  Nariai 

351 

87 

Japan 

X 

X 

SS 

water 

X 

vert 

i  dc 

x  : 

X  1  X 

! 

x  I  1  X 

i  i  X  1 

i 

l  0-75 

?  I  3-100 

!  1 

!  ? 

7-20 

EI-6enk,  et  al. 

352 

80 

U  of  Ml 

X 

X 

SS 

water 

X 

vert 

!  dc 

X  1 

X  1  X 

! 

:  !  x 

:  : 

44-75 

I  !  39-68 

:  1.2 

!  .16-1.4 

0-.26 

Muller-Steinhagen,  et  al. 

353 

88 

N  Zeal, Can 

X 

X 

SS 

H20, heptane 

X 

X 

vert 

X 

x  !  x 

X  1 

x  i 

:  : 

i 

32-151? 

L  1  7 

I  1-7.3 

(  0-.37 

.2-.  76 

Yin,  et  al. 

354 

83 

Can, Taiwan 

X 

X 

Inconel  600 

water 

x 

vert 

dc 

X 

X 

X 

i 

10-80 

I  !  273 

I  10-214 

!  0-2 

2 

Fleishman  b  Yuen 

355 

88 

0  of  CA 

X? 

X 

copper 

water 

X 

horiz? 

X 

X 

! 

X 

X 

? 

*  ? 

}  7 

i  . 15-.6 

.0018-.  12 

radial  flow  thru  wire  aesh 

Passes  b  Gentile 

356 

88 

France 

X 

X 

Inconel 

Freon-113 

vert 

X? 

X 

X 

X 

X 

1 

5-40 

1  I  6.25 

1  1.5-4 

!  0-.6 

.4-1.3 

transition  evaluation 

Lu  b  Jia 

177 

88 

China 

X 

X 

thin  copper  layer 

water 

X 

X 

X 

hori: 

x 

X 

X 

! 

X 

X 

5-12? 

I  ’  72.3 

!  1. 8-2.4 

I  .058-.  077 

.82-1.29 

void  fraction  measurement 

Blujcenkrantz 

357 

80 

Argentina 

X 

X 

X 

7 

Aezirc® 

water 

X 

? 

X 

X 

X 

X 

X 

! 

0-160 

1  !  232 

!  68-151 

i  ? 

.,4-4.65 

Boyd 

358 

83 

FV  A4M  U. 

X 

X 

water 

X 

X 

horiz 

dc 

X 

X  I  I 

1 

X 

0-74 

E  {  96.6 

!  7.7 

!  0-42 

4.4-41 

Lee  b  Si wan 

359 

89 

Motorola, IMf 

X 

X 

nichrome, platinum 

fluoro  FC-72 

horiz 

dc 

X 

X 

X 

X 

X  I 

i  X 

13-46 

1  !  .05-.  65 

i  1.3-2. 0 

!  0-1.7 

2. 5-8. 3 

Celata,  et  al. 

129 

89 

Italy 

X 

X 

X 

SS 

R-12 

vert 

X 

X 

X 

X 

X 

X  ■ 

i 

X 

!  0-23 

1  (  153-299 

1  12.5-27.5 

i  0-.12 

1-1.5 

transient  &  s.s.  heating 

Boyd 

174 

89 

PV  AbH  U. 

X 

X 

Aazirc 

water 

x 

X 

horiz 

dc 

X 

X 

! 

X 

183 

E  i  96.6 

:  16.6 

:  0-42 

4.4-32 

Jung,  et  al. 

360 

89 

UHD.NBS 

X 

X 

X 

SS 

refrig  aixs 

horiz 

dc 

X 

X 

1 

•upto  x*,95  E  !  444 

i  2.6-4 

!  .01— ,045 

KaBSa 

.25-.72 

Jung,  et  al. 

361 

89 

HID, NIST 

X 

X 

5S 

refrig  aixs 

horiz 

dc 

X 

1 

up  to  ««.?  E  !  444 

!  pr=0.08 

1  .01-. 045 

.016-.  045 

.25-.72 

follow-on  to  previous  expt 

6u,  et  al. 

110 

89 

U  of  ky,  m 

X 

X 

copper 

fluoro  FC-72 

horiz 

X 

X 

X 

X 

!  x 

i 

,5-38 

L  !  1 

:  i 

!  0-1.4 

1-4 

straight  b  curved  flow 

Polonskii,  et  al. 

362 

90 

USSR 

7 

X 

SS 

water 

X 

i  horiz 

X? 

X 

X 

!  1 

X 

? 

1  ? 

\  100-140 

1  .2-1.3 

j 

.5-2 

Hasan,  et  al,. 

363 

90 

AZ  State.EPRl 

X 

X 

SS 

Freon-113  1 

1 

vert  1 

1  dc 

! 

x  :  x 

X 

:  i 

1 

1  20-37 

I  !  121 

!  2.4-3. 1 

•01—. 15  i 

.58-1.1 

364 

90 

UK 

X 

X 

copper 

Freon-113  J  ! 

s 

vert 

t 

!  1 

!  ! 

i 

X 

0 

L  1  ? 

1  1 

i  0-. 039 

.1 

Sivagnanai  b  Varaa 

365 

90 

India 

X 

X 

Pt 

var /water  {  1 

S 

vert 

dc 

I 

x  I 

X  1 

! 

1  ! 

! 

1 

10-40 

A  !  1650 

1  ? 

i  .2-3 

.15-1 

binary  aixtures 

EI-6enk,  et  al. 

366 

90 

U  of  HI 

7 

1  i  X 

SS 

water  1  ! 

X 

vert 

dc 

t  1 

x  : 

! 

:  i 

:  : 

!  x 

7 

1  71 

!  ? 

i  ? 

.013-1.52 

Shi«,  et  ah 

49 

90 

Canada/Korea 

x  : 

X 

SS,SS/aluainua 

water 

X 

vert 

X? 

1  X 

X 

X  i  X 

1 

i  : 

i  x  : 

: 

42-120 

I  !  75-94 

!  1.1-3. 5 

0-11.3  ! 

1-5.9 

finned  b  unfinned  geometries 

Dietrich,  et  ah 

367 

90 

New  Zealand 

X 

X 

X 

SS 

isobutanol 

vert 

dc 

!  x 

X  1 

i  i  x  \ 

X 

? 

l 

:  i 

:  0-.5 

.01-.4 

coiled  wires 

Vandervort,  et  ah 

368 

90 

RPi 

X 

X 

SS,  Inconel 

water 

X 

X 

X 

vert? 

dc 

l  X 

X 

X  '  X 

X  r 

X 

!  I  i 

X 

40-135 

E  !  2-25 

:  6,  12 

:  io-ioo 

5-40 

many  preaature  failures 

Yuncu,  et  ah 

369 

90 

U  of  Hiaai 

X 

X 

SS,  nichroae 

Freon-11 

vert 

X 

i 

I 

! 

X 

X 

X 

? 

!  81-86? 

:  l 

!  0-.02B 

0-.035 

surface  enhancements 

Erauer,  et  ah 

370 

90 

Gera  any 

X 

X 

SS 

R-12 

vert 

dc 

x 

X 

X 

X 

X 

10-75 

1  :  46,4 

!  10-33 

0-.12 

.5-3 

void  fraction  measurement 

Hosaka,  et  ah 

371 

90 

Japan 

X 

X 

SS 

Freon-113 

vert 

ac 

X 

X 

X 

X 

X 

50- B0 

?  !  50 

!  11-24 

!  0-6 

9.3-32 

Hihara  b  Saito 

372 

90 

Japan 

X 

X 

55 

R22.fttl4.ftix 

horiz 

dc 

X 

0 

OB 

:  6 

.008-. 026 

.2-.  35 

Hahne, .et  ah 

373 

90 

Germany 

X 

X 

copper 

R-12 

vert 

X 

X 

X 

X 

2.2-10.0  I  1  25 

!  7-8.1 

!  0— .045 

.2-2.25 

Dougherty,  et  al. 

374 

90 

Col  uadi  a  U. 

X 

X 

SS, inconel 

water 

x 

X 

vert 

dc 

1 

X 

X 

X 

x  ! 

X 

X 

X 

60-110 

I  !  85-158 

!  1.38,3.45 

1.26-3.16 

0-2.5 

liu  b  Kakac 

375 

91 

U  of  hiaai 

X 

X 

nichroae 

R-ll 

vert 

dc 

x 

x 

X 

7 

!  80.7 

•  7 

.022-. 055 

.007-. 02 

Ravindran,  et  ah 

376 

91 

India 

X 

X 

X 

SS 

water 

X 

horiz 

ac 

X 

X 

K 

X 

4-10 

I  !  27-92 

;  i.i,  i.s 

1  .0085-. 034 

.017-, 38 

Boyd 

377 

91 

PV  AIM  U. 

X 

X 

copper 

water 

X 

X 

horiz 

dc 

1 

X 

x  ! 

1 

X 

129 

1  i  49 

!  4.5 

0-11 

1.56-8.55 

Feng  b  Johann  sen 

378 

91 

6eraany 

X 

X 

tonel 

water 

X 

X 

X 

vert 

X 

X 

X 

X 

X 

i 

X 

5-57 

1  5  5.2 

l  1.1-12 

!  .02-3.3 

.01-1 

transition  boiling 

Celata,  et  ah 

379 

91 

Italy 

X 

X 

X 

SS 

R-12- 

vert 

X 

X 

X 

X 

X 

x  ! 

X 

X 

0-23 

1  1  153-299 

!  12,5-27.5 

!  0-.12 

; 

1-1.5 

trans  b  s.s.  heating, press, flow 

Kocamustifaogullari  b  Hang 

380 

91 

U  of  HI 

X 

- 

water 

X 

horiz 

X 

X 

? 

1  306 

I  ? 

I 

3.75-3.7 

void  fract  aeas,  unheated  tube 

Hasan,  Roy  b  Kalra 

381 

91  1AZ  St  U/EPRI 

X  1 

X  i 

1 

SS 

R-113 

! 

i 

1 

vert 

dc 

X 

X 

X 

! 

x  !  x  !  20-37 

I  i  124 

1  1.8-2. 5 

•  .04-.  14  i 

.50-. 9 

vapor  characteristics  measured 

Balino  b  Canverti 

302 

91 

Argentina 

!  x 

!  X 

SS 

Freon-12 

! 

1 

horiz 

dc 

X 

X 

i 

:  i 

:  : 

X 

0 

I  !  23.4 

!  10.3 

0-.32 

.24-1.6 

eccentricity  effects 

Akhanda  b  Janes 

51 

91 

Bangladesh/UK 

x  ! 

X 

SS 

water 

X 

X 

X 

1 

■  hori: 

ac 

X 

X 

i  X 

X 

5-30 

I  !  13.8 

:  i 

1  0-.55 

.2-1.1 

finned  b  unfinned  geometries 

Boyd 

383 

91 

PV  AIM  U. 

X  1 

X 

zirconium-copper 

witer  (j) 

X 

X 

horiz  ! 

i  dc  ! 

X  1  X 

i  i 

188 

I  1  96.6 

!  16.6 

i  0-36 

4.4-31.5 

Gyulmamedov,  et  ah 

3S4 

91 

USSR 

X 

x  !  x 

1  SS 

water/TrilonF 

X 

X 

vert 

1  ac 

X 

X 

X 

X 

:  : 

X 

x*-.7to.7  I  !  150 

!  8-160 

1  0-2.4? 

.1-1.5 

chelating  agent  effect 

6alloway  b  Mudawar 

385 

92 

Purdue 

X 

X 

X 

copper 

fluo-o  FC-72 

X 

circuafr 

x  ! 

X 

X 

X 

X 

!  x 

x  ! 

i  ; 

X 

35-50 

L  !  1-4 

(  1-1.7 

i  0-1.65 

7 

centrifugal  force  effect 

Dykas  b  Jensen 

386 

92 

RPI 

X 

X 

X 

SS 

Freon-113 

X 

vert 

x  1  dc 

1 

X 

X 

X 

X 

1 

x=-.04to.7  L  1  10.4 

I  1.5,  5 

!  0-.33 

.05-5 

leroux  b  Jensen 

387 

92 

RFT 

X 

X 

X 

X 

SS 

Frecn-113 

X 

vert 

dc 

X 

X 

X 

X 

0-8 

1  !  10.4 

i  1.5,  5 

i  0-.33 

.05-. 5 

Celata,  et  ah 

132 

92 

Italy 

X 

X 

X 

SS 

R-12 

vert 

X 

X 

X 

X 

X 

x  i 

X 

0-23 

1  !  153,298 

!  12-27.7 

:  .026-.  12 

.4-1.6 

trans  b  s.s.  power, press, flow 

Ckozco  b  Chen 

389 

92 

U  of  IL 

7 

sphere 

copper 

Freon-113 

vert 

X 

X 

X 

X 

i  X 

X 

X 

X 

? 

! 

!  1 

i  0-.55 

I 

transient  incl.  film  boiling 

Sadasivan  b  Lienhard 

60 

92 

U  of  Houston 

X 

X 

X 

SS,  Nichroae 

var, water 

X 

vert 

ac,dc 

X 

X 

:  x 

X 

X 

0-5 

L  !  33-267 

1  ? 

:  0-7.5 

.1-3 

acceleration,  wettability  effet 

Panchal,  et  ah 

389 

92 

ArgonneNL,U!L 

X 

X 

X 

aluminum 

R-ll 

vert 

X 

X 

? 

!  272 

!  ? 

1  ? 

? 

spirally  fluted  b  plain  tubes 

Willingham  l  Hudawar 

390 

92 

Purdue 

X 

X 

electronic  chip 

FC-72 

X 

horiz 

X 

x 

x 

X 

X 

X 

3-36 

I  1  3.6-13.3 

:  1.36 

:  o-.b 

.13-4 

flow  over  electronic  chips 

Y and er vort,  et  ah 

391 

92 

RPI 

X 

X 

SS 

water 

X 

X 

vert? 

dc 

X 

x 

X 

X 

x  ! 

X 

X 

0-150 

£  !  1-100 

1  1-22 

!  0-104 

5-50 

Vandervort,  et  ah 

80 

92 

RPi 

X 

X 

SS.Ni, brass, Inconel 

water 

X 

X 

vert 

dc 

X 

■ 

■ 

X 

X  ! 

X  1  X 

X 

X 

X 

X 

10.5-142  E  !  1.66-49. 

1  1.31-22.0  1  0-123.0 

3-42.7 

void  fraction  measurement 

NOTES! 


Ill  coolant  treatment  dis 

distilled 

(63  adot 

heat  flux 

dei 

deionized 

dem 

demineralized 

[73  mdot 

mass  flux 

deg 

degassed 

gas 

gassed 

183  G 

■ass  velocity 

untr 

untreated 

nwtd 

no  water  treatment  discussed 

•  193 

where  a  question  mark  directly  follows  an  item  or  value  it  is  uncertain  but  probable 

121  dTsub  subcooling,  Tsat— Tbul k 

!  inlet,  L  local,  E  exit,  A  average  (of  inlet  b  outlet),  ?  unknown  location  and/or  value 
(in  some  cases  listed  as  inlet,  the  inlet  bulk  temp  b  exit  pressure  were  used  to  compute  dTsub) 


131  Vel  velocity 

141  L/D  for  tube  D  =  Dh 

for  circular  annulus  D  =  Dh  =  2(a-b)  a  =  outer  radius,  b  *  inner  radius 
for  rect  duct  or  channel  D  -  Dh  =  2ab/la+b)  a  =  width,  b  =  height 
for  ribbon  D  =  H,  W  =  width 

C5J  Dh  hydraulic  diaaeter  (see  note  4) 
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AEDC-TR-93-3 


NOMENCLATURE 

A  Surface  area,  ft2 

a  Centrifugal  acceleration,  ft/hr2 

Bo  Boiling  number,  =  q"/G  ifg 

c,  cp  Specific  heat,  Btu/lbm  °F 

CD  Drag  coefficient  for  a  sphere 

CHF  Critical  heat  flux,  Btu/ft2  hr 

D  Diameter,  ft 

d  Bubble  departure  diameter,  in. 

DNB  Departure  from  nucleate  boiling 

erfc  Complementary  error  function 

f  Friction  factor  as  given  in  Eq.  (5) 

G  Mass  velocity,  =  pV  except  as  defined  in  Eq.  (12),  lbm/ft2  hr 

g  Acceleration  of  gravity,  ft/hr2 

gc  Gravitational  constant,  ft/hr2 

h  Heat-transfer  coefficient,  Btu/ft2  hr  °F 

i  Enthalpy,  Btu/lbm 

ifg  Latent  heat  of  vaporization,  Btu/lbm 

J  Energy  conversion  factor,  =  778  ft-lbf/Btu 

k  Thermal  conductivity,  Btu/ft  hr  °F 
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AEDC-TR-93-3 


L  Length,  ft 

mv  Vapor  rate,  lbm/hr 

Nu  Nusselt  number,  =  hD/k 

p  Pressure,  psia 

Pr  Prandtl  number,  =  juc/k 

q",q/A  Heat  flux,  Btu/ft2  hr,  except  as  noted 

qVV  Internal  heat  generation  rate  per  unit  volume 

R  Ideal  gas  constant 

r  Characteristic  radius,  ft 

R'  Laplace  number,  =  r[g(gf  -  qv)/o]W2 

Re  Reynolds  number,  =  qY D/fi 

T  Temperature,  °F 

t  Wall  thickness  (Fig.  35) 

t^2  Time  to  reach  half  of  initial  flow,  sec 

ttr  Flow  rate  transit  time,  sec 

V  Velocity,  ft/sec 

x  Vapor  quality 

Z  Length,  as  defined  in  Figs.  21  and  24 

SYMBOLS 

a  Thermal  diffusivity,  ft2/hr 

AC  Excess  concentration  of  the  volatile  component  in  the  vapor,  percent  (see  Fig.  17) 
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Aij  Enthalpy  difference  due  to  inlet  subcooling,  =  iiocal  -  ijn,  Btu/lbm 

ATbuik  Complete  temperature  difference,  =  Tw  -  Tb,  °F 

ATsat  Superheat,  —  Tvv  -  Tsat,  °F 

ATSUb  Subcooling,  =  Tsat  -  Tb,  °F  (may  be  inlet,  outlet,  local,  or  average) 

/ 1  Viscosity,  lbm/ft  hr 

q  Density,  lbm/ft3 

a  Surface  tension,  lbf/ft 

r  Shear  stress,  lbf/ft2  [except  as  defined  in  Eqs.  (38)  and  (39)] 

4>m  Specific  mass  flow  rate  in  Eq.  (39),  kg/m2  sec 

SUBSCRIPTS 

avg  Average 

B  Fully  developed  nucleate  boiling 

Bi  Extension  of  the  fully  developed  nucleate  boiling  curve  to  incipience  point 

b  Bulk  fluid 

CHF  Critical  heat  flux  conditions 

conv  Pure  convection  conditions 

cr  Critical  heat  flux  conditions 

D  Diameter 

De  Equivalent  or  hydraulic  diameter 

e  Equivalent  or  hydraulic 
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f 

Film 

FC 

Forced  convection 

h 

Hydraulic 

htr 

Heater 

I 

Inner 

i,in 

Inlet 

i 

Liquid 

LB 

Subcooled  boiling  incipience 

m 

Mean 

O 

Outer 

o 

Outlet 

r 

Reduced 

sat 

Saturation 

ss 

Steady-state  conditions 

tr 

Transient  conditions 

ts 

Test  section 

V 

Vapor 

w 

Wall 
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